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CHAPTER 1 


CLASSIFICATION AND FORMATION OF PARTIAL 
DIFFERENTIAL EQUATIONS 


1.1. Partial Differential Equations (P.D.E.) 

An equation containing one or more partial derivative of an unknown 
function of two or more independent variables is known as a partiaL 
differential equations, we list the following examples of partial differential 
equations: \\X 






Order of A Partial Differential Equation 

The order of a partial differential equation is defined as order of the highest 
partial derivative occurring in the partial differential equation. 

In above examples of partial differential equations (l),(3), and (4) are of 
the first order, (5) and (6) is of the second order and (2) is of the third 
order. / 

Degree of A Partial Differential Equation 

The degree of a partial differential equation is the degree of the highest order 
derivative which occurs in it after the equation has been rationalized, i.e., 
made free from radicals and fractions so far as derivative are concerned. 

In above examples of partial differential equations (l),(2),(3), and (4), 
(6) are of first degree while (5) is of second degree. 

Notations: 

When we consider the case of two independent variable we usually assume 
them to be x, y and assume z to be the dependent variable. We adopt the 
following notations throughout the study of partial differential equations 




Sometimes the partial differentiations are -also denoted by making use of 

.pi ? '’1. du du d 2 u ’ 5 2 u j 

suffixes. Thus we write u r =—, « =—, u„=—r, =-and so on. 
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1*4. ^Classification of First Order Partial Differential Equations 

1.4.1. Linehr Equation 

A first order partial differential equation is said to be a linear equation if it is 
linear in p,q and z ,le., if it is of the form 

P(x,y)p+Q(x,y)q = R(x,y)z+S(x,y) 

Examples: yp-xq = xyz+x and p + q = z+xy 

1.4.2. Semi-Linear Equation ; i* " 

A first order partial differential equation is said to be a semi-linear equation 
if it is linear in p and q and the coefficient of P and. 0 the functions of 
x and y only i.e. if it is of the form \ 

P{x,y) P +Q{x,y)q = R{x,y,z) 


Examples: e x p- yxq -xz and yp + xq = —— 


Quasi Linear Equation ,' nNx^, 

A first order partial differential equation is $aidto v be a'quasi-linear equation 
if it is linear in p and q , i.e. if it is of the 


Classification and F ormation of Partial 
Differentia! Equations 



P(x,y,z)p+Q(x,y,z)q = R(x,y,z). ^ Njw 

x 2 + z 2 1 p - xyq = z 3 x+y 2 ^nH'^jp^Fy 2 zp—xy 

Non-Linear Equation 

Partial differential equations of the form f(x,y,z,p,q)=0 which do not 
come under the above three types are said to be non-linear equations. 

* if. ‘ ' ■ f,. ■ ' ’ V . i .. .. 

Example: pq = z , p 2 +q 2 =l and x 2 p 2 +y 2 q 2 = z 2 are all non-linear 
partial differential equations. t i ! 

^ < ! '/]•■ . U*.V* ; : // • ' 

Formation of Partial Differential Equation 
By Elimination of Arbitrary Constants' 

Consider an equation F(x,y,z,a,b) = 0 ...(1) 

Where a and b denote arbitrary constants. Let z be regarded as function of 
two independent variables X and y . Differentiating (1) w. r. t x and y 
partially, we get , .r 

dF 3F dF dF . , 0 , 

—+ p—and —+ ?—= 0 ...(2) 

ox dz oy oz 

Eliminating two constants a and b from three equations of (1) and (2), we 
shall obtain an equation of the form f ;(x,y,z f p,q) = 0 which is partial 
differential equation of the first order. 

In a similar manner it can be shown that if there are more arbitrary constants 
than the number of independent variables, the above procedure of 
elimination will give rise to partial differential equations of higher order 
than the first. ! ~ " ' r'*. 




argra 




mil 

































Example: Form tjie partial differential equation by eliminating arbitrary the 
constants a &b from z=ax+by+ab. 

Solution: Given z=ax+by+ab ...(l) 

Differentiating equation (l) partially with respect to x and y we obtain 


f-‘-« ' -< 3 > 

Substituting p and q for a and b in equation (l) , \ye get the required 
partial differential equation as z = px+qy+pq 


r 


Three Situations May Arise Y \ 

Situation I: When the number of arbitrary constanf^mgss than the number x 
of independent variables, then the elimination of arbitary'bonstants usually 
gives rise to more than one partial d^ffer^tial^ii^te^o^Srder one. 

For Example: Consider z = ax + y...( 1) 


Where a is the only arbitrary constant andare two independent 
variables. \ x > 


Differentiating (l) partially w. r .t. V. 


f'lr* 


Differentiating (1) partially w. r. 




Eliminating a between (l) and (2) 


dz) ... 

z= w + " 

Since (3) does not contain arbitrary constant, so (3) is also partial 
differential under consideration. Thus, we get two partial differential 
equations (3) and (4). /' 

Situation II: When the number of arbitrary constants is equal to the number 
of independent variables, then the elimination of arbitrary constants shall 
give rise to a unique partial differential equation of order one. 

For Example: Eliminate a f and h from az+6 = a 2 x+y ...(1) 

Differentiating (1) partially w. r. t. x and y , we get 


= a 2 ...(2) a — =1 


Eliminating a from (2) and (3) , we have —.— = 1, which is the unique 

, dx dy 

partial differential equation of order one. 

Situation III: When the number of arbitrary constants is greater than the 
number of independent variables, then the elimination of arbitrary:constants 
leads to a partial differential equation of order usually greater than one. . 


ms 

i 






^Sar^^tu^I^asi^p LI.^ New 8588844789 

^ .— _ _ E-mail: info^dipsaradei^^ . - 

























Art ISO 9001 i'2008 Certified Institute 




dx 1 

d 2 z 


For Example: Eliminate a,b and cfrom 
z=ax+Sy-¥'cxy 

Differentiating (1) partially w. r. t x and y , we have 


1/4 T 

— = o + cx 
dy 

Differentiating (2) partially w. r. t. xand (3)w. r. t. y, we have 

: K ■ 


Differentiating (2) w.r. t y , we have 


Now, (2) and (3) 


dz) J dz) A ^ 

=> x — =ax + cxy and y — =by\cxy 
\dxj yfy) \> 

r dz\ (dA l 

x — + y — = ax + by + cxy+cxy 

I An I _I 


W'UkJ • ,, iM .,my US i„,(ii.,»i(6)) ...(7) 

Thus, we get three partial differential equation given by(4), (5) and (7) 
which are all of order two. 

1.5.2. By Elimination of Arbitrary Function 4> 

Let «|»(u,v)=0, where u and v are functions of x, y and z. ...(1) 

We treat z as dependent variable and x andy as independent variables so 
that 

dz dz dy . , dx . ... 

— -p, — -q, — = 0 and ^- = 0. ...(2) 

dx dy dx dy 

Differentiating (1) partially w.r. t. x, we get 


dm du dx + du dy du dz' d<j> (dv dx dv dy dv dz 

duydx dx dy dx dz dx) dvl^dx dx dy dx dz dx 

dd>(du du\. d<t>(dv .. dv\ n ; / 

— — +p — +— —+/>— =0 - - 

du\dx dz j Svlcbc dz) o i- '1: ; ;; 
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a: : !ip+2 

du _ W dz 
d£ m + du' 
dv l dx P dz, 


Similarly, differentiating (1) partially w.r.t 'y ', we get 




du = & 

<¥ f du du' 


Eliminating <j> with the help of (3) & (4), we get 


dv dv'l (dv dv 
lk + P dzJ _ {dy + q dz 
du dw'j du du 
dx P dz J [dy + q & 


du duYdv dv^j (du 

Pp+Qq=R 
Where 

* 

p_du dv du dv du dv du 

dy dz dz dy’ dz dx 

Thus we obtain a linear partial differential equation of first order and of first 
degree, i.e. in p and q . -rv- / 

Example: Form a partial differential .equation by eliminating the arbitrary 
function (j) from <j){x + y + z, x 2 + y 1 - z 2 ) = 0 what is the order of this partial 
differential equation? 

Solution: Given <f>{x+y+z, x 2 + y 2 -z 2 j-0 ...(1) 

Let u =x+y+z arid v=£ +y 2 -z 2 ...(2) 

Then (1) becomes <J>(m,v) = 0 ...(3) 

Differentiating (3)w. r. t x partially, we get 
d<f>(du du'\ d<j>(dv 

— —+ p— +— —+p— =0 • ...( 4 ) 

0«^3r dz J dz J 

from (2), 

du , du du . dv . dv dv 

dx dy dz dx dy dz 

from (4)and (5), -v' ' ' 
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j -2(x-pz) 

td£\ (l + />) 

U'J 

Again, differentiating (3)w.r.t, y partially, we get 

dd>(du d<j>(dv dv'j 

— — +q — +— — +q — =0 
du\dy dz J 3v\<5y dz) 

* (1+,)+2 (f) (> ’“ z?)=0 ’ by(5) 

(d£\ 

ydtj -2 (y-qz) 

(d£\ '(l + «) 


Eliminating (f) from (6) and (7), we obtain N \. 

{x-^pz)_ = {y-qz) V ''\\\\ 

(i+p) (i+ 9 ) 

or (1 + q){x-pz) = (\ + p)[y - qz )or (y+Ap^^q = x-y 

which is the desired partial differential equation offifst order. 

Example : Form the PDE by eliminating ure s TtTbitrary function from 
z=f(x+it)+g(x-it), where i=jA 

Solution : V . >>■ 

Given z=/(*+#)+g(*-if) >Vw ' ...(1) 

Differentiating equation (l) twice partially with respect to x and/, we get 
j- = f'{x+it)+g'(x^it) t? * 

CyC ■ ■ ■ vi 


g= r( , +i ,)**•(>-«). , y ■ -<*) 

Here, /' indicates derivative of / with respect to (x+it) and g' indicates 
derivative of g with respect to (x - it ). Also, we have 1 

?L=i[f'(x+U)rg'{x-itj\ 1 

^ = ~{r( x + i( ) + s"{x- it)) :(3 ) 

From Equation (2) and (3), we at once, find that 

d 2 z th ... 

^ » ■ , ~< 4 > 
which is the required PDE. .. 

A'iu/e: We observed that by eliminating arbitrary functions, we always : 

produce quasi-linear partial differential equations, only. However, we get 
both quasi-linear as well as non-linear partial differential equations when we 

ftliminatp arhitrarv p.nnstflnts t \ • nSfi - ■ 


eliminate arbitrary constants. 
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CHAPTER 2 ' 


LINEAR PARTIAL DIFFERENTIAL EQUATION 
OF ORDER ONE x 




Lagrange’s Method 

Theorem: The general solution of the quasi-linear partial differential 
equation 

Pp + Qq-R 

is = 0 , 


where 41 is an arbitrary function and u(x,y,z)=q 


independent solutions of — = — = — 
P Q R 


Where q &c 2 are arbitrary constants and afle^st oite-of u&v must contain 


Working Rule For Solving Pp+Qq = RB\ Lagrange’s Method 


Step 1: Put the given linear partial differehti^requation of the first order in 
the standard form 

Pp + Qq = R -(I) 

Step 2: Write down Lagrange’s auxiliary equations for (1) namely, 

dx_cty_dz_ /«\ 

P Q R ■ ' 

Step 3: Solve (2) by using the well known methods. Let u(x,y,z)=q and 
v[x,y,z) =c 2 be two linearly independent solutions of equation ( 2 ). 

Step 4: The general solution (or integral) of (1) can be written in one of the 
following three equivalent forms: 

<)>(«, v) = 0 orw = <|>(v) or v=<)>(«) . 

Four Types of Problems Based on Pp+Qq = R 
Type-I 

dx dy _dz , , 

~P~~Q~~R ’ 

Suppose that one of the variables is either absent or cancels out from any 
two fractions of given equation (1). Then an integral can be obtained by 
usual methods. The same method can be repeated with another set of two 
fractions of given equation ( 1 ). n t r; : l 


— 
































' 

: 


M Example: Solve* * -— p + xzq = y z . 
x 


f 2 

y z o 

Solution: Given -— p+xzq = y 

X ■" 


The Lagrange’s auxiliary equations for (1) are 


dx dy _ dz 
(y 2 z/x) ** y 2 

Taking the first two fractions of (2), we have 
x 2 dx = y 2 dy or x 2 dx - y 2 dy = 0, V. 

integrating(3),( 1 /3)x 3 -(l/3)y 3 =(l/3)c, orx 3 -y 3 = 

V'Nw X 

Next, taking the first and the last fractions of (2) we 
xy 2 dx = y 2 zdz or xdx-zdz = Q 


xy 2 dx-y 2 zdz or xdx-zdz = 0 . f ... '.. ■/•••(5) 

.Integrating(5), (l/ 2 )x 2 -{l/ 2 )z 2 =)l'/ 2 j^ % oi^^W^ . ...( 6 ) . y: 

From (4) and ( 6 ), the required general integral i\ y\ 

/ » » 2 2\ L-«s. ■ (tv- 

<|>(x -y ,x -z j = 0 , <|> being an arbitrary frmcttoq^ 

Type-11 \i \T" : , 

-c>or~ /iv 

Suppose that one of integral of equation (1) is known by using Rule-1 and •' 
also suppose that another integral cannot be obtained by using Rule-1. Then . 

one integral known to us is used to find another integral as shown in /£jtX.-V - - 

following solved example. y-L; ■ 

Note that in the second integral, the constant of integration of first integral : L ’• 

should be removed later. 


Example : Solve p +3 q = 5z+tan [y - 3x) 

Solution: Given /7+3^ = 5z+tan(y-3x) 

The Lagrange’s subsidiary equations of (1) are 

dx dy _ dz 
1 3 5z+tan(y-3x) 

Taking the first two fractions, dy-3dx = 0 

Integrating (3), y - 3x = q 

Using (4), from (2) we get — = ——— 

1 5z + tanc ( 

Integrating (5), x-(l/5)log(5z+tanCj) = (l/5)c 2 
or 5x - log [5z+tan (y - 3x)] = c 2 , using (4) 


5x-log[5z +tan (y - 3x)] = <|>(y- 3 x), where <j) is an arbitrary function. 
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Type^III J 

. dx _ dy _dz, 

T~~q~~r 


Let I\,Q\ and R\ be functions of x,y and z . Then by a well known 
principle of algebra, each fraction in ( 1 ) will be equal to 

^c+Qidy+^dz 

P X P+Q X Q+R X R "V 

IfijP+Q 0 +^i? = O, then we know that the numerator of ( 2 ) is also zero. 
This gives P x dx+Q x dy+R x dz = 0 which can be integrated to give 
u x (x,y,z)=c x . This method may be repeated to 
u 2 (x,y,z) = c 2 .P X ,Q X ,R X are called multipliers. As a 

be constants also. Sometimes only one integral 
multipliers. In such cases second integral should b 
1 or Rule-2. 

Example : Solve j (b-c)l a}yzp + {(c-a)/Z>} zx<f 

■ ^ 4 \ 

Solution : Given -c)/«J yzp + |(c- a)/b}%xq 

i 

The Lagrange’s auxiliary equations of (1) a 

adx _ bdy _ cdz 
( b-c)yz (< c-a)zx (a-b)xy ^ 

Choosing x,y,z as multipliers, each fr^ion^ 

axdx + bydy+czdz _ axdx -f -bydy + czdz 

xyz\[b-c)+{c-a)+(a-bj\ 0 

axdx+bydy+czdz = 0 or laxdx + 2bydy+2czdz = 0 

Integrating, ax 2 +by 2 +cz 2 =q 


Again, choosing ax,by 9 cz as multipliers, each fraction of (2) 

a xdx + b ydy+c zdz 
0 


a 2 xdx + b 2 ydy + c 2 zdz 


xyz[a(b-c) + b(c-a) + c(d-b)] 

Integrating, a 2 x 2 + b 2 y 2 + c 2 z 2 = c 2 
From (3) and (4) the required general solution is given by 
<j>(ax 2 +by 2 +cz 2 ,a 2 x 2 +b 2 y 2 +c 2 z 2 j = 0, where <j> is an arbitrary function. 

Type-IV 

dx dy _dz 


Let and be functions of x~y and z, Then by a well known 
principle of algebra, each fraction of (1) will be equal to 

P x dx+Q x dy+R x dz 

P x PyQ x Q+R x R .a a;;-;, 
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Suppose" the numerator of (2) is exact differential of the denominator of (2). 
Then (2) can be combined with a suitable fraction in (1) to give an integral. 
However, in some problems, another set of multipliers P 2 ,Q 2 and R 2 are so 
chosen that the fraction 

P 2 dx + Q 2 dy + R 2 dz 

P 2 P+QiQ + R 2 R V ^ 




* 

■. ' ■ 


is such that its numerator is exact differential of denominator. Fractions (2) 
and (3) are then combined to given an integral. This method may be repeated 
in some problems to get another integral. Sometimes only integral is 
possible by using this Rule-4. In such cases second integral should be 
obtained by using Rule-1 or Rule-2 or Rule-3. 

Example: Solve (y+z)p+(z+x)q=x+y , ' 


9® v:--' m 

•V>3 4 - 7 >< v 


Solution : Here the Lagrange’s auxiliary equations are,. 


y + z z + x x + y 

V 

Choosing 1,—1,0 as multipliers, each fraction ot(. 

\ 

dx-dy _ d(x-y) ^ \ 

(y + z)-(z + x) -(x-y) 


Again, choosing 0,1,—1 as multipliers, each-fraction of (1) 

\ ' 

dy-dz d(y-z) _ V~~ 

(z+x)-(x+y) -(y-z) \^\i~ 

Finally, choosing 1,1,1 as multipliers, each'fraction of (1) 

_ dx+dy + dz _d(x + y + z) 

(y + z) + (z+x) + (x+y) 2(x + y + z) 


(2), (3) and (4) 


d(x-y) _d(y~z) _ d(x + y + z) 
-(x-y) -(y-z) 2 (x + y + z) 


Taking the first two fractions of (5), we have 
d(x-y) _d(y-z) 


. 


•. ft?-.' 


...(5) 


Integrating, log(x-y) = log(y-z) + logq 
or log{(x-y)/(y-z)} = logc, or {x-y)l{y-z) = c x 

Taking the first and the third fractions of (5), we have 

2 d(x-y) | d(x + y + z) =Q 
(x-y) x+y + z 

Integrating, 21og(x-y)+log(x+y+z) = logc 2 

Or (x-y) 2 (x+y+z) = c 2 

From (6) and (7), the required general solution is 

(|>{(x-y) 2 (x+y+z), (x-y)/(y-z)] = 0, <|) being an arbitrary function. 


-( 6 ) m 


•••« fa; 
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Cauchy Problem For First Order Quasi-Lineal* PDE or Integral 
Surfaces Passing Through A Given Curve 

In the last section we obtained general integral Pp + Qq = R . We shall now 
present two methods of using such a general solution for getting in the 
integral surface which passes through a given curve. 


be the given equation. Let its auxiliary equations give the following two 
independent solutions. 

u(x,y,z) = q and v(x,y,z)=c 2 ...(2) 

Suppose we wish to obtain the integral surface whichVas^esNhrough the 
curve whose equation in parametric form is given by N nT s -' 


Where t is a parameter. Then (2) may be expressed a: 




We eliminate single parameter t from the equal 
involving q and c 2 . k " s -v 


1 ) and get a relation lip- 


pm--! 


Finally, we replace q and c 2 with heto-of^^and obtain the required 
integral surface. s j 

Example : Find the integral surface oft^inearpartial differential equations 
x[y 2 +z\p-y[x 1 +z\q =(x 2 -y 2 \z wmciT‘“tontains the straight line 


Solution: Given xly +z 


Lagrange’s auxiliary equations of (1) are 


Solving we get 


Taking t as parameter, the given equation of the straight line x+y = 0, 
can be put in parametric form. 


Using (4), (3) may be re-written as 


These give 2(-Cj)-2=c 2 or 2cj+c 2 +2=0. ...| 

Putting values of q and q from (3) in (5), the desired integral surface is 
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Existence and Uniqueness of Integral Surface Passing Through a 

Given Curve ■ 

- 

Given the Partial differential equation Pp + Qq = R with the initial curve 

r:*b(0.tf>(0’ r o(0 —0) 

„ . O . .. T , ^(Wo. Z o)^0(Wo^o) 

Unique Solution: If-- * - -7 - 


then (1) has unique solution 


No-Solution: If 


solution 


Infinite Solution: If 


(1) has infinite solutions. \ >s \ ^ 

Example'. Consider the Cauchy’s problem < u~'-u y = 2 check whether 
equation has unique solution passing thro.ugh' tl}e'Gurve (25,5,2s) 

Solution: GivenP(x 0 (s),y 0 (s),z 0 (s)) =\ N 2(x 0 (^Vo ( s )> z o (*)) = ~ l and 


Now we define A 


1 + 2 = 3 * 0 (unique solution) 


Example: u x -u y =u passes through (s,s, sin s) then number of solution of 
this partial differential equation 0 < s < 1 

Solution: The given partial differential u x -u y =u ...(1) 

Compare (l)with Pp + Qq-R 
P = l,Q = l,R = u,x 0 =s,y 0 =s,u 0 =sins 


coss 


sins 


cos 5 


Hence (1) has no solution, equation (1) has zero number of solutions, 
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2.4 v> Surfai:& Offlogohal to a Given System of Surfaces 


represent a system of surfaces, where C is parameter. Suppose we wish to 
obtain a system of surfaces which cut each of (1) at right angles. Then the 
direction ratios of the normal at the point (x,y,z) to (1) which passes 

through that point are 

* 1 dx 8y dz 

Let the surface z = $(x,y) ...(2) 

cut each surface of (1) at right angles. Then the normal at (x,y,z) to (2) has 

direction ratios—,— , -Ii.e. p,q,~ 1 . Since normalsVt (x,y^-'to (1) and 
dx dy 

(2) are at right angles, we have 


which is of the form Pp + Qq = R. t 

Conversely, we easily verify that any solution of&) is orthogonal surface of 


Example : Find the surface which intersectKjhe surfaces of the system 
z(x+y)=c(bz+\) orthogonally and which passes through the circle 

x 2 +y 2 =l,z = l. 

Solution'. The given system of surfaces is'g'ivemby 


mmm. 


The required orthogonal surface is solution of 


Lagrange’s auxiliary equations for (2) are 




2(3z + 1) z(3z + 1) x + y 

Taking the first two fractions of (3), we get 

dx-dy = 0 So that x-y = C x 


Choosing x,y,-z(3z+l) as multipliers, each fraction of (3) 


xdx+ydy - z(3z+1 )dz 


xdx + ydy- lz L dz - zdz =0 
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Integrating (l / 2) jc 2 + (1 / 2) y 2 - 3 x 


imam 


Or x 2 +y - 2r - z = C 2 

Hence any surface which is orthogonal to (1) has equation of the form 
x 2 + y 2 -2z 3 - z 2 = y) <t> being an arbitrary function 


we must choose <j>(x-y) = -2. Thus, the required particular surface is 
x 2 + y 2 — 2z 2 — z 2 = —2 

The Linear Partial Differential Equation With n Independent 
Variables 

Let Xj, x>>v ..x n be the n independent variables and let p[~J—^PP 2 = ~~ > 

Q z \ ^ 

,.p n = -, where z is the dependent variable . CqnsiderJhegeneral linear 

partial differential equation with n ^independent vafiables/' v - . 

P l p l +P 2 p 2 +...+P n p n =R, -0) 

Where P\, P 2 ,..., P n are function ^v^3h-*2v--> x „. Let 

Mj =q, =C 2 ,.., u n =c n be any n independent v int'egrals of the auxiliary 

equations. 


Then the general solution of (1) is given , « 2 ,..., «„) = 0 

Note that the above procedure is generalization of Lagrange’s method. 
Example : Solve +x 3 x t p 2 +X\hPi +x ] x 2 X 3 =0 

Solution : Re-writing the given equation in standard form 
*2%Pl +^ X \Pl + XplPi =-X l X 1 Xy 
The auxiliary equations of (2) are 


x 2 x 3 x 3 x, x,x 2 x ( x 2 x 3 

Taking the first and the fourth fractions of (3), we get 

X\dx\ +dz =0sothat x 2 + 2z = c x 

Taking 1 st and 2 nd fractions of (3), we get 

Xydx^ —Xrjdxr^ so that x 2 —x 2 2 = c 2 

Finally, 2 nd and 3 rd fraction of (3) give 

Xjdxj = Xidxi, so that x 2 2 - x 3 2 = c 3 

Hence the required general integral is 

<t>(xj 2 +2z, x, 2 - x 2 2 , x 2 2 -x 3 2 j = 0,(J) being an arbitr 
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CHAPTER-3 

NON-LINEAR PARTIAL DIFFERENTIAL 
EQUATIONS 

3.1. Types of Solutions 

A solution or integral of a differential equation is a relation between the 
variables, by means of which and the derivatives obtaihed there from the 

V. \ \ 

equation is satisfied. Let us now discuss various classes .of integrals of a 
partial differential equation of order one. i. 


3.1.1. Complete Integral (C.I.) or Complete Solution,(C. S.) - ' 


Let us consider a relation (f> ( x, y, z, a, b) = 0 


\ w 


In which x,y,z are variables such that , z -is dependent on x and y 
Differentiating (1)partially w. r. t x and.y- 


respectively, we obtain 


(F F F 

the matrix M= ° ” ^ 

/b F xb F yb 


is two. 


d</> d<f> . d<f> d<j> 

—+ —n = 0 and—+—<7 = 0, 
d.x dz dy dz 


Since there are two arbitrary constants (namely a and b ) connected by the 
above three equations, these can be eliminated and there will appear a 
relation ofthe from f(x,y,z,p,q) = 0 ...(3) 

Which is a partial differential equation of order one. 

Suppose, that (l)has been derived from (3), by using some method; then 
the integral (1) , which has as many arbitrary constants as there are 
independent variables, is called the complete integral of (3). 


A two-parameter family of solutions 

z = F(x,y,a,b) ...(2) 

is called a ‘complete integral’ of (1) if in the region considered, the rank of 


V --X.ir 

& 


' 
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3.1.2. General Integral (G.I.) or General Solution (G.S.) 

' (d?u) du. . 

Assume that in k —r - — ip, one of the constants is a function of the 
\ox ) dt 

other, say b- F(a),then (1) becomes <j)[x,y,z,a,F[a)^ = 0. ...(5) 


Now (5) represents one of the families of surfaces given by the system (1). 
As before, the equation of the envelope of the family of surfaces given by 
(5) must also satisfy (3). Again the equation so obtained will be distinct 
from that of the envelope of the surfaces, and it is not a particular integral. It 
is known as the general integral and is obtained by eliminating a between 


<j>(x,y,z,a,b) = 0, |^ = 0 and |^ = 0. -(4) 

oa ob 

The relation between x,y and z so obtained is called the singular integral. In 
general, it is distinct from the complete integral. However, in exceptional 
cases it may be contained in the complete integral, that is, singular integral 
may be obtained by giving particular values to the constants in the complete 
integral. Since other relations may appear in the process of getting the 
singular integral, it is necessary to test that the equation of singular integral 
satisfies the given differential equation. 

3.2. Method of Getting Singular Integral Directly From The Partial 
Differential Equation of First Order 

Let the given partial differential equation be / (x, y, z, p, q ) = 0 .. .(1) 

whose complete integral is of the form §(x,y\z,a,b) = 0 ...(2) 


where 'a and 'b 1 are arbitrary constants. 

The singular integral of (1) is obtained by eliminating 'a 1 and 'o'between 
equation (2) 


and ^ = 0 


</>lx,y,z,a,F(a)) = 0. and = , ...(6) 

da 

Since other relations may appear in the process of 'getting ; the singular 
integral, it is necessary to test that the equation of general integral'satisfies 
the given differential equation. \ “ ' 

Particular Integral: A particular integral of (3). is obtained by giving 

particular values to a and b in (1) which is the complefeintegral of (3) 

\ \ \ v 

3.1.3. Singular Integral (S. I.) or Singular Solution (S.S.) 

We known that the locus of all the points who§e co-drdinates along with the 
values of p and q satisfy (3) , represent 'the clqubly infinite system of 
surfaces given by(l). The system is doubly infinite, since there are two 
constants a and b and each of these can take anjnfmite number of values. 
Since the envelope of all the surfaces given by ('l)Ts touched at each of its 
points by some one of these surfaces, file coordinates of any point on the 
envelope along with the values of p and ^belonging to the envelope at that 

point must also satisfy (3). Hence we conclude that the equation of the 
envelope is a solution of (3). The envelope of the surface given by (2) is 
obtained by eliminating a and b between the equations 
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The values of z, /?, q derived from (2) when substituted in (1) will reduce it 
into an identity and the substitution of the values of p and q (but not of z ) 
will in general render (1) equivalent to the integral equation. By using this 
substitution p and q are replaced by functions of x,y,z,a and ft in (1). It 
follows that the singular integral is given by (1) and the equations obtained 
on differentiating (1) partially w. r. t. ' a ' and 'b ', namely the equations 


*0, (5) & (6) hold if 


showing that there exists a functional relation betvyrierr p and q which does 
not contain ’ a 1 and 1 b' explicitly. Let this functional relation'b© 

ty(p,q) = 0 . V \\ \ ' ...(7) 

If both the constants ’a’and 'ft'occur in p and q (which does not always 
happen), then (7) shows that one of them is a ftinqtion df the other and the 
equations using them give general integral whicffismot now required. 

Equation (5) arid (6) are also true if ', x > 


Elimination of p and q from (1), (7) and (8) will yield a relation between 
x,y,z free from 'a'and 'ft'. If this relation satisfies the given differential 
equation (1), it must be the singular integral. 

Example: Consider 

f(x,y,z,p i q) = z-px-qy-p 2 -q 2 =0 ...(15) 

The two-parameter family of planes 

z = F(x,y,a,b) = ax + by + a 2 +b 2 , ...(16) 

is a complete integral of (15), since the matrix 




y + 2bO 1 


Let us now take b = J(l-a 


ax + 


On eliminating a , we get 
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This is a particular solution of the given partial differential equation. 
If we take b = a, then z = ax+ay+2a 2 


— = 0 =>x+ y = -4a . 
da 

On eliminating a, the envelope is 
8 z = -(x + y) 2 . 

This is another particular solution of the given p.d.e. 

Now from equation(l6), the conditions F a = 0 and F b = 0 become 




F a = {x +2a) = 0 X>J\ 7 ) 

F b =(y + 2b) = 0, 

Respectively. On eliminating a and b between equation'^), (17) and 
(18) we obtain the singular integral'.. \ \ % 




Which is a paraboloid of revolution. \ 

Note: The singular integral can also be obtained 'directly by eliminating p 
and ? between (15) and 

f P =-x-2p = 0, V ...(20) 


f=-y-2q = 0. 


3.3. Compatible System of First-Order Equations 

Consider the first order partial differential equations 

f{x,y,z,p,q) = 0 ; ^ 

and g(x,y,z,p,q) = 0 ...(2) 

Equation (1) and (2) are known as compatible if they have atleast one 
common solution. 

To find condition for (l)and (2) to be compatible. 

Let J - Jacobian of / and g = * 0. ...(3) 

d{P,q) 

Then (1) and (2) can be solved to obtain the explicit expressions for p and q 
given by 

p = <f>(x,y,z) and q = y/(x,y,z). ...(4) 

The condition that the pair of equations (1) and (2) should be compatible 
reduces then to the condition that the system of equations (4) should be 
completely integrable, i.e., that the equation 
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dz = pdx+qcty or </> dx+y/dy-dz=0, using (4) ■. ! ...(5) 

should be integrable (5) is integrable if 

V oz ) \ dz ) \dy dx) 

Which is equivalent to + <j>^~ = — + (//— ...(6) 

dx r dz dy T dz v 

Substituting from equations (4) in (1) and differentiating w. r. t. 'x'and 'z' 
respectively, we get 

df df d(j> df dy/ \ 

—+——- + ——— = 0 \ . \. , ...(7) 

dx dp dx dq dx V X. \ 





md o 

dz dp dz dq dz 


From (7) and (8), + f f...(9) 

ox dz op v ox oz ) ~xzq\.QX dz J 

\ • ; 

Similarly (2) yields \ ' \ 


^4;b : - 

dx dz dp dx dz) dq\ dx '- ""dz —, 


Solving (9)and(10),^ + ^4 ...(11) 

dx dz J ( d(xj)J \ di(z,p) J 

Again, substituting from equations (4) in (1) and differentiating w. r. t. 'y\ 
and 'z' and proceeding as before, we obtain 

f-HS-SSl 

Substituting from equations (ll)and (12) in (l)and replacing <fi,y/by p,q 
respectively, we obtain 

\%L*LMsA or [r ,g] =0> ...(is) 

J[d(x,p) d(z,p) J J { d(y,q) d(z,q) J 

Whe re [/, 2 1 ...04) 

d(x,p) S (z,p) d(y,q) d(z ,q) 

A Particular Case 

Theorem: To show that first order partial differential equations p ^ P(x,y ) 

dP 

and q - Q(x,y) are compatible if and only if — =—. 

dy dx 

dz dz 

Proof: Given— =p-P(x,y) and— = q-Q(x,y) ...(1) 

dx dy 

f dz ^ ( dz ^ 

Since dz = \ — <ix+ — dy-pdx + qdy, ...(2) 

UJ l dy) 
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It follows that the given partial differential equations (1) are compatible if 
and only if the single differential equation is integrable. 

dz-Pdx+Qdy ...(3) <m 

m 

Since Pand Q are functions of two variables x and y, hence Pdx+Qdy is 

an exact differential if and only if — =—.Therefore (3)is integrable if 

dydx 

and only if — = —. I 

dy dx 

Ifl 

dP dO v,;“- 

Remark 1. If— = —,then the system of two given partial differential 

dy dx \ . 

equations (1) is compatible and hence these will possess a common solution. 

V\ 0 

dP dO \ 

Remark 2. If— * —,then the system of two given partial differential 

dy dx . \*\._ •" 

equations (1) is not compatible and hence these''equations possess no 
solution. V v ‘ ''O’N ' 

Example: Show that the equations xp=yq^M^z:^xp + yq) = 2xy are 
compatible and solve them ■■ ^ 

Solution. Let f(x,y,z,p,q) = xp-yq = 0 \ N ^ ...(1) 

V w * Vw If 

andg(x,y,z,p,q) = z(xp + yq)-2xy = Q \ ^ ...(2) 

W a/I 





d(f,g) _ dx dp _ p x _ 2 \> 
“ d(x,p) dg dg zp-2y xz 
dx dp 


d(/,g) _dz dp _ O x 
d(z,p) dg_ dg xp + yq ‘ xz 
dz dp 


= -x p- xyq, 


d(/,g) _ dy dq = -q -y = ._ 2 

d(y,q) dg_ dg_ zq-2x zy ^ 
dy dq 

v v\ 

d{f,g) dz dq 0 -y 

and—^—-= = - =xyp+yq 

d(z,p) dg_ dg xp+yq zy\ • 

dz dq 

rj . i d(f,g) , d(f,g) Af’g) . ^f’S) „ ^ 2 0 22 

= -xp(xp + yq) + yq(xp + yq) = -(xp -yq)(xp + yq) = 0, using (1) 

Hence (1) and (2) are compatible. 






>; ■ \: ■ , v' • v ■ 

v< >• • ■■ 

4 + >< - \i ; f- . 

yy. •: 

- • '.* * '■ i •’ ‘V V- •-• ••: • ‘ -» *: 

:./.v • 


/ - •• . . i 
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Solving (1) and (2) for p and q , p = — and 


Using (3)in dz = pdx+qdy,wehwe dz = \^-]dx+ -\dy or zdz = d(xy). 


z c 

Integrating, — = xy + — or z 2 = 2xy + c, where c is an arbitrary constant. 

3.4. Charpit’s Method (General Method of Solving Partial Differential 
Equations of Order One But of Any Degree) 

Let the given partial equation differential of first order and non-linear in p 
and q be f(x,y,z,p,q )=0 V - v ...(1) 

We know that dz-pdx + qdy. •••(2)'' 

The next step consists in finding another relation F^y^'pq) = 0 ...(3) 

such that when the values of p and q obtairie'dj^^qlving (1) and (3), are 
substituted in (2), it becomes integrable. The integration of (2) will give the 
complete integral of (1). * '* 

In order to obtain (3), differentiate partially(l) and.(3) with respect to x 
and y and get 


\ v ' 


dx dz dp dx dq dx 

dF dF dF dp dF dq 

dx dz dp dx dq dx 

dy dz dp dy dq dy 

dF dF dF dp dF dq 

and — +— q +-——+—— = 0. 
dy dz dp dy dq dy 


Eliminating — from (4) and (5), we get 
dx 

'V + # p+ Z£m.J§L + W + 5F^W =0 

K dx dz dq dx) dp dz P dq dx)dp 

[dxdp dx dp) dz dp dz dp) P [dq dp dqdp)dx~ 

dq 

Similarly, eliminating -Y from (6) and (7), we get 
dy 

'VW_dFdf\(cfdF_dF_df) J^dF_dF^]dp_ = Q (g) 

K dy dq dy dq) [dz dq dz dqj q + [dp dq dp dqjdy " A ' 

Since — = —— = the last term in (8) is the same as that in (9), except 

dx dxdy dy . * 

for a minus sign and hence they cancel on adding (8) and (9). 
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Therefore, adding (8) and (9) and rearranging the terms, we obtain 


...do) 


dx dzjcp dz Jdq \ dp dqjdz 


l l dq)dy 


This is a linear equation of the first order to obtain the desired function 
F .Integral of (10) is obtained by solving the auxiliary equations. 


§L Lj£) (§r 

8x J l dz) 


dx _ dy _dF ~.. 

Wdl~dL~^' j 


dq J fy 


Since any of the integrals of (11) will satisfy (10), an integral of (11) which 
involves p or q (or both) will serve along with the givetr-qquation to find p 


and q .In practice, however, we shall select the simplest-integral. 


Note: In what follows we shall use the following standard notations: 


<L =f <L = f QL-f 0L = f 

dx Jx ’ dy Jy ' dz dp Jp 'dq /?i 


Therefore, Charpit’s auxiliary equations (11) may-bere-written as 


dp _ dp _ dz _ dx j dy\ dF* 

L+pfz fy+vfz -pf P -qf q ~f P 

\ \ s '- 

3.4,1. Working Rule While Using Charpit’s Method^ . 


Step 1. Transfer all terms of the given equation to L.H.S. and denote the 


entire expression by /. 


df df 


—,—..etc .occurring in step 2 and put these in Charpit’s 

dx dy 


equations (1 l)or (12). 


Sept 4 . After simplifying the step 3 , select two proper fractions so that the 
resulting integral may come out to be the simplest relation involving at least 
one of p and q . 


Step 5. The simplest relation of step 4 is solved along with the given 
equation to determine p and q . Put these values of p and q in 

dz = pdx + qdy which on integration gives the complete integral of the 
given equation. 

The Singular and General integrals may be obtained in the usual manner. 

Remark: Sometimes Charpit’s equations give rise to p-a and q = b ,where 
a and b are constants. In such cases, putting p = a and q = b in the given 
equation will give the required complete integral. 

Example: Find a complete integral of px+qy = pq. 





Step 2. Write down the Charpit’s auxiliary-equations (11) or (12). 

Step 3. Using the value of / in step 1 write down the values of 


Solution: Here given equation is f(x,y,z,p,q) = px+qy -pq = 0 

Charpit’s auxiliary equations are 

dp _ dq _ dz 

( dfldx)+p(df/dz) ( df/dy)+q(df/dz ) -p{dfd dp)-q{df /dq) 




- V-V’V* 
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4> rfy 


dp _ dq 


-( 2 ) || ? 


• • 


-O' - ?) -p(x-q)-<i(y~p) p+p - 0 ?+?-0 fjg|§| 

Taking the last two fractions of (2), (l/p)rf/> = (l/ 

Integrating log p =log q+ log a or p = aq. ...(3) |§p|?i 

Substituting this value of p in (1), we have 

aqx + qy-aq 2 =0 or 0 } = «+^, as ...(4) 0^i0t 

.’.From(3)and(4), <7 = (ax + y)/tf and p = ax+y v \ \ ...(5) 

Putting these values of p and q in dz = pdx+qdy, we j^tN^O ? 

rfz = (ox + y )^+[(ax + y) / a] or arfz = (ax + y)(adx'+dy)^ X ' ; V ' 

or adz = (ax+y)d(ax+y) = udu , wl\ere«&syw+y\' 

integrating, az = (l/2)u 2 +6 = (l/2)(ax+y) 2 +'b. ' ; . 

. \ Vy :• . r-- 

Which is a complete integral, a and 6 being-arbitrary constants. 

\ \ *' 

Example: Find the complete integral of the partial differential equation e. 

(p 2 + q 2 }x = pz and reduce the solutioh\which\passes through the curve 

x =°- z2 ‘ 4y - @fe i&m 






\ ■ - 

Solution: Let /( x,y,q,p,q ) =(p 2 + q 2 jx-^pz = 0. 


- 0 ) 


Chapit’s auxiliary equations are 


dp _ _ Jz _dx _dy 9 X -.. - . 

fx + Pfz fy+Pfz -Pf p -tfi ~f P ~f t 


giving ^ by (1) or 2 pdp + 2qdp = 0. 

q (~pq) ■ 


Integrating, p 2 + q 2 = a , where a is an arbitrary constant. 

2 / \ __ 

Solving (l) and (2), p=— and q= — x^(z 2 -aV). 

, a 2 xdx aJ(z 2 -a 2 x 2 )dy z dz-a 2 xdx , 
••• * =-/*& + qdy = — + - J-t -or n ---- = «fy. 

Putting z 2 - a 2 x 2 = t so that l{zdz- a 2 x dx} = dt, we get 


1 y fA - 

—= dt = a dy or - xt 2 dt=ady 

2 ft U J 


Integrating, t 2 =ay + b or ^[z 2 -a 2 x 2 j =ay + b, as t-z 2 -a 2 x 2 

Or z 2 -a 1 x 2 = {ay+bf or z 2 = a 2 x 2 +(ay+6) 2 " 

Which is complete integral. 


. .v ,y/ ■ / 

...( 2 ) 

■' ■' ■ • 

. - . . • 








pit: 1 

?> • >• ’“-v 

i^vK - ' -' i: - 

'A.;..;'., : 
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Non-Linear Partial Differential Equation 


THe pafametnc equations of the given curve x = 0, z 2 = Ay are given by 
jc= 0, y = r, z = 2t ■ ...(5) 

Therefore the intersections of (1) and (2) are determined by 

At 2 = {at 2 + £>) or aV +2[ab-2)t 2 +b 2 =0 ...(6) 

Equation (6) has equal roots if its discriminant=0, i.e., if 

A(ab-2) 2 -Aa 2 b 2 =0 or a 2 b 2 = 1 so that b = — 

a 

Hence from (4), the appropriate one parameter sub-system is given by 


z 2 -a 2 x 2 + ay + — or a 4 (jc 2 + y 2 ) + a 2 (2y - z 2 ) +1M), 

which is a quadratic equation in parameter V. Therefore, v this has form ifex. i 
envelope surface ^ 

(2y-z 2 ) 2 -4(x 2 + y 2 ) = 0or (2y ...(7) 

The desired solution is given by the function z v de finedby equation (7). 

Example: Find a complete, singular -Nand\]general integrals of 

(p 2 +q 2 )y = qz. \ ^ 

\ Vs* 

Solution: Here the given equation is \ ■ 

f(x,y,z,p,q) = (p 2 +q 2 )y-qz = 0. \ ...(1) 



■'■ ; r >: : ^‘Y? : :V■ 


Charpit’s auxiliary equations are 

dp _ dq dz _ dx _ dy 

f X + Pfz fy+Pfz -Pfp-qfl -fp -/« 




dp _dq 


dx dy 


j ~ - - 

-pq p -2 p y+qz-2q y -2 py -2 qy+z 
Taking the first two fractions, we get 2pdp + 2qdp =0 so that 


p 2 +q 2 =a 2 


Using (3), (1) gives a 2 y = qz or q -——. 

z 

Putting this value of q is (3), we get 

P = J/ “V) = ] f 2 -[^r-] 

Now putting these values of p and q in dz = pdx + qdy, we have 
. a [j 2 TTT, , a 2 ydy zdz-a 2 ydy 


dz = “J(z 2 -a 2 y 2 )dx + ^dy or f 


" 2 ^-«v) 

Integrating, (z 2 -a 2 y 2 ) 2 =ax + b or z 2 - a 2 y 2 =[ax + bf , 

Which is a required complete integral, a,b being arbitrary constants. 


: x . v ; 

'! T-; -. V *, ■ •• s " ,• . ; 

: fe'g&ifii.iSZt i -■ J - ’>'V •'' \ 

• ' 
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Singular integral. Differentiating (4) partially w. r. t. a and b , we have 


Q = 2ay 2 +2(ax+h)x ...(5) 

andO = 2(ax + Z>). ...(6) 

Eliminating a and b between (4), (5) and (6), we get z = 0 which clearly 
satisfies (1) and hence it is the singular integral. 

General Integral. Replacing b by <j>(a) in (4), we get 

z 2 - a 2 y 2 = [ax + <|>(a)] 2 . ...(7) 

Differentiating (7) partially w.r.t. 

a, -2 ay 2 =2[ax + cf)(a)J.[x + <j»'(a)J. ...(8) 

General integral is obtained by eliminating a from (7) and '(8) . % 


by (1) 


Putting q=a in (1), we get p= 


2 x(z-ay) 
(x 2 -a) 


Putting values p and q in dz = pdx+qdy, we get 

2x(z-ay) dz-ady 2 xdx 

dz = —~—;—-dx + ady or -— = —-. 

x -a z-ay x -a 




Example: Find a complete and isingdlhK.. integrals of 

2xz-px 2 -2qxy + pq = 0. \ \ \ 

Solution: Here given equation is \ 

' ' 1 

f{ x >y,z,p,q)=2xz-px 2 -2qxy+pq-0 \‘...(1) 

S \ \ 

Charpit’s auxiliary equations are 

dp dq _ dz _ dx _ ‘dy 

dx dz dy dz dp dq dp 'dq 

dp dq dx dy dz , 

or -—-— = -^- = — -= —-— = — -, by (1) 

2z-2qy 0 x -q 2xy-p px + 2xyq-2pq 

The second fraction gives dq= 0 so that q = a 


; 


' . 

itv-u 

:&v- : . 

•Cv- J • 

.#;A. : . v rr : - - . 

« ■ 

%Wi- • : : 

; •• ' 

i'.V . - 
• • 


Integrating, log(z-qy) = log(x 2 -aWlogfe ' 

or z-ay = b(x 2 -a) or z = ay + b(x 2 -a), ...(2) 

which is the complete integral, a and b being arbitrary constants. - •• 

Differentiating (2) partially with respect to a and b , we get 

0=y-6and 0=x 2 -a ...(3) 

Solving (3) for a and b, a=x 2 and b = y. ...(4) 

Substituting the values of a and'6, given by (4) in (2), we get z = x 2 y, 1 ||‘ ,y * *“ - 
which is the required singular integral. 
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Example: Find the complete integral, general integral and the singular 
integral of 2(z+xp+yq) = yp 2 

Solution: Given equation is f(x,y,z,p,q)=2(z+px+qy)-yp 2 =0 ...(1) 

Charpit’s auxiliary equation are 

dp dq _ dz _ dp _ dp 

fx+pf, f y +( ifz -pfp-qf, ~f P -f t 


2p + 2p 2q-p 2 + 2q -p[2x-2yp)-qx2y -(2x-2 yp) -2 y 

by (1) 

dp dy , dp —2dy 

Taking the first and the last fractions, — = or X-= 

4 P -2 y 

Integrating, log /> + 2 log y = log a or py 2 =a. , 'X, -(2) 

Solving (l)and (2) for p and q, /*=•“■-and r • 

>’ \ > >’ 2y 


, , , a , z ax a \ 

dz = pdx + qdy = —dx + -- + —+- rfy\-, 

y y y 2 \ 




; 



' - -Jj' 




Multiplying both sides by y and re-arranging, we-get 

\ \ 

♦ zrf,) - 4 

It J 2y \ \ 

or d[yz)-ad.^- -~^-y~ i dy= 0. 

7*) f a 1 \ 

Integrating, yz-a — + —- =b, 

yJ v 4 t J 

where a, b being arbitrary constants. 

General integral: Replacing b by <j)(a)in (3), we get 

\ f 2 v 
x a , v 

*"‘ww* w 


f \ f 

X Cl 

Differentiating (4)partially w. r. t. 'a',—- + —j = <|)'(a) ...(5) 

ItJ 12 y J ' ' 

Then the general integral is obtained by eliminating a from (4) and (5) 
Singular integral. Differentiating (3) partially w. r. t. 'a' and ’ft’, we get 

--+(^7=0 ...(6) 
w Ut 

0 = 1 - ...(7) 

Relation (7) is absurd and hence there is no singular solution of the given 
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3.4.2. Special Methods of Solutions Applicable To Certain Standard Forms: 

We now cottsiderequations in which p and q occur other than in the first 
degree, that is non-linear equations. We have already discussed the Charpit’s 
method. We now discuss four standard forms to which many equations can 
be reduced, and for which a complete integral can be obtained by inspection 
or by other shorter methods. 

Standard Form I: Only p and q present 

Under this standard form, we consider equations of the form 

f{p>q)= (i -0) 

Charpit’s auxiliary equations are 

dp dq _ dz _ dx _ dy _dF 

L+pL fy+pf, -pfp-qf, ~f P ~f q 0 \, 

X N *"x. 

y ^ \” 

dp dq m \\. \ \ 

giving -f = ~, by (1) w \ 

0 0 v X 

l \X v ‘ 

Taking the first ratio, dp = Q so that p—constants,.say ...(2) 

% \ X N 

Substituting in (1), we get f(a,q) = 0, giving^ Xcoftstant = 6, say, ...(3) 

X \ \ 

V \ \ 

Where b is such that f(a,b) = 0, X, — ...(4) 

Then, dz = pdx + qdy-adx + bdy, using.,(2) and (3J“ 

Integrating, z=ax+by+c , y ...(5) 

I, N. Xy 

Where c is an arbitrary constant, (5). together with (4) give the required 
solution. \s^. 

Now solving (4) for b , suppose we obtain b = F(a), say. 

Putting this value of b in (5), the complete integral of (1) is 

z = ax + yF[a) + c, ...(6) 

Which contains two arbitrary constants a and c which are equal to the 
number of independent variables, namely x and y. 

The singular integral of (1) is obtained by eliminating a and c between the 
complete integral (6) and the equations obtained by differentiating (6) 
partially w. r. t. a and c ; i.e., between 

z =ax + yF(a) + c, 0 -x + yF'[a) and 


Since the last equation in (7) is meaningless, we conclude that the equations 
of standard form I have no singular solution. 

In order to find the general integral of(l), we first take c = <f>{a) in (6), <j> 
being an arbitrary function and obtain z = ax + yF{a) + <f>(a). ...(8) 

Now, we differentiate (8) partially with respect to a and get 

^'0 = x + yF'(a) + 0'(a). ...(9) 

,f§£ Eliminating a between (8)and(9), we get the general solution of (1). 
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Non-Linear Partial Differential Equation 


• Remark: Sometimes change of variables can be employed to transform a 
given equation to standard form I. 

SOLVED EXAMPLE 

Example: Find a complete integral of the equation Jp+Jq = 1. 

Solution: The given PDE is of the form f{p,q) = 0. Therefore, let us assume 
the solution in the form z = ax+by+c 

Where Va+V<i> =1 or 6 = (l-\/aj 2 

Hence, the complete integral is found to be z = ax+(\-'Ja ) j L y + c. 

Example: Solve pq = k, where A'is a constant. 

•> V\ 

Solution: Given that pq-k. 

Since (1) is of the form f(p,q) = 0, its solution is z^'ax 1 + c , ...(2) 


Where ab-k 


or b = —\ on putting a for x p and b for q in (1) 

a \ \ ^— 


From (2) ,the complete integral is z = ax + {k hi^ + c; ...(3) 

\ \ 

Which contains two arbitrary constants a and.c . N'— 

For singular solution, differentiating (3) partialfy with respect to a and c , 

( k \ !'ST""" 

we get 0 = x - — y and 0 = 1. Biit OM.is' absurd. Hence there is no 
singular solution of (1). 

To find the general solution, put c = ^(a) in (3). Then, we get 


= ax + \— y + ^(a). 
\aj 


Differentiating (4) partially with respect to 'a' , we get 


O = *-(yJj' + 0'( a )- -(5) 

Eliminating a from (4) and (5), we get the required general solution. 
Example: Solve p 1 + q 2 =m 2 , where m is a constant. 

Solution: Given that p 1 +q 2 =m 2 , ...(1) 

Since (1) is of the form f(p,q)~ 0, its solution is 
z = ax + by + c, ...(2) 

Where p 2 +q 2 =m 2 , or h = (/n 2 -n 2 ) 2 , on putting a for p and b for q 

in (1). 

.-. Form(2), the complete integral is z = ax-\-y{m 2 -a 2 ) 2 +c, ...(3) 

Which contains two arbitrary constants a and c . 
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:». i , V. : : 


For singular solution, differentiating (3) partially with respect to a and c , 


we get 0 = x- - - — — r and 0 = 1. But 0 = 1 is absurd. Hence there is no 

( 2 2 \\ HH 

\m -ay 

singular solution of (1). 

To find the general solution, putc = <J>(a) in(3). Then, we get 
z = ax + y(m 2 -a 2 ) 2 +<|>(a). 

Differentiating (4) partially with respect to 'a', we get 


...(4) 


' 

: 

’t-- - • : ; 

. 


0 = a 


\ 


\ 

\ \ X 

\\V 

\ 


...(5) 


XX •‘V ; ■ 
! . 


(m 2 -a 2 ) 2 

Eliminating a from (4) and (5), we get the required general solution. 

X>.. “■ 

Standard Form II: Clairaut Equation. - \ 

A first order partial differential equation is said'to be.of -Glariaut form if it 
can be written in the form 


! 


z = px + qy + f(p,q). 

Let F(x,y,z,p,q)=px+qy+f(p,q)-z \ 

Charpit’s auxiliary equations are 
dp dq 


\ W 
\ 


NX 



\ ^ 


...( 1 ) 

...( 2 ) 




dx dy 


t dx J X & 

dp dq 
or — = —= 


r dF' 


dy 


, dF , 
+ «'&' 


dz 


/ aF ^ Tar-X 


K d Pj 


~<1 


dF 

\ d <lJ 

dx 


dF_ _dF 
dp dq 

dy 


0 0 I 



f a/ ] 


(df ) 


(df\ 

1 

*8 

1 

S 

1 

^3 

kPP-j 

-q 

\p9j 

-x- 

K d P) 

-y- 

A) 


,by(l) 


Then, first and second fractions => dp = 0 and dq- 0 => p = a and q=b. 

Substituting these values in (1), the complete integral is 
z = ax + by + f[a,b ) 

Remark 1: Observe that the complete integral of (1) is obtained by merely 
replacing p and q by a and b respectively. Singular and general integrals 
can be obtained by usual methods. 

Remark 2: Sometimes change of variables can be employed to transform a 
given equation to standard form II. 

Example: Find the complete integral of the equation 
z = px+qy + F+ p 2 +q 2 

Solution: The given PDE is in the Clairaut’s form Hence, its complete 
integral is 


z = ax+by+Jl+a 2 +b 2 . 


t 




• I ■ 

XX y 

I; 

-y” ^ 

X • 

i NX 

|: NX 

■ 


. 


ipN 


' 

- 


I 

' 


' ' ' \ :■ 

*•' v.-- • : V. • 
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Example: Find the complete and singular integrals of the 
Z=*px + qy*]og(pq) 



Solution: The complete integral is z = ax + by + log(«Z>) 

or z = ax J rby + \oga + \ogb, a, b being arbitrary constants ...1 

Differentiating (1) partially with respect to a and b , we get 

1 N 1 ^ 1 1 

0 = jc+ — and 0 = y + — so that a= —and 6 = — ...( 

\aj \b J x y 

Eliminating a and b from (1) and (2), the required singular integral is 


z = -1 -1 + log — or z = -2 - log (xy ). « "'X \ 

yXy) \\\ 

Example: Find the complete and singular, -integrals of " the 

z = px + qy-ljpq. V \X 


Solution: The complete integral is z =,ax + by - 2<fah. ...( 

S, s 

Differentiating (1) partially with respect to a and 'A , we,get 

0 = x —and 0 = y-(—]== so that'** = M"and y-.f^-. ...(2) 

2y/ab l 2 'Jab) \ V b 


Now, using (1) x-z-x-{ax + by-2-fab^ = \ 


' a ialJ + 2 ^’ 


using (2) 


Similarly, using (1) 



-z = y-(ax + by-2^/ab),= . --a. — -bj— +2\[ab 

' A V h V n v h 



From (3) and (4), f 

Which is singular integral as it satisfies the given equation. 

Standard Form III: Only p , q and z are Present 

Under this standard form we consider differential equation of the form 

f(p,q,z)=0 

Charpit’s auxiliary equations are 

dp dq dz dx dy 


d ±) +p m m + (v 

dx) V dzJ Idyl 


V) Jar) iK -JL 


\dp) \dqJ dp dq 


dp _ dq 

(dfY Tdf 

p — q — 
A dz ) l dz 


dz _ dx _ dy 

dp J \dq J bp dq 


, using (1) 
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ill 111 

Taking the first two ratios, — dp- - dq 

\p) uj 

Integrating, q = ap, a being an arbitrary constant. 
Now, dz = p dx+q dy= p dx + ap dy, using (2) 
or dz-p (dx+a dy)= pd[x + ap)- p du. 

Where u = x + ay. 


Now, (3) => p =— and so by (2) q = ap - a — 
du \du 


Which is an ordinary differential equation of first order: Solving (3), we get 
z as a function of u . Complete integral is then obtained, by replacing «by 
[x + ay). '■*“ 

Working Rule For Solving Equations of The Form 

f(p,q,z) = 0. \ "'X.JC' ' -(1) 

. \ A 

Step 1. Let u = x + ay, where a is an arbitrary.constant. ...(2) 

\ 'N C. 

dz ■ (dz ^ 

Step II. Replace p and qby —and a -— respectively in (l)and solve 

du , ' 

the resulting ordinary differential equafioft of first order by usual methods. 

V \ ^ 

Step III. Replace u by x +ay in the solution,obtained in step II. 

Remark 1: Sometimes change of variables can be employed to reduce a 
given equation in the standard form III. 

Remark 2: Singular and general integrals are obtained by well known 
methods. 

SOLVED EXAMPLE 

Example: Find the complete integral of p(\+q) = qz . 

Solution: Let us assume the solution in the form 

u=x+ay ,where a is an arbitrary constant. 

T , dz dz 

Then,p =—, q = a— 
du du 

Substituting these values in the given PDE, we get 

dz (, dz \ dz 

— 1 + a — \ = az —. 
du\ du J du 


That is,a— = az-l or a — — = du 
du az-l 

On integration, we find ]n{az-l) = u + c = x+ay+c 
Which is the required complete integral. 
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Non-Linear Partial Differential Equation 


Example: Find a complete integral of 9(p 2 z + q 2 ^-4. yijhi 

Solution: Given equation is 9 [p 2 z + ? 2 ) = 4, —(1) 

which is of the form f(p,q,z) = 0. Let u = x + ay, where a is an arbitrary 0 

constant. Now, replacing p and qby -^and a|-^ respectively in(1), i}/ 

du \du J 

we get 

r dz v jdz) 2 ] ■ (dz y_ 4 

9 z — +a — =4or — -~7 -pr- 

\du) \du) J W M / 9[z + a J 

or rfu = ±f|lx(z + a 2 ) 2 <fe,separating variables u and x 










is*-* yiy/'z.-' 


Integrating, u + b = ± - x 

v 2 J 




\ £ - 

or M + ^=d;'(z-i.q 2 


or (« + 6) 2 =(z + a 2 ) 3 or (x + ay + 6) 2 = (z + a 2 ) 3 , as .^^x + ay 

V . . \ \\ 
which is a complete integral containing two arbitrary constants a and b. 

Example: Find complete and singular integrals 6 f^T + z 3 j = 9 z*pq. 

._ N. 

Solution: The given equation is of the\Torm f(p,q,z) = 0. Let 

| fa 

u = x + ay, a being an arbitrary constant. ..Replacing p by — and q by 

''v, ^ ^ du 

(dz \ Vv^,. 

a — in the given equation, we have 
du) 


4(l + z 3 ) = 9z 4 a — or ±- —7 

U») (, + z 3)i 


dz = 2du 


v;,-: 


or ± — x2t dt — 2du, putting 1 +z -t so that 3z dz — 2t dt 

\ 1 J 

Integrating, ±4at-u + b or ±Ja (l + z 3 ) 2 = x + ay + b 

or a(l +z 3 ^ = (x +ay+ b) 2 , —(0 

which is a complete integral containing two arbitrary constants a and b. 

Singular Integral: Differentiating (1) partially w. r. t. a and b by turn, we 

get py' ""?• v-& ypl: 

\ + z % = 2y(x + ay + b) —(2) Spv yA X+AM 

and 0 = 2 (x + ay + b). —(3) ^py.yi'- 

Eliminating a and b from (1),(2) and (3), the singular integral is 
1 + z 3 = 0 . -(4) 

From (4) , /> = —= 0 and ? = —= 0 . Thus these values of p and 4 
together with 1 + z 3 =0 satisfy the given equation. Hence 1 + z 3 =0 is the 

required singular integral. . 
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Non-Linear Partial Differential Equation 


Standard form IV: Equation of The Form /, (x, p) = fi (y, <?) • i.e. a form 
in which zdoes not appear and the terms containing xand pare on one 
side and those containing y and q on the other side. 

Let F(x,y,z,p,q) = f l (x,p)-f 2 (y,q) = 0. -0) 

Then Charpit’s auxiliary equations are 


clFV (5F) (3F) (8F 

dx) \dz) [dy J + ^dz 

dp _ dq _ dz 

r IT"j/r_ (5T i 

dx 8y M dp J 


dz_ _ _dx_ _ _dy_ 

~ p {dp) 9 [ai) ~ dp dq 


, by (!) 


Taking the first and the fourth ratios, we have \ 

dp f -^-1 dx 0 or df x = 0 ., _ V % 

ydp) \dx) v . s 

Integrating, f x =a, a being an arbitrary constant. \ \ 

.-.(1) => f l {x,p) = f 2 (y,q) = a. 

Now, (2) =>f l (x,p) = a and f 2 ( y,q ) -(3) 

* \ V 

Form (3), on solving for p and q respectively, we get 

\ \ ^ 

p = F x (x,a), say and q = F 2 (y,a), say \v_ ...(4) 

Substituting these values in dz = p dx + q dy, we get 
dz = F x [x,a)dx + F 2 [y,a)dy. 

Integrating, z = |fj (x,a) dx + |F 2 (y,a)dy + b, 

Which is a complete integral containing two arbitrary constants a and b 

Remark.l: Sometimes change of variables can be employed to reduce a 
given equation in the standard 

form IV. 

Remark.2: Singular and general integral are obtained by well known 
methods. 

SOLVED EXAMPLE 

Example: Find the complete integral of the PDE: p 2 y (l + x 2 j = qx 2 
Solution: The given PDE is of separable type and can be rewritten as 
P 2 (l + x 2 ) q 

——- = — = a (say), where a is an arbitrary constant. 

u 


4ax 

p=-== ) q = ay 
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Substituting these values of p and q in 


dz = pdx+qdy, 

We get dz =-^=dx + aydy 
dl + x 2 


On integration, we obtain z = \fayl + x 2 +^y 2 +b 

which is the complete integral of the given PDE. 

Example: Find a complete integral of x(l + y)p = y{\ + x)q. 

Solution: Separating p and x from q and y , the given, equation Reduces to 


. 

M = M v\ 

( l+x ) ( 1 + t) ACx. 

Equating each side to an arbitrary Constant a, we have 


yu t i i A jv > 

= a -a so that p-a - \and 




Putting these values of p and q in dz - p dx + q'ty we get 




dy or dz = a\— +'k [dx + a — +1 dy. 

^ \y J 


o—- v -» , v —-v-'g xy + x + y) + b, 

which is a complete integral containing two arbitrary constants a and b. 

Curves & Surfaces in Three Dimensions 

A point in 3 -dimensions connected by relation /(. x,y , z) = 0 ...(l) 

Lies on surface => equation (1) is an equation of surface S 

* If a point with coordinates (x,y,z) lies on two surfaces ^ & s 2 given by 

f{x,y,z) = 0 &g(x,y,z) = 0 ...(2) 

Then the surfaces s { & s 2 intersect in curve C. 

Therefore locus of point, whose coordinates satisfy a pair of relations of type 
( 2 ) is a curve in space. 

* Any three equations of the form x = f x (f), y = f 2 (t) & z = / 3 (t), .. .(3) 

in which t is a continuous variable may be regarded as parameter equation 
of a curve. 

* If we take length of the curve as measured from some fixed point as 
parameter and P any point on 

the curve x = x(i), y = >’(•?), z = z(s) ...(4) 

We can find that direction cosines of the tangent to the curves (4) at point P 
are given as 
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dx dy dz 
ds ’ ds' ds 


* Assume that curve (4) lies on the surface s with equation F(x,y,z) = 0 
Therefore point P(x(s),y(s), z(s)) lies on surface s 

i.e. F(x(s),y(*),z(s)) = 0 ...( 6 ) it 


dF dx dF dy dF dz 

—.— +-- +-= 0 . 

dx ds dy ds dz ds ' ' 




From equation (5) & (7) it is clear that tangent to curve C at p is 
perpendicular to the line whose direction ratios are " \ 

'dF dF dF) H 

.dx'dy'dz) ; \>- v \ ,■ 


Obviously these are d.r’s of normal to the surface at X . " - 

V ^ \ 'N \ \ 

Example: Surface S:z = f(x,y) ' \\ ' 

=> F(x,y,z) = f(x,y)-z = 0 \ \\ 

■ " 

=> Fx = p, Fy = q, Fz = -1 \ V '"'s. N 


x N 

\ 

v>— 


Hence d.c’s of normal to, surface at (x, y, z) 


^p 2 +q 2 +\ yjp 2 +q 2 +1 ^p Z +q 2 + 1 J Ipy 

* Using (5) & (7) we can write the equation of tangent plane at point 
P(x,y,z) to the surface S { with equation F(x, y,z) = 0 is 

(*-x)f*(r-,)f t( z-,)£., ...(9) || 

Where x,y,z are any other point on the tangent plane. Similarly equation of 
tangent plane to the surface S 2 with equation G(x, y, z) = 0 at the point 
P(x, y, z) is W&U 

{X-x)^-+(Y-y)^-+(Z-z)~- = o ...(10) 

* The intersection of tangent planes given by (9) & (10) is the tangent at p 
to the curve C which is intersection of surface S i &S 1 . From (9) & (10) 
equation of tangent line is 

X-x _ Y-x z-x 

dF dG dFdG dFdG dFdG 'dFdG dF dF W&M 

dy dz dz dy dz dx dx dz dx dy dy dx '$$$$*■ 

=> direction ratios of tangent line are Sfeft! 


- 1 


d(F,G) d(F,G) d(F,G) 
d(y,z) ’ d(z,x) ’ d(x,y ) 


5 


mi?-.: 
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3.6. Cauchy's Method of Characteristics for Solving Non-Linear 


P.D.E. 





, 


Let non- linear P.D.E is: FI x,y, z,— I - 0 ...( 1 ) 

•x/;'•< V'XV'< 

The plane passing through point P(x 0 ,y 0 ,z 0 ) with its normal parallel to the 
direction n defined by direction ratios is uniquely specified by 

set of numbers D(x 0 ,y 0 ,z 0 ,p 0 ,q 0 ) -> called plane element of space). 

Conversely any such set of 5 real number defines a plane is 3-D. In 
particular a plane element ( x 0 ,y 0 ,z 0 ,p 0 ,q 0 ) whose components satisfy on in¬ 
equation V S \0 

F(x,y,z,p,q) = 0 •'••(?) . 

is called an integral element of equation( 2 ) at (x 0 ,/( H z 0 )'.x. ; . . . 

It is possible to solve (2) to obtain an expression \T“x 
q = G(x,y,z,p) -( 3 ) 

Keeping x 0 ,y 0 &z 0 fixed and varying p, we obtain a set of plane elements 
{x 0 ,V 0 ,z 0 ,p, ^( jc 6 »To> z o>p)} which depend pn'single parameter p. As p 
varies, we obtain a set of plane elements, all of which pass through point P 
and which therefore envelope a cone withVertexV, the cone so generated is 
called elementary cone of equation ( 2 )'at.poiht^_ 

Consider a surface 5 with equation z^ (x,y)'-~ ...(4) v 

If g(x,y), g x (x,y) & g y {x,y) are continuous in certain region R of xy 

plane, then tangent plane at each point of s determines a plane element of YYY, Y;y ; Y ■ ■ 

type A'' 

{*o Do ),g x (* 0 . v 0 )»£, ( x o > V 0 )} •••(5) 


called the tangent element of surface j at point {x 0 ,y 0 ,g(x 0 ,y 0 )} 

Theorem: A necessary & sufficient condition that a surface be an integral 
surface of a P.D.E. is that at each point its tangent element should touch the 
elementary cone of the equation. Y Y 

Proof: Let surface s be z = g(x,y) ...(l) 

A curve c with parameter equation ...(2) : 

lies on surface (1) if z(*) = g(x(r), y(f)) Vte/ (say) 

If P 0 is a point on the curve corresponding to value r 0 , then the direction -Y Y • ' - : ; 

ratios of the tangent line P 0 P x are {x'(r 0 ), y'(t 0 ), z'(t 0 )} . This direction will /'.-AY 

be 1 to thedirection (p 0 ,q 0 - 1) if z'(t 0 ) = p 0 x' 0 (t 0 )+q 0 y' 0 (t Q ) cW Y- 

So, any set {x(t)>.F(0’ z (0>P(')>f(0} -( 3 ) 

of 5 real number satisfying the condition 

•P P PPI 

z'(t) = p{t)x'(t) + q(t)y'{t) ...(4) 

Defines a strip at the point (x,y,z) of curve c. . 
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If such a strip is also an integral element of equation 

F(x,y,z,p,q) = 0 ...(5) 

We say that it is an integral strip of equation (5) 

Therefore set of from (3) in an integral strip of equation ( 2 ) provided they 
satisfy conditions (4) and conditions. 

F {x(t),y(t),z(t),p(t),q(t)) = ti ...(6) 

•Vfe/. 

If at each point the curve (2) touches a generator of elementary cone, we say 
that the corresponding strip in a characteristic strip the point 
(x + dx,y + dy,z + dz) lies in the tangent plane to the elementary cone at P if 

dz = pdx + qdy V\ \ ~~~ 

\ \ \ ^ 

where p, q satisfy (5) \ v ^ 

\ ^ V"n 

Differentiating (7) w. r. t. p we get v - \\ \ ' 

• 

0 = dx +—dy , \ s \ ...( 8 ) 


from (5) 


dF dFdq n 
dp dq dp 


Solving (7), (8) & (9) for mtio?''"of dy, dz to dx we have 

dx dy dz _ . . 

"' (10) 

=> Along a characteristic strip, x'(t),y'(t),z'(t) must be proportional to 
Fp, Fq & pFp + qFq respectively. 

If we choose parameter t in such a way that 

x'(t) = Fp] 

/«)-4 •<"' 

z'(t) = pFp + qFq ...(12) 

Along a characteristic strip p in a function of t so that p’(t) = + 

w dx dt dy dt 

_dp_dF_ + dp dF 
dx dp dy dq 


dp_dF_ + dq_dF_ 
dx dp dx dq 


Since 

dy dx 


Differentiating equation (5) w.r.t. * we get 
dF dF dF dp dF dq A 

—+—p + —-—+-- = 0 

dx dz dp dx dq dx 

so that on a characteristic strip ? 
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pV) = -( F x+pFz) ...(13) 

Similarly we may obtain 

q'(t) = -(F y+ qF z ) ...(14) 

Now we have a system of following O.D.E. to determine characteristic strip. 


x'{t) = Fp,y'{t) = F r z'(t) = pF p +qF q ,p'(t) = -(F x +pF z )} 

&q'(t) = -(F y+ qF z ) J "‘ (15) | 

these are called characteristic equations of D.E. (5) 

Theorem: Along every characteristic strip of equation Ffay,z,p,q) = 0 the 

V *\v v 

function F(x,y,z,p,q) is a constant. \'\. -1 

Proof: Along a characteristic strip we have \ >> n f 

j'F{x(t),y(t),z(t),p(t),q(t)) = F x X' +~F y y + F z z\+ F p p' + F q q' 

= F x F p + F y F q + F z ( P F p +qF q )-F p (F x + P F Z )- F^F y + qF z ) 


=> F(x,y,z,p,q) = k a constant, along the strip.'-. - 

, \ 1 & 

Theorem: If a characteristic strip contains atleast one integral element of 

F(x,y,z,p,q) = 0, it is an integral strip of.equation F(x,y,z,z x ,z y } = 0 


Proof: Let solution of PDE, F{x,y,z,z x , z y ) = 0 ...(l) |f- 

Passes through curve F whose equations are 

x = 6(v),y = <t>(v),z = Y(v) ...(2) 

then in the solution x = x(p 0 ,q 0 ,x 0 ,y 0 ,z 0 ,t 0 ,t) etc ...(3) 

of the characteristic equation (15) we may take 
x o = 9( v )>Fo =( K V )> z o = v( v ) 

as initial values of x,y,z. The corresponding initial values of p 0 , q 0 are 
determined by 

V'( v ) = A> e '( v Wof( v ) 

F{Q(v),ty(v),y(v),p 0 ,q 0 } = 0 

If we substitute the values of x 0 ,y 0 ,z 0 ,p 0 ,q 0 and the appropriate value of t 0 
in equation (3), we find that x,y, z, can be expressed in terms of two 
parameters t & v to give 

x = X(v,t),y = Y(v,t) & z = z(v,t) ...(4) 

Eliminating v, t from these three equation, we get a relation 
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v(x,y,z) = 0 

which is equation of integral surface of equation ( 1 ) through curve r. 
Example: Find solution of equation z = ^p 2 +q 2 ^+(p-x)(q-y) 

Which passes through the x - axis 

Solution: Here F(x,y,z,p,q) = z-^p 2 -^q 2 -(p-x)(q-y) = 0 ...(l) 

Initial values are: x 0 =v,y 0 = 0, z 0 = 0 ,p 0 = 0, q 0 = 2v, / 0 =0 characteristic 
equation are 


*** r? 

- = F p =-p-q + y 



dt q 
dz 


F a =-q-p + X 


= pFp + qF q = P(-P-q + y) + q (-#- M*) = ~p^q : 2p<? + py + qx ...(A) 

\ \ \ 

-s \ 

\ 

= -F x -pF z =-q + y-p \ \\ ...(5) 


-F v ~qF z =-p+x-q 


\ ‘X. 

\ N- 


From (2) & (5) => — = — => x = p'+x< \ 

dt dt \ \ v 

\ v "' 

at t = 0 , x 0 = v, p 0 = 0 => c, = v v ' s — 

Hence x = p+v...( 7) 

Similarly from (3) & ( 6 ) we have y-q-2v 

From (4) we have + — = +q + p-y+ p+q-x-p-q+y 

dtdtdt 

=> J t i-P~<l + x ) = P + <l-x 

=> p + q-x = ve~‘ 

Similarly ^(p+q-y) = -p-<]+.y => p + q-y = 2ve~ t 

from(7)&(8) p+q = x+y+v 
So equation (9) & (10) becomes 


y + v = ve~ l 


y=v(e~ t - 1) 

& x + v = 2ve~ l 


x = v\2e~* -1 

Therefore p = 2v(e~‘ - 

l) & 

q = v(e~ l +lj 

Therefore (4) becomes 



^ = 2v(e'-l)(- 2 v<f') 

+ V^£ 

r,+i )(- vc '') 


mrnmm 
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z =—v e +3v e +c 2 
2 


5 , ^ 

at t 0 =0, z 0 =0, => c 2 =--v +3v =+— 

(4) z = — v 2 e~ 2( + 3 v 2 e~ ( + — 

2 2 

From (11) we have 
x = v (2 e~ l -lj 

>' = v(e w -l) 

=> x - y - ve~ l 

j; = ve~*-v = x-'y- v =>v = 2 ^-x 


_/ x-y 

x-v = ve => e =—— 
2>>-x 


Therefore (12) becomes 

/ \2 5 fx-^) 2 (x-J>) 

z={2y - x) r^ +3 idf 


2 \?>{x-yf+6(x-y){2y-x) + (2y-x) 


= ^5x 2 + 5/ - 10xy +18y - 6 x 2 -12/ + 4/ + x 2 - 4xyj 
^ z = L[-3y 2 + 4yx] = ^y[4x-3y] 

Example: Find the characteristics of equation pq = z & determine the 
integral surface which passes through parabola x = 0,y 2 = z. 


Solution: Given P.D.E. F[x,y,z,p,q) =z-pq=0 

- 0 ) 

2 V 

Initial value is x 0 = 0 , t 0 = 0 , y 0 = v, z 0 = v , p 0 = -, q 0 = 2 v 

characteristic 

equations are 


S3n 

1 

II 

O 

II 

•SI'S 

...( 2 ) 

II 

II 

1 

...(3) 

^- = pF q +qF q =-2pq 

...(4) 
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Hence 


Solving (9) & (10) we get p 


From equation (4) 


Hence z = v e 


From equation (7), (8) & (11) Q we Have; x = 2v( e 1 -1 ) & y 


Therefore equation (12) becomes 

.. ( 4 p -*) 2 .,( 4 y +*) 2 


Example: Determine the characteristic strip of equation z -■ 
the integral surface passing through parabola x 2 +4 z = 0, y 


Solution: Given P.D.E. 


Initial values are, t ( 


Characteristic are 


. • 
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How the choose initial values: 


Given that y 0 = 0 


Let x 0 = v then x 2 + 4 z =0 gives, z 0 = ——, by 

differentiating x +4z = 0 , we get, p 0 =- — and 
substituting these values in given P.D.E we get 
initial value qd=- 7 = ■ 

4i 




y = -2p 2 +2q 2 
dt 

= pF p +qFq 


~~F X - P F Z = -p 


= -F y -qF q =-q 


Here — = 2— =>x = 2p+c, => x ~,2(p^v) . \ \ \ 

dt dt 1 v 7 \. 

— = -2— =>y = -2q + c 2 => y = -2q + V2 v \ ' 

* dt 2 


Also K\- : 

^(p + ?) = -(p + 9) => p + q=c 2 e-‘ \ 

Initial values: => 1+-2L = C3 => c 3 = l(V2-l) 

/> + 9=^(V2-l)e‘ ( ...(4) 

A.lso = -/>+? = -(/>-?) 

a/ 

-/ V V V ( r~ A 

=> P-q = c 4 e => ---— = c 4 ^ C4=_ 2^ 2+ ' 

P-9 = “(>^ +1 )«" , -( 5 ) 

Solving (4) & (5) we get 

2p = V < [V2-l-V2-l] => P = ~e~‘ 
and 2? = -V‘[V2-1+ V 2 + 1 ] => 

Substituting in equation (2) & (3) we get 

x = 2p + 2v = -ve~ l + 2v = v^2-e -7 j ...(6) 

y = -2# + y[2v = -^2ve~ t +^2v = V2v^l-e _1 j ...(7) 


■-K'fi-') 
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Using these values we have 


Solving (6) & (7) for v & e 1 we get 


Using these value in (8) solution becomes 

s 

/ rr \ 2 / _ rr \2 


it satisfies the given P.D.E. 

Example: Given P.D.E. F(x,y,z,p) = p 2 + q 2 -4z = 0 
Initial surface z = x 2 + 1 


Solution: Initial values are t n = 0. 


Characteristics equation are: 


9999183434 & 9899161734,8588844789 
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Linear Partial Differential Equations with Constant 
Coefficients and Equations with variable Coefficients 


LINEAR PARTIAL DIFFERENTIAL EQUATIONS 
WITH CONSTANT COEFFICIENTS AND »f 
EQUATIONS WITH VARIABLE COEFFICIENTS 

4.1 Homogeneous and Non-Homogeneous Linear Equations With A.\' 

Constant Coefficients 

'Y 3 £-TiSVjr. ‘ ‘ ' 

\ . > • - • ’• 1 V >' ; , • -• 

A partial differential equation in which the dependent Variable and its 
derivates appear only in the first degree and are not multipliedtogether, their 
coefficients being constants or functions of xand v^-is'-toown as a linear, 
partial differential equation. The general form of such an equation is 


. d n z d n z A d n z 

A)-t"A—;—I— i-Aj— 

dx n dx n ~ x Qy dy n 




a«-i 7 s/i-l 

dxT 1 V' 1 "V- 1 J 

+...+ M 0 ^ + M l ^ + N 0 z = f( x ,^^S ...( 1 ) 

dx dy) 

Where the coefficients A^,A b ....A n ,BQ,B b . ,B n _ b M 0 ,M\ N 0 are 

constants or functions of xand y. UAQ,A b ...A n ,B 0 ,B h . ,B n _ b M 0 ,M { 

and N 0 are all constants, then (1) is called a linear partial differential 
equation with constant coefficients. 


"V- 1 


- v 


dd V- 

For convenience — and — will be denoted by D( or D r ) and D'(ovD.,\ 

dx dy { ’ $01 

respectively. Then (1) can be rewritten as Mpff 


y-rjA;. 


[(AqD" +A l D n ~ l D’+ ... +A n D’ n ) [B 0 D nA + B l D n ~ 2 D '+... + 5„_ 1 D , ”“ 1 ) 


+...+ (M 0 D + Af l D') + N 0 ]z = f(x,y), 
Or briefly, F(D,D')z = f{x,y) 


...(3) §1 


When all the derivatives appearing in (1) are of the same order, then the ; ' 

resulting equation is called a linear homogeneous partial differential equation 

with constant coefficients and it is then of the ; - 


Form (AA>" + A x D nA D' + . . +A n D' n )z= f{x,y) 


...(4) 


On the other hand, when all the derivatives in (1) are not of the same order, A; 

then it is called a non-homogeneous linear partial differential equation with ■ 

constant coefficients. 
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Coefficients and Equations with variable Coefficients 


4.2 Working Rule For Finding Complementary Function of Linear 
Homogeneous Partial Differential Eq uatlon With Constant 
Coefficients 


Step l:Put the given equation in standard form 

(. A 0 D n +A l D n - l D'+... + A n D' n) jz = f(x,y ) . -0) 

Step 2: Replacing D by m and D by 1 in the coefficient of z ,we obtain 
auxiliary equation ( A.E .) for (l) as 

A 0 m n +A l m"~ 1 + . + A n =0(2) 

Step 3: Solve (2) for m, following cases will arise. \ 

Case (i): Let m = (distinct roots) V \ ■ ..... 

\. 'v \ Vv x.. > 

Then, C.F. = <t>i(v + m,x) + <|>2 (y + m 2 x) + ...- + < t’«C v + / ”« v )> v - 
Where are arbitrary functions \ . W 

Case (ii): Let D = m (repeated n times). Then corresponding to these roots 

C. F. then =<J)| (y + mx) + x<j> 2 (y + mx ) +.. -\ : + x (/ +mx ) 

Case (iii): Corresponding to a non-repeated fact° r D on L.H.S. of (2) the 
part of C.F. is taken as <j>(y). \ 

Case (iv): Corresponding to a repeated factor D on L.H.S. of (2) the part 
of C.F .is taken as s 'v—„ 

<l>i(y) +42 (t) +^ 2 <I ) 3(t) +. 

Case (v): Corresponding to a non-repeated factor D on L.H.S. of (2), the 
part of C. F. is taken as <t>(x). 

Case (vi): Corresponding to a repeated factor D on L.H.S. of (2) the part 
of C.F. is taken as 

(j>i (x) + y<t >2 (x) + y^<|)3(*) +. ... + y m 

Alternative working Rule For Finding C. F. 

Let the given differential equation be F[D,D)z - f{x,y). Factorize 
F(D,D r ) into linear factors of the form (bD-aD'). Then we use the 
following results: 

(i) Corresponding to each non-repeated factor (bD - aD ), the part of C.F. 
is taken as §(by + ax). 

(ii) Corresponding to a repeated factor (bD - a D) , the part of C.F. is 

taken as {py + ax ) + x<j) 2 (by + ax) + I 








x 2 §T > (by + ax) + ....+x n> { § m (by + ax). 


: 5 *- 


A ^3 \ u s / - rm 

(iii) Corresponding to a non-repeated factor D, P art of C.F. is taken as '..y\h tl’.Ci'iCi 
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(iv) Corresponding to a repeated factor D m , the part of C.F. is taken as 

<I>1 O')+42 (>')+* 2 <h (v)+• -x mA bm (>•) 

(v) Corresponding to a non-repeated factor D', part of C.F. is taken as 

<K4 

(vi) Corresponding to a repeated factor D ' m , the part of C.F is taken as 

<t>i (*)+yh (*)+ y\ (*)+. +y m ~\ (4 

Example: Solve 2r + 5s + 2t = 0 

Solution: Now, r= ~ = D 2 z,s = 4iL =£>D'z and t = ^-4r = D' 2 z . 

I dr J [dxdy) -dx 2 

Hence the given x 

. , ’’snO 

Equation can be re-written as 1 2D 2 +5DD'+2D ,2 \z±0'~-^ \ ' 

Its auxiliary equation is 2m + 5m + 2 = 0 or (2m± 1 2) = 0. 

V\ V s * 

So m — —1 12,-2 and hence the general soluti6n.qf'(li,is 

\ ' s \Jn' ~ 

z = §\(y-xl2) + <}> 2 (>■ — 2x ),(jjj and <t» 2 beingarbitraty.functions . ...(2) 

. X-— 

Let § x (y-xl 2) = <^ {(l / 2)(2_y - *)] = y j (X - (2)becomes 

z = \j/ x (2y - x) + (j) 2 (y - 2*), \|/j and <j> 2 Being arbitrary functions. 


(~A \ ,4 V~ 

Example: Solve — \ - —\ =0 \ \,x 
IcbrJ Idy I v. 


t 4 I ftv 4 


Solution: Re-writing the given equation is (D 4 -D A )z= 0 ...(1) 

Its auxiliary equation is m 4 -1 = 0 or [m 2 - lj^m 2 + lj = 0 

So that m = 1,-1,/,-/ where i = \f-l 
Hence the general solution of (1) is 

*■■=■ 4>1 {y - x) + <t »2 {y + x) + «t >3 (y + ix) + <t> 4 (y-ix), 

Where <t»i,<t> 2 ,<t» 3 and <J > 4 are arbitrary functions . 

Method of Finding Particular Integrating (P.I.) of a Linear 
Homogeneous Partial Differential Equation 

F(D,D')z = f(x,y). 

Method-I. When f(x,y ) is of the form <t>(a*+6y). 

If F{D,D') be homogeneous function of D and D' of degree n, the 

1 . tnl / i \ ; 


f(d,di 


tsy 1 ’ (ax+by) 


F(a,b ) 


<|>(ar + hy). 


Provided F(a,b)* 0, b^ittg the nth derivative of <|) w.r.t .ax + by as a 
whole. 
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Linear Partial Differential Equations with Constant 
Coefficients and Equations with variable Coefficients 


Exceptional case: When F(a,b) = 0, then the above theorem does not hold 
good. In such a case the new method is based on the following theorem. Note 
that F(a, b) = 0 if and only if (bD - aD') is a factor F(D,D'). 


I tv ,-.i 


-<b (ax + by) = -<J>(ax + by). 

(bD-aDf ^ b n n\ 

Working Rule For Finding Particular Integral Where 

f( x >y) =<)»(<« +b y) 

Formula (i): When F(a,b)* 0 and F(D,D') is a homogeneous function of 
degree n, then we have 

P. I. = ---<b(ax + by) =—-— [f....f f(v)dv dv. dv ..where v = ax + by 

f(d,d') k ’ F(a,b) }i J . \ 

note that R.H.S. contains a multiple integral of «th order.' ' ./' N; - 


-<b( ax + by) = -<|> (ax 

(bD-aDf 7 b n n\ V \ ' 


\ \ \ 

Example: Solve {d 2 + 3 DD ’+2 D' 2> jz = x + y \ \\ 

v \ \ 

Solution: The auxiliary equation of the giveftequation is 


m 2 + 3m + 2 = 0 giving m = -1,-2. 


\ 


C. F. = <)>! [y - x) + <j> 2 (y - 2x) ,^ ^'bein^arbitrary functions. 
1 

Now, P.I. = — -+ 

D 2 + 3DD'+2D' 2 v 


l 2 +3.1.1 + 2.1 2 


|| v dv dv, Where v = . 


- —| — Jv = (l/6)x^v 3 /6j = (l/36)(x + y) 3 . 

6 l 2 J 

Or z = § l (y-x) + $ 2 (y- 2x) + (l/36)(x + y) 3 . 

1/2 

Example:r + s-2t = (2x + y) 

Solution: The given equation can be re-written as 

^4 + --2^4 = (2x + y) 1/2 or (D 2 +DD'-2D ,2 )z = (2x + y) m 
dx dxdy dy 

Its auxiliary equation is m 2 + m-2 = 0 so that m = 1,-2 

/. C.F. = 4> t (y + x) + <j) 2 (y - 2x) ,<J>j, <j> 2 being arbitrary functions . 


Now PJ. = —---(2x + y ) 1 

D 2 + DD'-2D' 2 v 

= ———- , ff v n dv dv , where v = 2x+ 

2 2 +2.1-21 2 JJ 
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Linear Partial Differential Equations with Constant 
Coefficients and Equations with variable Coefficients 


=—f-v 3/2 Jv= —v s/2 = — (2 x + >0 5/2 
4 J 3 . 4 3 5 15 

Hence the required general solution is y = C.F. + P.7. 
or z = <t>! (y + x) + $2 ( y - 2x) +(1 / 15)(2x-+ yf 2 . 

Example : (D 2 -3DD' + 2D' 2 )z = e 2x ~ y +e x+y + cos (x + 2y) 

Solution: Here auxiliary equation is m 2 - 3m + 2 = 0 so that m = 1,2. 

C. F. = 4>i (>’ + x) -r 4> 2 (y + 2 x),4>i,<}) 2 being arbitrary functions. ...(1) 
Now, P. I. corresponds to e 2x ~ y \ " x \ 


D 2 -3DD ,l +2D' 2 
v=2x-y 


o lx ~y 


2 2 -3.2(-l) + 2(-l) •; 


— \e v dv = (\!\2)e v = {\l\2)e 2x+2y 


Again, P.I. corresponding to e x+y 


D 2 -3DD'+2D' 2 


e x+y = 


D-D'lD-iD' >| 


where 


—y-j-JeVvi, where v = x usiri^ formula (i) 


D-D' (D-D') 1 


? x +y - -r a x+ y 


Finally, P.I. corresponding to cos(x + 2y) 

= —z ---r-cos(x + 2y)= --—-—ffcosvt/v dv , where 

D 2 -3DD'+2D' 2 v ’ l 2 - 3.1.2+ 2.2 2 JJ 

v = x + 2y ...(4) 

=-fsinv^v = -- cos v = ——cos (x + 2y). 

3 J 3 3 V ’ 

From (1), (2), (3) and (4) the required solution is z = C.F. + P.I. 

or z = 4*!(y + x) + § 2 {y + 2x) + -^e 2x ~ y -xe x+y --^cos(x + 2y). 

Method II. When /(x,y) is of The Form x m y n or A Rational Integral 
Algebraic Function of x and y. 

1 D' 

If n<m, should be expanded in powers of — whereas if 

m < n,—— - should be expanded in powers of —. 

f(D,D) 1 1 D- 
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Linear Partial Differential Equations with Constant 
Coefficients and Equations with variable Coefficients 




Example: Solve (D 2 - a 2 D 2 )z = x or ^4 _ a 2 —4 = 

dx 2 8y 2 

Solution: Here auxiliary equation is so that m= a,-a 

:. C.F. = <}> 1 (y + ax) + <j> 2 (y - ax) , 4>!, 4*2 being arbitrary functions ...(1) 



Now, PJ. = 


D 2 -a 2 D' 2 D 2 [\-a 2 (D a /D 2 )] D 


1 f, ^ D' 2 

x=—r 1 -a —- 
D 2 { D 2 


1 a 2 D' 2 1 x 3 

— 1 + =— + ... x = —x = — 

D 2 [ D 2 D 2 6 


Alternatively, we can compute P.I. as follows: 
1 1 


P.I. =- 


-a 2 D' 2 -a 2 D' 2 1 ~(D 2 /a 2 D' 2 ) x v a 2 ©ll[ a 2 D' 


'.w x' x 




a a U' j a u a\ A \.\ 

\ 

Hence the required general solution is z = ClFx-t A P. 

z = <t>i {y + a*) + <i >2 (y - a*) + * 3 / 6. usiri^'O) ah^2.) • 

or z = <j>[ (y + ax) + <t> 2 (y - ax) - (xy 2, j 7|-2a^J;^sing (1) and (3) 

A General Method of Finding the Particular Integral of Linear 
Homogeneous Equation with Constant Coefficients 

PL ~ (D-m l D , )(D-m 2 D , )...(D-m H D , )^ X,y ^ 

We shall use one of the following formulas: 

Formula I:--—- f(x, y) = \f{x,c- mx)dx where c-y + mx ...(2) 

D - mD' 

Formula II: — - f(x,y ) = f f(x,c + mx)dx, where c = y-mx ...(3) 

D+mD \ 

Hence in order to evaluate P.I.(\), we apply (2) or (3) depending on the 
factor D-mD 1 and D + mD'. Note that result (3) can be obtained from (2) ’ 
by replacing mby-m. 

Example: Solve — + — = sin x. 

dx dy 

Solution: Re-writing, the given equation is (D + D')z - sin x ...(1) 

C.F. = <|>(y-x), Where <|> is an arbitrary function 


and P.I. -= sin x = f sin xdx = -cosx 

D + D 1 J 

Hence the required solution is z = C.F.+P. I = <j)(_y - x) -cosx 
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Linear Partial Differential Equations with Constant 
Coefficients and Equations with variable Coefficients 


Coefficients: 

A linear partial differential equation which is not homogeneous is called a 
non-homogeneous linear equation. 

Consider a linear equation with constant coefficients 

F(D,D')z = f(x,y) ...(2) 

If (1) is a homogeneous equation, then F(D, D') can always be resolved into 
linear factors. On the other hand if (1) is not homogeneous then F(D,D') 
cannot always be resolved into linear factors. 

Hence, we classify non-homogeneous linear equations into the following 
types: * \ \ 



. \ 

(i). When F(D,D') can be resolved into factors each of which is first degree 
' D'. For example consider x^lhe'-. equation, 


in 


(D 1 - D' 2 + 2D +1 )z = x 2 y 3 i.e. (D + D' + 1)0<D' + l)z.= x 2 y 


2 3 


(ii).When F(D,D') cannot be res61yd3“into. line'af'fqrnK,(a£) + bD' + c) For 
example, consider equation 

(2D 4 -3D 2 D' + D' 2 )z = 0 i.e. (2D 2 - D')0 2 -D)z = 0. 


Working rule for finding C.F. of non-homogeftepds equations which are 
reducible into linear partial differential equatiQn,,with constant coefficients, 
namely, F(D, D')z = 0. 


v 


Method I When F(D,D') can be expressed fts product of factors linear in 
D and D'. Then we apply the followin^qjes: 

(i) Corresponding to each non-repeated factor (bD - aD' - c), the part of 
C.F. is taken as e( CJc/ *)<j>(6y + ax), if b *0. 


(ii) Corresponding to a repeated factor (bD-aD'-c ) r , the part of C.F. 1 
is taken as (assuming that 6 * 0 


e (cr/i)[<^ ^f}y+Qx(+x§ 2 (by+ax)+ x 2 (j> 3 (6y+ax) +. +x r l $ r (by + ax)]. 


r-h 


(iii) As a particular case of (i) with c = 0, corresponding to each non- 
repeated factor (bD- aD'), the part of C.F. is taken as §(by + ax). 

(iv) As a particular case of (ii) with c = 0, corresponding to each repeated 

factor (bD-aD') r , the part of C.F. is taken as 
(pi (by + ax) + x <(>2 (by + ax) + x 2 <t> 3 (by + ax) + . + jc r—1 <p r (by + ax). 

(v) As a particular case of (iii) with a = 0,6 = 1, corresponding to each 
non-repeated factor D, the part of C.F. is taken as <{>(y). 

(vi) As a particular case of (iv) with a = 0,6 = 1, corresponding to each 
repeated factor D r , the part of C.F. is 


<|>i(y) + x<|> 2 (y) + x 2 <|>3(y) 


+_+ X 


r -1 


Mt)- 


(vii) Corresponding to each non-repeated linear factor (bD-aD'-c), the 
part of C.F. is taken as e~( c y ,a )§(by + ox), if a 0 


■ . ■ 


H§ ■; ■ 




■ 

■ 


■ 
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Linear Partial Differential Equations with Constant 
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• (viii) Corresponding to each non-repeated factor («£)' + c ) 5 the part of c.F. 
is e^) /a i(ar), when a* 0. 


(ix) Corresponding to repeated factor (aD 1 + c) r the part of C.F is 

e‘(°') /a [(j) 1 (ax) + y$ 2 (ax) + y 2 <h (ax) +.+ y r ~\|>, (ax)] 

(x) As a particular case of (viii), with c = 0, a = \, corresponding to 
repeated factor D', the part of C.F.=(j), (x). 

(xi) An a particular case of (ix), with c = 0, a = l, corresponding to 

repeated factor D ' r , the part of C.F. is taken as 

(*) ■+ y$i W+ y\ W ■+.+/~‘i (*) • 

Method II: When F(D,D') Cannot Be Resolved Into 
Z>and D' 

N. ' x 

In such case we apply trial solution as explained the following-example. 
Example: Solve (2D 4 - 3D 2 D' + D' 2 )z^V. N 
Solution: The given equation can be re-written as \ \ 

\ v_- 

(2D 2 - D')(D 2 - D')z = 0 \ V \ S -O) 

Consider (2 D 2 -D')z= 0 \ S’' ...(2) 

Let z = Ae hx+ky be a trial solution of (•])'.'then'we*have 

\ 

D 2 z = Ah 2 e hx * iy and D'z = Ake 1 ** 19 . Puttingthese values in (2), we get 
A(2h 2 -k)e hx+ky =0 so that lh 2 -k= 0 or k= 2 h 2 . 

Hence the most general solution of (2) is z = ^ Ae la+2l,2y .. .(3) 

Next, consider (D 2 -D')z = 0 ...(4) 

Let z = A'e hx+ky be a trial solution of (1). Then we have 

D 2 z = A'h' 2 e h x+k y and D'z = A'k’e h x+k y . Putting these values in (4), we 
get 


A'(h' 2 -k')e h ' x+k ' y = 0 So that h u -k' = 0 or k' = h' 


S h * 

A'e 

From (3) and (5), the most general solution of (1) is 


h'x+h'^ y 




t'e h x+h y ,A,A’,h,h',k,k 'being arbitrary constants. 


Methods of Finding P.I. of Non-Hontogeneous Equations with Constant 
Coefficients. 

F(D,D')z = f(x,y) 

CaseI: When f (x,y) = e ax+by and F(a,b)ytO. 
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m,m F(a,b) 

Thus in this case we replace D by a and D' by b 
Case II: When /(x,y) = sin(ax + £iy)or cos (ax + b 


Then, P. I. 


or cosax + 


F(D,D') 


Which is evaluated by putting D 
provided the denominator is non-zero. 

Case III: When f[x, y ) = x m y n 


Then P. I. =- x m y n = [F(D,D')]~ l x m y n , \ 

F(D,D') 1 V ■ y ’ 

Case IV: When /(x,y) = Ve ax+by , when V is a function of x'and y 
The P. I. =--- Ve myhy = L ^ 

Example : Solve (£> z -D a +D-D')z = e M \ 'Njx" 

\ ’ \ 

Solution: The given equation can be re-writteqas 

[(£> - D')(D + D') + (D- D')]z = e 2 *^ or (B~ D)(B + D' +1 )z = e lx+3y 

_ \ 'Ss''*— 

C.F. =(j) 1 (y + x) + e ^^2 (y ~x),<j > lH <|>2 being arbitrary functions. 


F(D,D') 


,2x+3y \ 


and P. I. =-_ 

(D-D%D + D'+l) (X>4X2 + 3 + l) 

Hence the required general solution is x = CF. + P. I. 


Solution: P. I. corresponding to e 


P. I. Corresponding to 5 cos 


cos(x + 2 y) = 5 


cos(x + 2y) 


cos(x + 2y) 


cos(x + 2y) = -5 


(Z) 1 — l)cos(x + 2y) 


Required P.I 


2sin(x + 2 y) - cos(x + 2y) 
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Partial Differential Equations with Variable Coefficients 
Reducible to Equations w ith Constant Coefficients 

We propose to discuss the method of solving the so-called Euler-Cauchy type 
partial differential equations of the form 


having variable coefficients in particular form (namely the term 


is multiplied by y 


multiplied by x‘ 


multiplied by x r y"~ r and so on. If D 


written as 


X V ' 

such - v -^equalipns are: 

y 2 D' 2 +xD- yb '•*= log\\ 

Method of Reducing Euler-Cauchy Typ'e.Equatibn to A Linear Partial 
Differential Equation with Constant Coefficients' 

Consider Euler-Cauchy type equation \ N 

(a 0 x n D n +a,j t"-' yD"-'D' + a 2 x"- 2 y 2 bVb'\\. + a n y n D' n +..)z = f(x,y) 


Examples 
(x‘D‘ - y 


Where D 


Define two new variables u 
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and so on. Similarly, we have 


Linear Partial Differential Equations with Constant 
Coefficients and Equations with variable Coefficients 


fSIP 


yD' = D' x , y z D a = D{(D(- 1) y 3 D' 3 - A'(A'-1)(A'-2 
and so on. Also, we have xyDD' - D X D[. 


and x m y n D m D' n = 0,(0,-l)...(0, -m + l)0,'(o,'-ll.Jo,' -n + l). ...(10) 


Using the substitutions (2) and results (4), (5), (6), (7), (8), (9) and (10), the 
given equation (1) reduces to an equation having constant coefficients and 
now it can easily be solved by the methods already discussed for 
homogeneous and non-homogeneous linear equations with constant 
coefficients. Finally, with help of (2), the solution is obtained in terms of old 
variables * and y. \ N. \ 

4.4.2. Working Rule For Solving Euler-Cauchy Type .partial -Differential 
Equation ■„ 

(a 0 x n D n + a x x nA yD nA D'+ a 2 x n ~ 2 y 2 D n ~ 2 D ' 2 +...+ a„>-"D' n +v: .,)z^/(x,y) (1) 

\ \ \ \ % 

Step-1. Introduce two new variables u and v \ ' 

\ 

x-e u and y-e\ i.e., m = log x and v = log y \ \ ...(2) 


Step-2. We have O = — and 

dx 


D’ = —.''Also.let D X =—,D X = — 
. dydu dv 


Step-3. Use the following results in (1) \ 

. 

xD = D l ,yD' = D{ ,x 2 D 2 =D X (D X -%y 2 D' 2 = D{ (D x '-l) 1 ^ 

x 3 0 3 = A(A -1)(A -2),y 3 O l3 =D{ (D 1 , -i)(A’-2)andsoonj 

Step-4. Using (2) and (3) in (1) , we obtain the following linear partial 
differential (homogeneous or non-homogeneous) 

Iu n" . u . u n'-2n'2 . . I. n'» . )_ rA\ 


b 0 D," + b x D nA D[ + b 2 D?- 2 D[ 2 +... + b n D[ n +... z = g(u, v) 


Step-5. If (4) is homogeneous linear partial differential equation, then it is 
solved with help of previous method. Again, if (4) is non-homogeneous 
linear partial differential equation, then it is solved with help of previous 
methods. 

Step-6. Using u = logx and v = logy in the solution obtained in step 5, we 
finally obtain the required solution in terms of the original variables x and y 

Example: Solve x r - 3 xys + 2y t + px + 2qy = x + 2y 
Solution: The given equation can be re-written as 

2 fd 2 z\ f d 2 z ^ o 2 T fdz^ fdz^ 

[dx ) {dxdy) {dy ) {dx) [dy) 

Or (x 2 Z) 2 - 3xyOO' + 2y 2 D' 2 +xD + 2yO')z = x + 2y. ...(1) 

Let x = e", y = e\ so that u = log x, v = logy. .,.(2) 

r) r) r) 

Also, let Os—, D' = — , ft s— and n' = _. 

dx dy 1 du 1 dv 
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(1)Becomes [B,(B, -1)-3B,B,'+ 2D[{D[ -l) + D,+ 2D ,']z = e u +2e v . 

Or [pi -2D X D[ + 2D [ 2 )z = e“ + 2e v or (D x -D[)(D l -2D[)z = e u + 2e v . 

C.F. = ^(v + w) + ^ 2 (v + 2w) = $(logy + logx) + ^ 2 (logy+ 21ogx) 

Or C.F. 4> } (logxy) + <t> 2 (logx 2 y) = f\(xy) + f 2 (x 2 y) where f x and / 2 are 
arbitrary functions. 

Also, P.1, = 

_1_ ( e u +2e')= -—_-_ e lJ ' +0v +2—---e a " +lv 

(q-q)(q-2q) v ’ (q-q)(q-2q) (q -q)(q -n) 

1 1 • '"'-X \ 

= --7-7 e u +2 7 - 77 - 7 e' = X + y, using (2) Vx.' 

(l-0)(l-0) (0 —1)(0 — 2) 7 6W 

'w, V" 

Hence the required general solution is z=C.F>+P.I. or. 

z = f l {xy) + f 2 (x 2 y) + x + y. 

Example: Solve x 2 r - y 2 t + px-qy- logx. \ \ s 

or (x 2 D 2 -y 2 D' 2 + xB-yB')z = logx. \ 

, '"*x v v—' 

Solution: Let x = e“, y = e v so that u = logx. v^lo^y. ...(1) 

Also, let D =—, D' = —, B, = — and B,' = — 
dx dy 1 du ■ 

\ \ 

Then the given equation ( x 2 D 2 -y 2 D' 2 + xD ^yB'jz = logx becomes 

[b 1 (b 1 -i)-b;(b;-i)+b,-b;]z=m or (b 2 -b; 2 )z=u 

or C.F. = <|) 1 (v-M) + (t) 2 (v + «) = (|) 1 (logy-logx) + <t) 2 (logy + logx) 


or C.F. = <j)| (logxy) + <j> 2 1 log— =/; (xy) +/ 2 — where f x and f 2 are 


arbitrary 


P.I. = - -ju = —7-- 

D 2 -D 2 ( n' 2 

1 1 D 2 1 — —U- 

U ' 1 D 2 


functions. an 

i (. b; 2 Y‘ i ( D' 2 ) 

U =—r 1— r U =—r- 1 H-7- + .... \U 


D\ \ B, 


D l A D 2 







= ~T m= 7 m3 =~(l°gx) 3 ,by (1). 
JJ, DO 


/. Required solution is z = /j(xp) + / 2 — + - x(logx) 3 ,/j,/ 2 being 

UJ v6; 

arbitrary functions. 

Example: Solve (x 2 D 2 -Ay 2 D' 2 -AyD' -l)z = x 2 y 2 log y. 

Solution: Let x = e“, y = e* so that u = logx, v = logy. ...(1) 

Also the given equation reduces to 

[B, (B, -1) - 4B; (b; -1) - 4B,' -1] z = e V v v 

or (B 2 -B,-4B; 2 -l)z = e 2 “ +2v v. ...(2) 
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REGENERATING LOGICS fll. 


Linear Partial Differential Equations with Constant 
Coefficients and Equations with variable Coefficients 


Here (Z), 2 - Z), - 4Z),' 2 -1) cannot be resolved into linear factors in D, and D[ 
. To find C.F. corresponding to it, we consider the equation: 
(Z) 2 -Z) 1 -4Z); 2 -l)z = 0. ...(3) 

Let a trial solution of (3) be z = Ae hu+kv and 

D' 2 z = Ak 2 e hu+kv , D l z = Ahe hu+kv - ,Dyz = Ah Z e hu+kv ,..(4) 

Then (3) => A(h 2 -h~4k 2 -l)e hu+kv = 0 => h 2 -h-4k 2 -1-0. ...(5) 

C.F. of (2) = £ Ae hu+hv = ^(e")* =£>V 



D 2 x -D x -4D' 2 -\ 


g 2 “ +2v v = e 2 " +2v 


(D, + 2) , -(D,*2A(^2) ! -l 


_ ^ 2 u+ 2 v ____ 

D 2 + 3Z), - 4D,' 2 -16 D[ -15' 


\\ N 


(- l5 ) 4 1 4!!l * °.T D ;> iiL n, vf ilxfl, - f JrL A* 


=c 2u+2v —i+j —d[+—d ; 2 --D r - -Lz), 2 
(—15) L 115 15 5\’Sn 


<?*'(, 16 n , ) e 

(- 15)1 15 ‘ J (- 


(- 15)1 15 


(e“) 2 x (e v ) 2 (l 6 — 15v) 


:x 2 y 2 (16-151ogy),using (1). 


1 \ 

The required general solutionis z = ^Ax h y k + —— xx 2 y 2 (l6-151ogy) 

\225J 

where h 2 -h- 4k 2 -1 = 0, and A , h and k are arbitrary constants. 
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Second Order Partial Differential Equation 


CHAPTER-5 

CLASSIFICATION OF SECOND ORDER PARTIAL 
DIFFERENTIAL EQUATIONS, CHARACTERISTIC 
CURVES AND REDUCTION TO CANONICAL 
FORMS 

Consider equation of the type Rr + Ss + Tt+f(x,y,z,p,q) = 0 ■ ...(1) 

\ '''X „\ 

Where R,S,T are continuous functions of x and y possessingxCpntinuous partial 

• N >N \ 

derivatives of as high an order as necessary. Laplace transfonnatioh pn*(l) consists 

of changing the independent variables x,y to new set of independent variables u,v 

v • 

where , nn* 

S v -^ 

u = u[x,y) and v = v(x,y) ^ \ •••(2) 

\ N\ 

are to be chosen so that the resulting equation ip independent variables u,v is 
transformed into one of three canonical forms which arepasilylntegrable. 

\ ^ 

From (2), we have 

\ \ 

_dz _ dz du + dz dv _ dz _ dz du + dz 


dx du dx dvdx 


dy du dy„ dv*dy \ 


d__dud_ dvd_ d__du_d_ dv_d_ 
dx dx du dx dv ’ dy dy du dy dv 


d l z d r dz\_(du d + dv d^fdudz + dvdz 
dx 2 dx\dx) V& du dx dvJvdx: dxdv 


du d fdu dz dv dzY dv d f du dz dv dz 
dx du l dx du dx dv J- dx dv t dx du dx dv 


_ d 2 zf duY ^ d 2 z du dv + d 2 zf dv^j 2 dz d 2 u dz d 2 v 

du 2 \dx) dudv dx dx dv 2 \dx) du dx 2 dvdx 2 

d 2 z d f dz(du d dv 5 V dz du dz dv^l , 

dx dy dx\dy) \dx du dxdvj^dudy dv dy J 

_d 2 z du du d 2 z r du dv + du dv > d 2 z dv dv dz d 2 u dz d 2 v 

du 2 dx dy dudvydx dy dy dx) dv 2 dx dy du dydx dv dydx 


...( 4 ) 


and t = 


d l z d f dz 


dy 1 dyydy) \dy du dydvjydydu dy dv 


du d + dv d du dz ^dv dz 


, by (3) and (4) 


_dud(dudz dv dz j dv d ( du dz dv dz v 

dy du ^dy du dy dv > dy dvydy du dydvj 

_d 2 z'du^ +2 d 2 z du dv + d 2 z r dvV dz d 2 u dz d 2 v 

du 2 {dy) dudv dy dy dv 2 [dy J dudy 2+ dvdy 2 
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' Putting these values of p,q,r,s,t in (1) simplifying, we get 

, d 2 z d 2 z d 2 z f dz Sz'l 

a 'Z2 +2B ^ +c T 2 +f \ =0, 

du dudv dv V du dv) 

f du^ du du 

Where ^ = — +S - +T — , 


dx) dxdy l dy 


n n dudv 1 J du dv dudv) „dudv 

B = R - +-S -+- +T -, 

dx dx 2 dxdy dydx) dy dy 


c=«f—f+s——+rf—T 


dx) dxdy l<9>> 




...( 6 ) 




...( 8 ) : 



* o!z dz ^ v \. 

and F u,v,z ,—,— is the transformed form of f(x,y,z,p;q-\-. ' \ 

du dv) v\ \ v - '"v ff# 

Now we shall find out u and v so that (5) reduces to. simplest possible 

form. The method of evaluation of desired values of i} and ^ becomes easv 

2 --\ '\ \\ 3 '&:%?; . 

when the discriminant S -ART of the quadratic equation 


rx 1 +sx+t=o 


\ w 


is everywhere either positive, negative orzero,und now we shall present '' 

these three cases separately. < — tiii- i : ' 

> — - \ •... : : 

2 \ ' ••••' *' - 

Case I: Let S - ART> 0. When this condjti.onjs satisfied, then the roots 

X x ,X 2 of the equation (9) are real and distinct. The coefficients of 

\ S'*. *5' ~ 

d 2 z d 2 z V. ■• r • 

—-and—- in the equation (5) will vanish iPwe choose «and v such that &>••':> 

du 1 dv 2 x .’ 


du „ [ du 


,‘r Mdydd 


dx l dy 


dv (dv 

dx dy 


...( 10 ) 




Since q is a root of (9), we have RX 2 + SX x +T = 0 


Using (10), (6) gives A = (rX 2 + SX x +t) — = 0, by (12) 

\ ) rh; 


Again, since X 2 is a root of (9), we have RX 2 + SX 2 + T = 0 


Using (11), (8) gives C=\RX 2 2 + 


SX 2 +t) — =0 by(14) 
\dy) 


Re-writing (10), we have 


du) „ ( du 


dx) {dy 




...(ii) 

...( 12 ) | 

®tS -: 

•••(13) : 

•••( 14 ) ' V 

•••(is) 

-n ii; 


f=-r = j ...(17) ay awatswi 

Taking third fraction of (17), du = 0 so that u=c x , q being an arbitrary 

constant -(■*) ii§lt88 


constant 
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Taking first and second fractions of (17), 


+Xt — 0 


Let the solution of (19) be f\(x,y) = c 2 , c 2 being an arbitrary constant...(20) 
From (18) and (20), the general solution of (16) [i.e. (10)] is 
u = A(x,y) ...(21) 

Similarly, the general solution of (11) can be taken as 

v = f 2 (x,y) ...(22) 

We can easily verify that 


AC-B 2 =—I 


( 4 rt-S 2 )^-^ 

\ \dxdy dydx 


B 2 =-(s 2 -4RT 


\( dudv du dv 
\dxdy dy dx 


as A = C — 0 


...(23) 


Since S 2 -RT> 0, (23) shows that if > 0. Hence we may divide both sides 
of (5) by B . Then nothing that A = C- 0, (5) transforms of the form 

d 2 z f dz dz) \~... 

-< 24 > 

which is the canonical form of (1) in this base. V ‘ 

\ ) 

Case II :Let S -4RT=0. When thisxjcofrditic>n is satisfied, the roots ,/. 2 
of (9) are real and equal. We now take « exactly as in case I and take v to 
be an y function of x,y which is mdeperident of u . We have, as in case I, 

A=0. Also, since S 2 -4RT = 0 (23) shows that B 2 = 0 so that B= 0. 

Moreover in this case C* 0, otherwise v would be a function of u and 
consequently v would not be independent of u as already assumed. 

Putting A = 0, B = 0 and dividing by C, (5) transforms to the form 


= <t> K,V,Z,- 


...(25) 


Which is the canonical form of (1) in this case. 


Case III:Let S -4 RT <0. when this condition is satisfied, the roots X x ,X 2 
of (9) are complex. Hence is case III is formally the same as case I. 
Therefore, proceeding as in case I, we find that (1) reduces to (24) but that 
variables u,v instead of being real are now complex conjugates. To obtain a 
real canonical form we make further transformation u = a+i'P and v = a-i'P 
so that 

f lV . „ i(v-u) 

a =k ( m+v )>P=^T~ -(26) 


XT dz dz da dz 3(3 If dz dz ) 

du dadu dpdu 2(,da dp J M ’ 

, dz dz da dz dp 1 f dz . dz ) 

dv dadv dp du 2^ da dp j y 


...(27) 


...(28) 
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Second Order Partial Differential Equation 


I 


dz d f 8z\ if d . 3^1 f dz . dz) u a 

—- = T- — r-i7r.:r+'^ , by (27) and (28) 


dudv du^dv) 2\da 3p/2(da dp 


1 \ d ( dz .dz) . d f dz .dz 




‘4[da[da + '3pJ 'dpUa^dp 


\[&z_ -dh__.fz_ .d^z 

4 ( da 2 dadp ‘ dpda + * gp 2 


z 1f d l z d 2 z 


dudv 4l v 5 a 2 + dp 2 


d 2 z _ d 2 z 
3a3p 3(3 3a 


...(29) 


Putting u = a+/p, v = a-/p and using (27), (28) and (29), (24) reduces to 
+ ~~ = y f a, p, z—, —\ which is the canonical form of (1). in this 


da 2 dp 2 


da dp 


Remark: Second order partial differential equations-of the type (1) are 
classified by their canonical forms. An equation oft he form (1) is said to be 

(i) Hyperbolic if S 2 -4RT>0 \ 

\ ' \ 

(ii) Parabolic if S 2 -4RT = 0 ^ ~ 

(iii) Elliptic if S 2 -4RT <0 . , \ x 


Example: Reduce ^-4=x 2 


to canonical form. 


Solution: Re-writing, the given equation becomes r -x 2 t = 0. 
Comparing (1) with Rr + Ss + Tt + f(x,y,z,p,q) = 0, we have 


R = 1, S = 0, T = - 


Now, the X -quadratic RX 2 +SX+T -0 gives 
X 2 -x 2 =0 so that X = ±x. 

.-. Here X, = x and X* = -x 


(Real and distinct roots) 


Hence ^+A,, = 0 and — + X~> =0 become 
dx dx 1 

—+x = 0 and — -x = 0 
dx dx 

Integrating these, y+(l/2)x 2 =c x and y-(l/2)x 2 =c 2 . 

Hence in order to reduce (i) to canonical form, we change x,y to u,v by 
taking. 


- ■” ‘ 

■ 


u = y + (l/2)x 2 and v = y-(l/2)x 2 . ...(2) 

dz dz du dz dv dz dz . ' 

dx du dx dv dx du dv 

dz dzdu dzdv dz dz . ... 

^ duty dvdy du dv m0-Wm§- 
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oz_d(dz}d\ (dz 5 T & ^ 

faydcj Sc{ \du dv )J 3rv& dv) [ft/ dv 

d fdz dz\du^d(dz dz') dv + dz dz 
du\du dvjdx dv\du dv) dx du dv 


d f dz dz 


(dz dz 
+ 1.- 


d 2 z ^ d 22 + d 2 z +dz dz 
du 2 dudv dv 2 j du dv 


d 2 z d dz ( d d V dz _ d 2 z v z + v z 
an / dy 2 fy\dy) dv){du dv) du 2 dudv dv 2 

Putting these values of r and t in (1), we get 


•; • 




d 2 z . d l z d l z ) dz dz 


X \^du 2 dudv dV J du dv 

d 2 z 1 f dz dz \ d 


rsl ^2 *0 

d Z \ 2 9 Z i — 2 

& 2 dudv 


f dz 

dz' 

&z 1 { 

r\-r — 

dz 

'•'dz ' 

v du 

dv ) 

U1 — 

dudy A(il-v){ 

■ 9 < 

nS 


Which is the required canonical form of the givencquafion. 


d 2 z ( d 2 z \' d 2 z ■ 

Example: Reduce the equation — r + 2 —+-^-0 to canonical form 

dx ydxtyj'-ty^ 

and hence solve it. ~ 'V.. 

Solution: Re-writing the given equation, we get r + 25 +1 = 0 ...(1) 

Comparing (1) with Rr +■ & + Tt + / (x,y,z,p,q) = Q, here 

R = 1, S - 2, T -1 so that S 2 - ART = 0, showing that (1) is parabolic. 

The X -quadrate equation reduces to l 2 + 21 +1 = 0 so that A = -1,-1 
(equal roots). 

rTLrt nrtwafiwnnrlmrr oliorn nti o an noti An ip f 1 — O ar //v At i — A 


The corresponding characteristic equation is 


-1 = 0 or dx-dy = 0 


Integrating, x - y = c, c being an arbitrary constant. 

Choose u-x-y and v = x + y, ...(2) 

Where we have chosen v = x + y in such a manner that u and v are 
independent functions as verified below. 

du du 

a, ^ = a,av_a<av =ll+1U2#0 

5(x,^) dv dv cbc dy dy dx 
dx dy 

dz dz du dz dv dz dz . ... 

Now,p = —= using (2) ...(3) 

dx du dx dv dx du dv 

dz dz du dz dv dz dz . ... ... 

q = ~dy~~du 'dy + d^'dy~~du + dv’ USmg ^ 

_ /,•, j ,a\ 9 d d d d d ... 

From (3) and (4) , — = — + — and — = — + — ...(5) 

dx du dv dy du dv 


. 

: - . 
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'dz d [ cz i d d\(dz 

r = — r = — -T- = —+—— + —, by (3) and (5) 
Sc 2 dx\dx) dv)\du dv) 

d (dz dz'l d (dz ^ Sz'l d 2 z <3 2 z + d 2 z 
du\du dvJ dvydu dv) du 2 dudv dv 2 

d 2 z d f dz) d dV dz fiz^ , , 

dy dy\dy) du dv)\ du dv) 

d_( dz dz'\ d( dz _ | dz]_ tfz 2 g 2 v ^ g 2 z 

du 1 du dv) 3vl du dv) du 2 dudv dv 2 


and j = — = — — H —+— + — by (4 

dxdy dxydy) ySw dv)\ du dv J 

_d( to dz} d ( dz dz} d 2 z d 2 z 

du V du dv) du dv) du 2 dv 2 x 

Using (6), (7) and (8) in (1), the required canonic^ldorh^is^ 


d 2 z d f dz^ 

—- = 0 or-=0 

dv 2 dvidvj 


To Find The Required Solution, integrating\^9)‘partiairy w. r. t. V, 
dz 

we get — = Mu), ^being an arbitrary fiinction.\ ...(1 

dv i 

Integrating (10)partially w. r.t. V, i^\<j>\u)dv + y/(u) = v<l)(u) + w(u) 


or z = (x + y)^(x-y) + y/(x-y), which is the desired solution, <f>, if/ being 
arbitrary functions. 

Example: Reduce the following partial differential equations to canonical 
forms: 

ps2 f p\2 ^ 

—j + x 2 — \ =0 or r + x 2 t = 0 

to 

Solution: Re-writing the given equations, we get r + x 2 t = 0 ...(1) 

Comparing (1) with Rr + Ss + Tt + f(x,y,z,p,q) = 0, here 

R = 1, 5 = 0, T = x 2 so that S 2 -4 RT = -4x 2 <0, x*0, showing that (1) is 
elliptic. 

The X -quadratic RX 2 + SX + T = 0 reduces to A 2 +x 2 =0 giving 
X = ix, -ix. 

The corresponding characteristic equation are given by 
— + ix = 0 and —~ix = 0 


2 \ 2 

X X 

Integrating, y+i — =c, and y-i — 

l 2 J 12 


u 2> l (x 2 ) 

Choose u = y + i — =a + i/?and v = y-i — = a-i/3, 
2 J 2 j 
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Where a 


are now two new independent variables 


Using (5) and (6) in (1) the required canonical, form.j,S 
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Heat, Wave and Laplace Equation 


CHAPTER-6 


HEAT, WAVE AND LAPLACE EQUATION 


Introduction 

In physical problems, we always seek a solution of the differential equations 
which satisfies some specified conditions known as the boundary conditions. 
The differential equation together with these boundary conditions, constitute 
a boundary value problem. 

In problems involving ordinary differential equations, we may first find the 
general solution and then determine the arbitrary constants fromthe initial 
values. But the same process is not applicable to problems Involving partial 
differential equations for the general solution of a partialdifferential equation 
contains arbitrary functions which are difficult to adjust so-as to satisfy the 
given boundary conditions. Most of th^ boundai^ walue'prqblems involving 
linear partial differential equations dan be solved bytlie following method. 

A solution which breaks up into a producLof functions each of which 
contains only one of the variables. The following qxplairis this method. 

Method of Separation of Variables \ s \ J 

Example: Solve by the method of separation of variables: 


d 2 z ^dz + dz 
8x 2 dx dy 




Suppose that the given equation has solutibn'tifthe form 
z(x,y) = X(x)Y(y) 

Where x is a function of * alone and r that of y alone. 
Substituting this value of z in the given equation, we have 

XT - 2X7 +AT = 0 Where X' =—,Y' = — etc. 

dx dy 


Separating the variables, we get 


X" - IX' 


Since x and y are independent variables, therefore, (ii) can only be true if 
each side is equal to the same constant, a (say). 


X" -2X 1 


= a , i.e. X"-2X'-aX = 0 


and— = a,i.e. Y'+aY = 0 
Y 


To solve the ordinary linear equation (m), the auxiliary equation is 
m 2 -2m-a=0, where m = l±^(l+a) 

/.The solution of (hj) is X = C| J i+ V( l+a )] j: + -J(\+a)] x 
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and the solution of (iv) is Y = c 3 e ay . 

Substituting these values of x and y in (/) we get 

i.e. z = 

Which is the required complete solution. 

Example: Using the method of separation of variables, solve — = 2 — + u 

dx 3t 

where u(x,0) = 6e' ix . 

Solution: Assume the solution u[x,t) = X(x)T(t) 

Substituting in the given equation, we have 
XT = 2XT'+XT or (X'-X)T = 2XT' 

X'-X T' 



V 

v ^ 

v- ' 


N>. x 


\ \ 


or 


- = — = k (say) 
2X T 


. 

• • V? 

fot' 

■ 

, ; • '} 


\ \ 

•> 

- 

\ 

l+jiitlc 


X'-X-2kX = 0Or— = \ + 2k 
X 

T' 

And — = k 
T 


Solving (/), log A = (l+2£)x+logc or 2( = c'e\ 

From (ii) log T = fo + logc'or T-c'e h 
Thus u(x,t) = XT = cc'e { ' +2k)x e kt 

Now 6e~ 3x = u (x,0) = cc'^ +2k ^ x 
cc' = 6 and 1 + 2k = -3 or k = -2 
Substituting these values in (iii ), we get 

u=6e~ ix e~ 2i i.e. u = 6e^ x+2t> which is the required solution. 

6.3. Some Important Equations 

1. Heat (Diffusion) Equation: V 2 n = | 

One-dimensional: — -t- = | — 

8x 2 ' 

~ . , d 2 u d 2 u 

dx 2 dy 2 


-(0 

...(h) 




(L) 


UJ 

[dt J 


m 

f 5u ) 

UJ 



(-) 

f du' 

UJ 

k dt , 


rpl 1 • • 1 V U V U U U 

Three-dimensional: —- + —- + —- = 

dx 2 dy 2 dz 2 


d 2 u d 2 u d 2 u n v 

Ua 


Note: Sometimes we write c 2 instead of k in heat equation. 


. C! ! y 


V- * . \ \ ‘ [ ■, 




‘ ■y-.Jv 






i 
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is of the form u(x,t)=X(x). T(t) 

Where x is a function of * alone and t is a function of t only. 

Substituting this in (1), we get 

xr=c 2 xr, i. Q .r'/x=r/c 2 T ...(2) 

Clearly the left side of (2) is a function of x only and the right side is a 
function of t alone. Since x and t are independent variables, (2) can hold 
good if each side is equal to a constant k (say). Then (2) leads to the 
ordinary differential equations. 

^--kX = 0 ,./3) 


and—-fc 2 r = 0 


Solving (3) and (4) we get 

(i) When k is positive and = p 2 , say: 

X = c l e px +c 2 e- pX ,T = c 3 e c2 P 2t ; 

(ii) When k is negative and =-p 2 , say: 


d l u 1 


Remark: 1. V 2 is called the Laplacian operator and is defined as 

V 2 =—t +—7 + — 7 - A function u which satisfies Laplace’s equation 
dx' dy' dz' 

is called the harmonic function. 


*1 F___ l 0 U 


1. TT V UUUV11U1V11U1 « ^1 a ^ 

8x 2 dy 2 U 


. d 2 u d 2 u d 2 u 
Three-dimensional: —- + —— + —- 

dx 2 dy 2 dz 2 


3. Laplace’s (or harmonic) Equation : V 2 « = 

Two-dimensional: = 0 

dx 2 dy 2 

ru j- • , d 2 u d 2 u d 2 u 

dx 2 By 2 dz 2 

Solution of the Heat Equation 


Assume that a solution of - 





















An ISO 9001 : 2008 Certified Institute 


X - c 4 cos px + c 5 sin px, T = c 6 e c p f ; 

(iii)When k is zero: 

X = c 7 x + c 8 , T = c 9 

Thus the various possible solution of the heat equation (1) are 


u-(c,e- ...(5) 
u = (c 4 cos px + c 5 sin px) c 6 e~ c p ' •••(6) 

u^x+c^O) ...(7) 

Of the three solutions, we have to choose that solution which is consistent 
with the physical nature of the problem. As we are dealing with problems on 
heat conduction, it must be a transient solution, i.e., « is todecreasewith the 
increase of time t . Accordingly, the solution given by (6), i.e. of the form • 


u = (C ( cos px + C 2 sin px) C 3 e e p * 


is the only suitable solution of the heat equation. \ ' 

Example: General solution of heat equation when botfisthe ends of a bar are 
kept at temperature zero and the initial temperature, i&jgreseribed. 

If both the ends of a bar of length a are attemperature zero and the initial 
temperature is to be prescribed function f(x) in the bar, then find the 
temperature at a subsequent time t. \ ^ 

or The faces jc = 0 and x = a of an Hnfihite slab are maintained at zero 
temperature. Given that the temperature' v ft-(-Jt, f) = f(x) at t = 0. Determine 
the temperature at a subsequent time t. 

Solution: Here the temperature u(x,t) in the given solid is governed by the 
... . ,(d 2 u'\ du 


one dimensional heat equation k\ —- = — ...(1) 

) dt 

Since the ends x = 0 and x = a are kept at zero temperature, the boundary 
conditions are 

w(0,t) = 0 and «(a,f)==0 ,Tot all t ...(2) 

The initial condition is given by w(x,0) = /(x),0<x<u . ...(3) 

Suppose that (1) has solutions of the form u (x,t) = X(x) T(t) ...(4) 

Where Xis a function of x alone and T that of t alone. 

Substituting this value of u in (1), we get 


k X"T=XT' o- 


x" _ r 

1 X ~kT 


Since xand t are independents variables, (5) can only be true if each side is 
equal to the same constant, say p . Hence (5) leads to the following 
equations: 

X"-fiX = 0 ...(6) 
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andT' = -pikT: ...(7) 

Using (2), (4) gives X (0) T(t) = 0 andA'(u) T(t) = 0 ...(8) 

•» 

Since T(t) = 0 leads to u = 0, so suppose that T(t) * 0. 

.-.Form(8), Af(0) = 0 and A r (a) = 0 ...(9) 

We now solve (6) under B.C. (9). Three cases arise. 

Case I. Let /u-0. Then solution of (6) is X(jc) = Ax + B ...(10) 

Using B.C.(9), (10)gives 0 = B, 0 = Aa + B so that A- B = 0. 

Hence = 0 so that u = 0, which does not satisfy (3). 

Case II. Let fi = A 2 , A ^ O.Then solution of (6) is 

X(x) = A^ JC +Be~ Xx 
Using B.C. (9), (11) gives 0 = A + B and 
0 = Ae aX +Be aX s 

Solving (12), A = B = 0 so that X(x) = 0 and 
Which does not satisfy (3). So we also rej 

Case III. Let ju - -A 2 ,A*0. Then soluti 
X(x) = AcosAx + BsmAx • 

UsingB.C.(9),(13)gives 0 -A and ( 

So sin/la = 0. we have taken 5*0, since otherwise X = 0so that u = 0 
Which does not satisfy (3). Solving the trigonometric equation sin Aa = 0, we 


Aa = nn so that A = —, where n = 1,2,. 

a 

Hence non-zero solutions X n (jc)of (6) are given by 


v ( \ d ' t nnx 

X n ( x ) = B n sin - 

a 


Using (14), (7) reduces to 


dT n 2 n 2 k 


dt, as jj. = -A 2 = — 


(i ■) (n 2 n 2 k) 

or 1 dT = -C]dt, as C„ 2 = - - 2 - 

\T) a 

Whose general solution is T n (f) = D n e c "‘ 


•v « n (x,t) = X n (x)T n (t) = E n Sin — e - c "\ 

V a 


axe solutions of (1), satisfying (2) . Here E n (= B n D n ) is another arbitrary 
constant. In order to obtain a solution also satisfying (3), we consider more 
general solution. 
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u{x,t)^u n {x,t) = ^E n sm 


nn jc^i 


Cit 


~r- 


.(18) 



n =1 


Substituting t = 0 in (18) and using (3), we get f{x) = Z £ » sin 
which is Fourier sine series. So the constants E n are given by 

2 f" 


(tinx'' 


V a J 


J-'- 

. 


E ■ =t1.7(*H 


f nnx'' 


dx, n- 1,2,3 


. 

< : v 


v “ J 


Hence (18) is the required solution where E n is given by (19). 

Working Rule For Solving Heat Equation When Both The Ends of A 
Bar of Length a are Kept At Temperature Zerp and The Initial 




■ 


Temperature /(x)is Prescribed. 


1 ' 


XX, \\ 

Step-1. Proceed as previous and prove that the solutionjof theheat equation 

. N v\\ " \ 

du . V-^ v >- 


^ d 2 u^ 


\ dx1 j 


dt 


\ K : \ \ 

subject to the boundary conditions u(0,t) = u (a,i] Mljjor all t 


\ 


- 0 ) 

...( 2 ) 


\ 


and the initial condition w(x,0) = f(x), 0 <{x <\...(3 j 


is given by u(x,t) = sin| \e c "‘ 

n=0 V a 


\ 


Where E„ =-£/(*) sin 


^ nnx^ 


dx, 


V “ J 


...(4) 

•••(5) 


and Cl - 




,.( 6 ) 


Step-II. Compare the given problem with (1), (2) and (3) and find particular 
values of k, a and /(*)• 

Step-III. Substitute the particular values of k, a and f(x) in (5) and (6) to 
get E n and C] 

Step IV. Substitute the values of coefficients E n and C 2 obtained in step III 
in (4) to arrive at the desired solution of the given boundary value problem. 


d 2 u (\ 


Example: Solve the one-dimensional diffusion equation —-= — 

dx k 


du 


dt 


m 


the range 0<x<2^, t>0 subject to the boundary conditions: 
M(x,0) = sin 3 x for 0<x<2;rand «(0.t) = «(2;r,t) = 0for t> 0 

Solution: Proceed upto equation (18) as in previous example by taking 
a -2n and f{x) = sin 3 x. Then, equation (18) for the present problem 
reduces to 


/ \ t" 2 * 1 ) 

' nx ' ——- 


n =1 


“(*’*) = IX sin — e 4 , n = 1,2,3,... 


V x J 


• 0 ) 


Putting t = 0in (i) and using given condition w(x,0) = sin 3 x, we get 




: •.= 


■ 

' -: ’ . 






. 7 : 


■ 




-■ . 

• •’ • •• 


■ 

‘life.; - ' 

45 * 1 * - :• . .. • 

v - V ' 




—* 



mcm 
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V E„ sin -— = sin 3 x = —sinx—sin3x 
S l 2 J 4 4 


.\sin3x = 3sinx-4sin 3 x=>sin 3 x = — x(3sinx-sin3x) 


(x) f3x) 5x 

Or^sin — +£ , 2 sinx-f-£ 3 sin — + £ 4 sin2x + £’ 5 sin — 
12 / 12 / 2 


(7x1 (31 (11 

+£ , 6 sin3x + £ 7 sin— — x sinx- — xsin3x ...(«) 

V2 ) {4J VV 

Equating the coefficients of the like terms on both sides of(n), we get 


E 2 =—, E 6 = — and E n = 0 when n* 2 or n* 6. Substituting these values 
in (/), the required solution is , 


(IV \ 

u ( x, t ) = E 2 sin xe~ b + E 6 sin 3 x e" 9h = - xpsin xe" K - sin ixe~ 9h ]. 

v 4 / ,\V. 

\ \\ 

General Solution of Heat Equation When Both The Ends of A Bar are 
Insulated and The Initial Temperature is Prescribed. 

If both the ends of a bar of length a are insulated and the initial temperature 
, \ \ ' 

/(x) is prescribed, then to find the temperature at a subsequent time /. 

v ' V N. \ 

\ 

Solution: Here the temperature u(x,t) in .the given bar is governed by one 
dimensional heat equation 

...a) 

^ dx 2 J dt 

Physical experiments show that the rate of heat flow is proportional to the 
dll 

gradient — of the temperature w(x,/). Hence if the ends x = 0 and x = a of 
the bar are insulated, so that no heat can flow through the ends, we have 
the boundary conditions: u x (0,/) = u x (aj) = 0 for all t ...(2) 

Also, the prescribed initial condition => «(x,0) = /(x) for all x ...(3) 

Suppose that (1) has solutions of the form u[x,t)-X[x)T(t), ...(4) 

Where X is a function of x alone and T that of t alone. 

Substituting this value of u in (1), we get 


kX"T = XT' or 
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Since x and t are independent variables, (5)can be true if each side is equal 
to the same constant, say fi. Thus, (5)yields. X" - nX = 0 ...(6) 

and V = fikT ...(7) 

Differentiating (4) partially w. r. t. V, we get u x (x,t) = X'(x)T{t) ...(8) 

Using (2), (8) gives r(0)r(t) = 0 and X'(a)T(t) = 0 ...(9) 

SinceT(t) = 0 leads to u- 0,so assume that T(t)&0. Hence, from (9), we 


JT(0) = 0 and X'(a) = 0 
Three cases arise: 

Case I. Let // = 0. Then solution of (6) is X (x) = 

Which yields X' ( x) = A \ ^ 

UsingB.C.(10), (12)gives A-0 .Then (11) 

Again, corresponding to /i = 0, (7) yields 

dT E 

— = 0 So that T =constant = —, say \ 
dt 2B n 

v \ 

Corresponding to /i = 0a solution \pf the given boundary value problem 
from (4) is given by \ 

u(x,t) = Bx —1 = — ...(13) 

y ’ \2B) 2 

Case II. Let |i = X 2 , X * 0 Then solution of (6) is 
X(x) = Ae Xx + Be kx ...(14) 

Which yields X' (x) = Ale** - BAe^ ...(15) 

Using B.C. (10), (15) gives 0 - A(A- 5) and 0 = A (Ae* - Be *). 

These equations yield A = B = 0 so that X[x) = 0 and hence u = 0, 

Which does not satisfy (3). So we reject //=/1 2 . 

Case III. Let // = -/l 2 , A^O. Then solution of (6) is 

X[x)-AcosAx + BsmAx ...(16) 

Which yieldsX'(x) = -^/lsin>lx+5^cosAx ...(16)' 


Using B.C. (10), (16)' gives 0 = 5/1 and 0 = -A As\nAa +BAcosAa. Thus, 
5 = 0 and A A sin A- 0. Since A*- 0 and A *0 for a non-trivial solution, we 
must have 

sin/la = 0 So that Aa-nn 
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Thus, A = —, n = 1,2,3,... 
a 


Hence non-zero solutions X n (x)of (6) are given by 


f riKX^ 


X n (x) = A n cos - , by (16) and (17) 

v a J 


Using (17), (7) reduces to 
dT n 2 n 2 k 


...(17) 


...(18) 



T a 

rn 


,2 n 2 n 2 
dt, as n - -A -- 


or 


V* 7 


Where C\ = 


dT = -C 2 dt 
( n 2 x 2 k ) 




...(19) 

...'(20). 
...( 21 ) 
...( 22 ) 

are solutions of (1). For n - 1,2,3,..., each^one these satisfy the boundary 
conditions (2). Here E n (= A n D n ) in another arbitrary constant. 

\ \ 

Thus, (13) and (23) constitute a set of infiriit@»,solutions of (1). To obtain a 

*. . 1 \ ****** 

solution also satisfying the initial condition (3), we consider a linear 
combination of these solutions. Hence the-^omplete solution of (1) may be 
taken in the following form. 


\ 


Solving (19), T n (t) = D n e<‘ 

• u„(x,0 = X„(X)r„(t) = £„cos[ |e“ c "'^'by (18)’and"(21) 


u(x,t) = ^- + ^u n (x,t)or u{x,t) = ~ + Y l E n cos—e" c "' ...(23) 

} 1 a 


n =1 


«=1 


Substituting t = 0 in (23) and using (3) , we have 

...(24) 


/w4+i*«—. 

^ n =1 Q, 


Which is Fourier cosine series. So the constants E 0 and E„ are given by 

E 0 =-\°f(x)dx and E n =-[V(x)cos— dx, n = 1,2,3,.(25) 

a Jo a Jo w a 

Hence (24) is the required solution where E 0 and E n are given by (25) 

Working Rule For Solving Heat Equation When Both The Ends of A 
Bar of Length a are Insulated and The Initial Temperature f(x) is 
Prescribed. 

Step-1. Solution of the heat equation 


^ d 2 u^ 
K dx 2 j 


du 

dt 


Subject to the boundary conditions u x (0,f) = u x ( a,t ) = 0 for all i 
and the initial condition «(x,0) = f(x), for all x 


...( 1 ) 

...( 2 ) 

•••(3) 


■ 

■ ' • 








r ' v *; : * r ' 









! ■ 
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is given by u{x,t) = — + cos^^-e c "' 

2 n =i a 




2r* , x , 2 , s Mix 

Where £ 0 = — / [x)dx, E n =-\ f(x)co$ - dx, n = 1,2,3,.(5) 

fl J0 a ^o a 

(«V 2 ^) 

andC " = a 2 — -( 6 ) 

Step-II. Compare the given problem with (1), (2) and (3) and find particular 
values of A, a and / (x). 

Step-III. Substitute the particular values of A', a and /(x) in (5) and (6) 

and calculate £ 0 ,£,and C n . ''“x s \ 2 • 

Step-IV. Substitute the values of coefficients E 0 ,E n and f obtained is step 

III in (4) to arrive at the desired solution of the given'boundary value x 2 

problem. \ . raggl ^ / 

* \j- £ V ?'-f f: f f % f 

(d^u\ du \ \ *' 

Example: Solve A —j =—for 0 <x < n,t >u,iff 0,/) = if/r,f) = 0 

. Ox dt \ x. •. ■ ' .."■ \ . 


and «(x,0) = sinx. 


Solution: We can prove that the solution of heat equation 

\ V— 

fd 2 u^_du , ' 

obf dt 

Subject to the boundary conditions u x (0 ,t)=u"[a,t) = 0for all t 
and the initial condition u(x,t) = /(x), 0<x<a 

• • , / \ -E’o V 1 r, nKX -ch 

is givenby M(x,t) = —+ ^£ n cos- e ■ 

2 n =i a 

Where £ 0 = —f f(x)dx, E n - — \ f(x)cos^-dx,n = \,2,3... 
a Jo a Jo a 


•( 2 ) 


and C: = 


(n 2 ;r 2 A) 


...(6) . 


Comparing the given boundary value problem with the boundary value 
problem given by (1),(2) and (3), we have k = k, a-n and /(xfsinx. 
Hence, from (5), we get 


E n - — f sinxclx =—[-cosxF =—x2 = —...(7) 

u 77- Jo v L J0 V V 


n n 


and E n = — J'*sinxcosrcxd[r = — |^|sin(n + l)x-sin(n-l)x| 


; ■ 


if cos(« + l)x cos(«-l)x T 1 (-if (-if l, 1 . 

n k + 1 n-1 n (n + l) n-1 n + l n -1 
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If n = 2m-\ (odd) with m = 1,2,3,...., then E n =0,by (8) 
If n = 2m (even) with m = 1,2,3,..., then from(8), we get 


Substituting the values of E 0 and E 2m given by (7) and (9)'in (4) 
required solution is nT''— 


6.5. Solution of The Wave Equation 


Assume that a solution of 


y(x, t)-A r (x)r(t) where 
only. 


x is a function of-* and r is a function of 


Substituting these in (1), we get XT'-cXT i.e 


Clearly the left side of (2) is a function of * only and the right side is a 
function of t only. Since * and t are independent variables (2) can hold 
good if each side is equal to a constant k (say). Then (2) leads to the 
ordinary differential equations: 


Solving (3) and (4), we get 

(i) When k is positive an - p 2 , say X = c x e px +c 2 e~ px \ I 


(ii) When k is negative and 
X = c 5 cos px + c 6 sin px\ T = c 7 cos cpt + c 8 sin cpt. 

(iii) When £ is zero X = c g x+c l0 \ T=c n t+c n . 

Thus the various possible solutions of wave equation (1) are 
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+c A e- c P‘) 

y=(c 5 cos px+c 6 sin px)(c 2 cos cpt+qsmcpt) 
y={ c 9 x+c \o){ c u t+c n) 

Of these three solutions, we have to choose that solution which is consistent 
with the physical nature of the problem. As we will be dealing with problems 
on vibrations, y must be a periodic function of x and t . Hence their solution 
must involve trigonometric terms. According the solution given by (6), i.e. of 
the form 

y = ( Cj cos px+Cj s in px) (Cj cos cpt + C 4 s in cp/) ...(8) 

is the only suitable solution of the wave equation? \ \0' 

\ v 

Example: A sting is stretched and fastened to two points / nj>art. Motion is 
started by displacing the string in the form y = a sum ray/ I'from which it is- 
released at time t = 0. Show that the displacement of any point at a distance 



x from one end at time t is given by 
y(x, f)=asin(ra: / /) cos(nct //) 


\\ ' 
\ 

\\ 


Solution: The vibration of the string is given by -^' = c 2 —\ 

\ dt'"~ dx 2 

1 NN*., 


...(/) 


V> 


As the end points of the string are fixed, for all time, 

y(o,/)=o v 

and y (/,*) = 0 

Since the initial transverse velocity of any point of the string is zero, 
Therefore, f — I =0 

Also y(x,0) = asin(px//) 

Now we have to solve (/) subject to the boundary conditions (ii) and (iii) 
and initial conditions (iv) and (v) . Since the vibration of the string is 
periodic, therefore, the solution of (i) is of the form 

y(x,t) = (C] cos px+Cj sin px)(Cj coscpt+C 4 sin cpr) 

By (ii ), y(0,r) = C 1 (C^coscpt+C 4 sintpf) = 0 

For this to be true for all time, Cj = 0 


-(«) 

..(Hi) 

-H 

...(v) 


...(W) 


-H 


and ~- = C 2 sin px{C 3 (-cp. sin cpr)+C 4 (cp.coscpt)} 


By (/v), 1 — 1 = C 2 sinpx.(C 4 cp) = 0, whence C 2 C 4 cpsin px = 0 

V ot J l=0 

If C 2 = 0, (vii) will lead to the trivial solution y(x,f)=0 


o--*. r_ - 

. • ■ 


>- ; 
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The only possibility is that C 4 - 0 

Thus (vii) becomes y(x,t)=C 2 C i smpxcoscpt 
By (iii), y(/,/) = QQsin/)/cosc^=0 for all t. 
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...{viii) ■ 

X 7 - 


Since C 2 and C 3 * 0, we have sin pi = 0. /. pi = n/r , i.e. p = nnl /, where 
n is an integer. •: 

Hence (/) reduces to y(;c,/) = C 2 C 3 sin^y^cos-;-/k ■C"X k' k^ 

[These are the solutions of (/) satisfying the boundary conditions. These 
functions are called the eigen function corresponding td' the eigen values 
l n =nn/l of the vibrating string. The set of values Is called is 

spectrum]. \ - -X* * • . . * 

Finally, imposing the last condition . (v), we have 

flTUC 71X \ \ ^ 

_y(at,0) = C 2 C 3 sin — = a sin — which will be satisfied by tanking C 2 C 3 = a 

and n-\ \ \'\ . 

\ \\ . 

Hence the required solution is y(x,t) = asin.^cos—S- 

Working Rule for Solving One-Dimensional Wave Equation Satisfying 
the Given Boundary and Initial Conditions 

\ ""N ' d 2 v ( n fd 2 y'] 

Step-I. The given solution of the wave equatidn —t = \~t x —t —(l) 

dx v c J dt J ■' 

subject to the boundary conditions y(0,t) = y(a,t) = 0 for all t ...(2) 

and the initial conditions y(x,0) = /(x), y, (x,0) = g(x) 0<x<a ...(3) 

. • , / > ff„ rmct _ . nnct\ . nxx ... 

is given by yyx,t j = ^ E n cos-+ F n sm- sin- ...(4) 

n=iV a a ) a y- 

where £■„= — [ f(x)sin^^-dx and F n --^— f g(x)sin-^^-Jx ...(5) 
a Jo a nnc Jo a 

Step-II. Compare the given problem with (1), (2) and (3) and find particular 

values of c, a, /(x) and g(x). 

Step-III. Substitute the particular values of a,/(x) and g(x) in (5) and 
compute E^miF,,. 

Step-IV. Substitute the values of E n and F n in (4) to arrive at the desired 
solution of the given boundary value problem. ‘ 

Working rule for solving One-Dimensional Wave equation when both V 

the ends of the string of length / are fixed and they initial velocity of the 
string is zero. 

First of all, as explained in working rule, prove that the solution of one- . . a. . '. V" • 

d 2 f i y a 2 \ 

dimensional wave equation—j = “j ~zr 

ox \c jya ) 
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subject to boundary conditions y(0,r) = y(/,0) = 0 for all t 


initial velocity = — = 0, 0 < x < / 

and V dt J l=0 ^ ...(Hi) 

initial deflection =y(x,0) = /(x),0 < x < / 
is y(x,t) = 2_, E n sin—cos—— — (*v) 

n=l ‘ ' 

Where E n = j^f(x)sm^-dx. 

Example: A string of length / has its ends x = 0 and x = l fixed. It is 

4Xjc(1-x) x X \ 

released from rest in the position y =-^-. Find an expression for the 

L \ •% v 

displacement of the string at any subsequent time. . \ \ 

The displacement function y(x,t) is the solution of theyyave equation 


subject to boundary conditions: y(0,<) = y(l,i) = 0 'for all t > 0. 

\ *\ 

and initial conditions, namely, ~ v 

Initial velocity=[ — | =0 for 0<x <I: 

UJ,.„ \ 


...3(a) 


and initial displacement = y(x,0) = /(x) ■ 


(4;rx(/-x)} 


Comparing the given problem given by (1),(2), 3(a) and 3(b) with the 
boundary value problem given by (;), (ii) and ( iii ), the solution of (1) 
satisfying the above boundary and initial conditions is given by 


/ n ~ . nnx nnct , ,, . 

b y ( lv ) 

n —1 * / 


Where E n = y£/(x)sin-^y^a!x, by (v) 

Substituting the value of /(x) given by 3(a) in (5), we get 


(2) ( 4A N )fi / ' 2 ^ • nxx j 

E = — x — Ix-x sin- ax. 

I l 2 J JoV ' l 


81 /. 2 \( / nnx\ /, . J l 2 . tiKx'\ , „n ( l 3 tmx 

=-r- Zx-x-cos-(/- 2x) —^-rSin- +(-2 -^rcos- 

/ 3 [ V \ tm l ) y \ nV / J V \nV l 

[Using the chain rule of integrating by parts] 
f 0, if n = 2m (even) with m = 1,2,3,... 

1(321) / (2m -1) 3 ^r 3 ,if n = 2m - l(odd) with m = 1,2,3,... 
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l 


Substituting the value of E n in (4), the required expression is 

x 32AA 1 . (2m-lW;t (2m-l)zct 

K ’ ^ S(2/n-l) 3 / / 

Solution of Laplace’s Equation 

d 2 u d 2 u 
—r +—r = ° 
dx 2 dy 2 

Let u(x,y) = X(x)Y(y) be a solution of (1). 

2 2 

Substituting it in (1), we get Y + X XL = 0 V 

<fc 2 dy 2 < 'X 






...(!) x 


or separating the variables , 


1 4 2 X 1 af 2 y 

* 4 x 2 _ ? 4 y 2 


Since * and y are independent variables, (2) can hold-good dhly if each side 
of (2) is equal to a constant k (say). Then (2)"leads to the drdinary differential 
equations \ X^X" 

V \ ^ 


^--kX = 0 and^L + *f = 0. 

dx 2 dy 2 


Solving these equations, we get 


(i) When k is positive .andx^s equal to p 2 ,say 


v 



X = c x e p +c 2 e p , Y - c 3 cos py + c 4 s^^ - 

(ii) When k is negative, and is equal to -p 2 , say 

X = c 5 cos px+c 6 sin px, Y = c 7 e py + c%e~ py 

(iii) When k is zero; X = c g x + c w , Y = c n y + c n . 

Thus’ the various possible solutions of (1) are xg 
u - [c x e px + c 2 e~ px j (c 3 cos py + c 4 sin py) ...(3) 

u = (c 5 cos px + c 6 sin px) [c 7 e py + c % e~ py j ...(4) 

w = (c 9 v+c 10 )(c 11 y+c 12 ) ...(5) |JX 

Of these we take that solution which is consistent with the given 
boundary conditions. XX 


1 6.6.1 General Solution of Two-Dimensional Laplace’s Equation 


d 2 u d 2 u 


Consider X + X = 0. 
dx 2 dy 2 

Suppose that (1) has solutions of the form u(x,y) = X[x)Y(y) 
Where X and Y are functions of x and y respectively. 


From (2) ^- = X'Y and ^L = XY". 
dx 2 dy 2 


(?) ' 


X-■ >X'- 


•••(3) fXv 

v/xX XX-. .•*.*•••• v.- v.- •;•*. -r '•/ 

&?X :.r 

lXX' : • .. 
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Since the L.H.S. of (3) depends only on xand the R.H.S depends only ony, 
both sides of (3) must be equal to same constant, say /U. This leads to two 
ordinary differential equations 

X” - nX = 0 and Y" + /iY = 0, ...(4) 

whose solutions depends only on the value of ju . Three cases arise: 

Case I. When // = 0. Then (4) reduces to X" = 0 and Y" = 0. 

Solving these X = A } x +5, and Y = C l y+D i 

Then a solution of (1) is «(x,y) =(T,x+5,)(C 1 y+Z) 1 ). ...(5) 

Case II. When ju = X 2 , i.e., positive. Here A * 0. Then (4) reduces to 

\ \ 

X"-X 2 X = 0 aad Y"+A 2 Y = 0. v X^X 

Solving these, X = A 2 e ix + B,e~ ls and Y = C 2 cos Xy + D 2 sin A)C\'~ 
or, alternatively, A" = .4, cosh Xx+B 2 sinh Xx, Y = C 2 cos Xy + D 2 sin Xy 
Then a solution of (1) is u{x, y) = ) (C 2 cqs XyA-D 2 sin Ay) ...(6) 

or alternatively, u(x,y)=(A 2 cosh Ax+5, sinh Xx) (C. cos Ay+ D 2 sinAy) ...(6)' 

Case III. When n = -A 2 , i.e., negative. Here A * 0;Then (4) reduces to 

V \ \ 

X"+X 2 X = 0 and Y"-X 2 Y = 0. , \ 

Solving these, X = A } cos Ax + B } sin Ax and Y =-C 2 e h + ■ 

N. v- 

or alternatively, X = A, cosAx + 5, sin Axand K=t) 3 cosh Ay+ D 3 sinh Ay 
Then a solution of (1) is \,^ 

u(x,y) = (A } cosAx+5 3 sinAx)(^C 3 e i> +D } e' Zy j. ...(6) 

or alternatively, «(x,y)=(^ cosAx+5 3 sinAx)(C 3 coshAy+D 3 sinh Ay) ...(7)' 

Out of the above mentioned three types of solutions (5), (6) and (7), we must 
select an appropriate solution which suits the physical nature of the problem 
and given boundary conditions. 

.6.2 Dirichlet Problem in A Rectangle 

1 st Case . The Dirichlet problem in a rectangle is defined as follows: 
d^u d^u 

Laplace’s equation: —- +—- = 0, 0<x<a, 0 <y<b ...(1) 

dx l dy 

Boundary conditions: n(0,y)=0, w(a, y) = 0 0<y<b ...(2) 

u(x,b)= 0, 0 <x<a ...3(a)andH(x,0)= /(x), 0<x<a ...3(b) 

We shall solve (1) subject to the boundary conditions (2), 3(a) and 3(b). 
Suppose (1) has a solution of the form 

u(x,y)=X(x)Y(y) ...(4) 

Substituting this value of u in (1), we get 


. rX \ ... -'X?-'‘>; 


X"Y+XY" = 0 or—= -—. 

X Y 

Since x and y are independent, variables each side of 
same constant, say fi. 


* X ‘- 1 &? -*4 v x‘<V.- : 

t X-X- XX X X vj & 

'•'•-V V . - i. '■ _ ■ ,• ;■ 


(5) must be equal to the 

. 
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Then (5) gives X"-fiX = 0 
andF"+//T = 0 
Using (2), (4) gives X(0)T(y) = 0 and X(a)T(y) = 0 
giving X(0)=0 and X(a)=0 ...(8) 

Where we have taken T(y)^0, since otherwise w = 0 which does not satisfy 
3(b). 

We now solve (6) under B.C. (8). Three cases arise: 

Case I. Let // = 0. Then solution of (6) is X(x} = Ax+B ...(9) 

Using B.C. (8),(9) gives 0 = B and 0 = Aa+B. These give A=JB= 0 so that 
X[x) =0. This leads to u = 0 which does not satisfy 3(6)>.So wfgreject /u = 0. 

Case II. Let. n = A 2 , A* O.Then solution of (6) is X’(x) -.Ae x {+Be~ xX ..'.(10)^ 
UsingB.C. (8),(l0)gives 0 -A+B and 0 = Ae aX A-Be~ aX '...(11) 
Solving (11), ^ = 5 = 0 so that X (x)-= 0 and hence - # = 0. So we reject fj = x\ 
Case III. Let /j = -A 1 , X * 0. Solution of (6) is N \\ ~ 

X(jc) = ^lcoS/lx+ J 8sin/lx ...(12) 

\ \ ' 

Using B.C. (8), (12) gives 0 -A and 0^4cb^^a + 5sin^ so that 
A = 0 and sin/la=0, , \. 

Where we have taken B* 0, since otherwise A r £x)=0and hence u = 0 which 
does not satisfy 3(b ). 



Now, sinAa = 0 => Xa-riK-. 


. nn t 0 . 

X = —, w = 1,2,3,... 
a 


Hence non-zero solution X n (x) of (6) are given by 


X ni x ) = B n sin 


nnx 




...(13) 

...(14) 

T = 0 ...(15) 

V J 

■nny 

Whose solution is Y ri (y) = C n e a +D n e a ...(16) 

Using3(a), (4) gives 0 = X(x)7(6) so that Y(b) = 0, where we have taken 
A"( jc) * 0, since otherwise u = 0 which does not satisfy 3(6). 

Again, T(6) = 0 => Y n (b) = 0 ...(17) 

Putting y = b in (16) and using (17), we have 


2 2 

XT- ,2 “HI 

Using p = -A. = —-— (7) reduces to Y" - 
a 1 


0 = C n e nKh/a + D n e~ nnh/a so that C n = - 


-\D.e 


...(18) 


Putting this value of C n given by (18) in (16), we get 
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Uy)=D n 


-nny nnb nny -nnb 

e a e “ -e a e “ 


~ or Y„{y) = D n 


e a -e a 




Y n (y) = 2 D n e a sinh 


:{b-y) 


as e -e - 2sinh# 


••• u n (x, y) = X n (x) Y n = F n sin (—1 sinh{ ^ 

a ) a 


are solutions of (1), satisfying (2) and 3(a). Here F n =2B n D n e a are new 


arbitrary constants. In order to obtain a solution also satisfying 3(/(), we 
consider more general solution, namely, V '\.. 

u (*> y) = Z(•*>J 7 ) = Z F n sin [~\ sinh { —~ a "^9) 


Putting y - 0 in (20) and using 3(b) , we have 


/(-) = £ F..mh(^ U = 

«=i \ a )\ V a 


W 


Which is the half range Fourier sine series of/(xjin (0,a). Hence, we get 


r, • , nnb 2 <•» , , . nnx 

rsinh-= — /(x)sin- ax 

a a Jo W a 


or F =- 


nnb ] Jo 


t a ,/ x . nnx 

J 0 /(*) sin - dx. 


\ V 

\V 


a 


...( 20 ) 


Example: Find the steady state temperature distribution in a rectangular plate 
of sides a and b insulated at the lateral surface and satisfying the boundary 
conditions 

u(0,y) = u(a,y) = 0 for 0 <y<b 

and u(x,b) = 0 and u{x,0) = x(a-x ) for 0 <x<a 

Solution: First proceed as in above article upto equation (20) 

For the present problem, w(x,0) = f(x) = x[a-x ), 0 <x<a 


nnb V° 


t a r/ \ • nnx . 

f (x sin- dx-- 

JO J W 


f fl / 2 \ • , 

(ax-jc jsin- dx 


nnb}* 0 


...( 21 ) 


/ 2\ a nnx / ^ \ 

[ax-x -cos- [a-2xj — 

' ' nrr n ' ' 


wry a 


2 \ / 3 \ 

a . wrx . x a wrx 

-o SU1 — + r 2 33 008 — 

nn a \nn a 


2_ 2 <* 3 (-If 2a 3 _ 4a 3 r . , 

3_3 ^ 3_3 ~ 3 __3 I ( U 


[Using chain rule of integration by parts] 


cosech- 


35 S Y4'$±-%\ Yi d 
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0, ifn=2m (even) andro=l,2,3,... 

jsa 2 / 7? (2m-l)Hcosech {(2m-l);r6/a], if« = 2m-l and m = 1,2,3,... 




' 

' 


Substituting the above value of F in (21), the required steady temperature 
«(^, 7 )is given by 

v ' “V:-; 

, , 8fl ! v- 1 . (2m-l)xx . , (2m-l)(b-y)x Jlm-\}xb 

u(x,y )=—V -r-sin--—sinh--—-— coseclU-— 

1 ^ £(2m-l) 3 « a a 

2 nd Case 

The steady state temperature distribution in a rectangular plate of sides a and 
b insulated at the lateral surface and satisfying the boundary conditions. 


u (0, y ) = u (a, y ) = 0 for 0 < y < b and u [x, 0) = 0 and u(x,b^- /(x)for 

v .. v ''•V ■ \ 

0 < x < a x \ ^ 

Solution: The heat flow in a body for a two-dimensional case is governed by 

d 2 u d 2 u f 1 Ydw'l “ s t ,\ \ , .. . 

the equation —r + — T = — — . When the steady sate condition 

ax- 2 dy 2 (c 2 JldtJ \ \\ 

prevails, temperature u is independent of t so that — = 0. Then, the above 

\ dt 

equation reduces to Laplace’s equation Y"—- 

fr+fr - 0 '•.'"'"0'" -d) 

dx 2 dy \ v "-- 

Given that w(0,j) = M(a,j) = 0 for0<}><6 ...(2) 

Also, given that w(x,0) = 0 for 0<x<a ...3(a) 

and u(x,6) = /(x) for 0<x<a ...3(6) 

Suppose (1) has a solution of the form 

u(x,y) = X(x)Y(y) ...(4) \t£ 

Substituting this value of u in (1), we get 


ytt ytr 

X"Y + XY" = 0 or — = -— 


•(5) 


Since x and t are independent variables, each side of (5) must be equal to the ifjd.rFi; 
constant, say ju . Then (5) gives X” - /iX = 0 -(6) ■ 7% . =; 

and Y" + yY = 0 -(7) 

Using (2),(4) gives A(0)f(y) = 0 and X[a)X(y) = 0 

giving X(0) = 0and X(a) = 0, ...(8) 

Where we have taken Y(y) * 0, since otherwise u = 0 which does not satisfy ^ ;, ; - 

m SgBftiil 

We now solve (6) under B.C. (8). Three cases arise: ]?fffY-YY 

Case I. Let /u-0. solution of (6) is Ar(x) = Ax+B ...(9) .i 




■ . *-• 
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Using B.C. (8),(9) gives 0 = 5 and 0 = Aa + B. these give A = B = 0so that 
X (x) = 0. This leads to u = 0 which does not satisfy (3 )b. So we reject 

>o. 

2 

Case II. Let p. = X ,X *OThen solution of (6) is 

X(x) = Ae Zx +Be~ Xx ...(10) 

Using B.C. (8), (10) gives A + B = 0 and 0 = e aX + Be aX ...(11) 

Solving (11), A = B = 0 so that Af(x) = 0 and hence u = 0 . 

So we reject // = A 2 . 

Case III. Let /J.- -A 2 , A*0. Then solution of (6) is V ' S \ 

X (x) = A cos Ax + B sin Ax V'\. ...(12) 

Using B.C. (8),(12) gives 0 = A and 0 = A cosAa + B sinlx/ 
so that A -0 and sin Aa = 0 x .. 

Where we have taken 5*0, since otherwise A'(x) = 0 and hence u = 0 
which does not satisfy 3(h). 

Now, sin/la = Ogives Xa = nn, n = 1,2,3 \ 


:.A=—, n = 1,2,3,... - \ 

a \ 

Hence non-zero solutions X n (x)of (6) are \ 


given by X n (x) = 5„ sin 

\ a 


Using // = -A 2 =—^-—, (7)becomes Y" -— 7 = 0 ...(15) 

a a J 

nny nny 

Whose general solution is Y n (y) = C n e “ +D n e a ...(16) 

Using 3(a), (4) gives 0 = AT(x) T(0)so that T(0) = 0, where we have taken 
X (x) * 0, for otherwise we shall get u = 0 which does not satisfy 3(h). 

But 7(0) = 0 

=> Y n {y )=0 ...(17) 

Putting y = 0 in (16) and using (17), we have 
0 = C n + D n so that D n = -C n . 

Then (6) reduces to 

( nny nny \ 


Y n (y)-C n e a -e a =2sinh 


as e e -e 6 = 2sinh5. :j ...(18) 

.-. M n (x,y) = X„(x)7„(} ; ) = £„sinf—Isinhf^l !j -09) 

\ a ) \ a ) | 

are solutions of (1), satisfying (2) and 3(a), we now considet more general 
solution given by |! 

;K 

;i 
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i e.A x ’y) = &-■■ ( x >y) =Z E - sin—sinh— 

n-\ n=\ a Q 

putting y = b in (20) and using 3(b ), we get 

, V’fc • • n7tX 

/W = g^.s,nh—jsm—, 

^bich is the half range Fourier sine series of f(x) in (0,a). 
Hence, we have 

, n/rb 2r* f , . nnb 
£sinh- r = -J 0 /Wsm —dx 


PPIP 

yaurownyot* 




a a 


rmb | Jo 


t'/Wsin^* 




Hence, (20) is the required solution wherein £ n is given by. (21). ' 

3 rd Case , \"x, N \\, 

d^u \. * 

The Laplace’s equation —- + — r = 0 subject to the following boundary 

11 dx dy \ \\ ‘.'i 

conditions: u(x,0) = u(x,b) = 0 for 0<x<a. n(0,y) = 0 and 

«(«,>’) = /(>’) for 0 < y < b. \ VX | 

Solution: Given Laplace’s equation is = o. ...(1) 

43 

Also given w(x,0) = n(x,h) = 0 for 0 <\< a,'...(2) 
u(0,>') =0for °^y^b. ...3(a) 

and u(a,y) = f(y) for 0<y<b. ...3(b) 

Let a solution of (1) be of the form u(x,y) = X(x)Y(y). ...(4) 

Using (4), (1) educes to Xy 

r . r + A -r = 0or (j)v = -[y]r'....(5) 

gince the L.H.S. of (5) depends only on x and the R.H.S. depends only on y, 
each side of (5) must be equal to the same constant, say fJ.. Then, (5) leads to 

= Q -(6) yX 

and Y" + \iY = 0 ...(7) 

Using (2), (4) gives 2f(x)T(0) = 0 and X(x)7(h) = 0 

So that T(0) = 0 and Y(b) = 0 ...(8) ' /: 

where we have taken X (x) t- 0, since otherwise u = 0 does not satisfy (3b). 

We now solve (7) under boundary conditions (8). Three cases arise: 

Case L Let /i = 0. Then, solution of (7) is Y(y) = Ay + 5. ...(9) 

Using B.C. (8), (9) gives 0 = B and 0 = Ab + B . These give A = B = 0 so 
that Y(y) s 0.This leads to u = 0 which does not satisfy 3(b). So we reject 

O X-3-.* >?:•* 

r _ 
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Case H. Let = (negative). Here A * 0. Then solution of (7) is 
Y(y) = Ae yX +Be yX ...(10) 

Using B.C. (8), (10) gives 0 = A + B and 0 = Ae bX + Be bX . ...(11) 

(11) => A = B - 0 so that X (x) = 0. This leads to u = 0 which does not 
satisfy 3(b). So we reject /u--A 2 . 

Case III. Let /t= A 2 (positive). Here 2*0. Solution of (7) is 
Y[y)-AcosAy + BsmAy. •••(12) 

Using B.C. (8), (12) gives 0 -A and 0 = ,4 cos 26+ 5 sin 26 so that ,4 = 0 
and sin 26 = 0, where we have taken 5* 0, since otherwise we shall get 
Y[y) = 0. This leads to u = 0 which does not satisfy 3(6),, . 

W \ 

YIK \ X \ X 

Now, sin Ab = 0 => Ab = nn => A = —, n = 1,2,3,„. 

b x \ 

Hence, non-zero solutions Y n (y) of (7) are. giveii by- \^s 



F„(y) = 5„sin —^ , n = 1,2,3,... 

V U ) 


\ \ \ 
\ ' N . 
X 


2 2 \ ’ 
Also, then p = 2 2 = —— hence, (6) reduces'to 
b 2 \ 


X = 0, n = 1,2,3,.(14r\ V 


nttx n,?x~" 


Whose solutions are X n (x) = C n e b +D n e b . ...(15) 

Using 3(a), (4) gives 0 = X(0)F(y)so that X(0) = 0, where we have taken 
F(y) * 0, since otherwise u = 0 which does not satisfy 3(6). 

Now X(0) = 0 =>2f n (0) = 0. Putting x = 0 in (15) and using X n (0) = 0, 
we get 

0 = C n + D n so that D n = -C„. Then (15) reduces to 


X n (x) = C„ |^e -e b J = 2C„sinh|^J, as e d -e 6 = 2sinh0 ...(16) 

From (4), (13) and (16), we see that a solution u n (x,y)of (1) is given by 

u n ( x >y) = x n { x ) Y n (f) = E „ sil ^f^r] sin f n = 1,2,3,.(17) 

b J 


Where E n {-B n C n ) are new arbitrary constants. In order to obtain a.solution 
also satisfying B.C. (36), we consider more general solution 

“(*>f) = IX {x,y) = Y, E n sinh^sin^ ...(18) 

n =1 / 2=1 O U 
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Putting x = a in (18) and using B.C. (3b), we have 

f{y) = Z sinh^sin^ 

«=l b b 

which is the half range Fourier sine series of / (y) in (0,6). Hence 
_ . ,nna 2 rb , . . nny , 

£ » smh X = i)o fW sm T>' 



o> b,=- 


V° 


j 0 V W sA 


Hence, (18) is the required solution wherein E n is given by (19) 


4 th Case 


d 2 u d 2 u 


The laplace’s equation —r- + —r- = 0 subject to the following boundary 
dx l dy l \ \ \ v 

condition n(x,0) = n(x,6) = 0 for 0 < x < a ' - ... 

\ \ \ 

u(0,y) = f(y) and w(n,y) = 0 for 0<y<b \ 

Solution: Given laplace’s equation is — %+—~:=-Q- 

Also given n (x, 0) = u (x, b) = 0 for 0^ x £%, ' 

“(°,y) = f(y) for 0<y<6 ...: 

and u[a,y)-0 for 0<y<6 ...2 

Let a solution of (1) be of the form 
u(x,y) = X(x)Y(y) 

Using (4), (1) reduces to 


...3(6) 


XY + AT = 0 or 


x" r 


Since the L.H.S of (5) depends only on x and the R.H.S depends only y, 
each side of (5) must be equal to the same constant, say p. Then, (5) leads to 

X"-yX = 0 ...(6) 

And F" + pF = 0 ...(7) 

Using (2), (4) gives X(x) F(0) = 0 and X(x) Y(b) = 0 

So that F(0) = 0 and F(6) = 0 ...(8) 

Where we have taken X(x) ^ 0, since otherwise u- 0 does not satisfy (36) 

We now solve (7) under boundary conditions(8). Three cases arise. 

Case I: Let p = 0. Then, solution of (7) is F(y) = Ay + B ...(9) 

Using B.C (8), (9) gives 0 = 5 and Ab + B = 0 
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This give A = B = 0 so that Y (y) = 0 

Thisleads to u = 0 which does not satisfy 3(6). So we rejectp = 0 
Case II: Let p = -X (negative). Here X * 0. Then solution of (7) is 
Y(y) = Ae Xy +Be~ Xy ...(10) 

Using 5.C(8),(10) gives 0 = zH- £ and 0 = Ae u + Be~ U ...(l l) 

(11) => ^ = 5 = 0 so that X (x) = 0. This leads to u = 0 which does not 
satisfy 3(6). so we reject p = -L 2 

Case III: Let p = A, 2 (positive). Here A.^0. Then solution of (7) is 

V \ V 

F(y) = i4cosX ! y+5sinA ! y '""^.(12) 

Using B.C (8), (12) gives 0 = A and 0 = £ sin 7.6 so that sinT.6 - 0, where 
we have taken £ * 0, since otherwise we shall get r(y)-=0. " - 

This leads to u = 0 which does not satisfy 3(6).. \. 

\ \ \ 

Now, sin Xb = 0 => Xb = nn \ \ \ 

V 

^X = ^,n = 1,2,... 

\ 7 ' — 

Hence, non zero solution Y n [y) of (7) are'given by 
T«= sinful; n = 1,2,3... \ ^V. -(13) 

2 2 

Also then p = X 2 = — . Hence (6) reduces to 
6 

r 2 2) 

X"- x = 0:,n = l2... ...(14) 

b 2 


whose solutions are X n (x) = C n e b +D n e b ...(15 

Using 3(6), (4) gives X (a) 7(y) = 0 so that X(a) = 0 where we have 
taken T(y) / 0, since otherwise u = 0 which does not satisfy 3(a). 

Now, X[a) = 0 => X n (a) = 0. Putting x = a in (15) and using X n (a) = 0. 
,we get 


C„e b +D n e b =0 




-D„ e b 


; 


Then(l5) reduces to 


■/ 

\^r..v :>• . ...: : ■; -' .2 

vL'. - - 
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we see that a solution u } 


are new arbitrary constants. In order to obtain 


Where E. 


solution also satisfying B.C 3( 0 ), we consider more geherarsolution 


Which is the half range fourier sine series of '/(^j , in (0,h). Hence 


Hence (18) is the required solution where in E n is given 
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Green’s Function 


CHAPTER-7 

GREEN’S FUNCTION 

7.1. Construction of Green’s Function 

Consider a differential equation of order n : 

L(u)=P 0 (. x)u n + P { (x)u n ~'+P 2 ( x )u n ~ 2 +... + P n ( x )u = 0, ...(1) 

Where the function /^(x), /}(,x), P 2 (x),., ./^(x) are, continuous on 

[a, b], P 0 [x)*0 on [a, 6] and the boundary conditions \ ^ 

y , X' v 'N Vw 

V k (u) = a k u(a) +a k l u'(a)+a k 2 u"(a)+ . +a k n ~ l u n ~ l (a) +b k u(b)+b k v!{b)+b k u{b) +?X,„ 

+b k n ~ x u n ~ x (b), (k = \,2,...,n) " -(2) 

V".\\ Y ' 

Where the linear forms V h V 2 ,...V n in 

_i \ " 

u(a ), u'(a),...,u n (a), n(b), u'(b),...,u n ~ l (b)'are linearly independent. 

N>- 

The homogeneous boundary value problem, (f);-{2) contains only a trivial 
solution u(x) = 0. . \v._ 

\ 

Green’s function of the boundary value prpblem (1), (2) is the function 
G(x, £,) constructed in x for any point £>a,< x < b satisfying the following 

u(x) = 0. 

(i) . G(x, 4) is continuous in x for fixed £, and has continuous derivatives 

with regard to x upto order (w-2) inclusive for a<x<b. 

(ii) . Its (n-l)th derivative with regard to x at the point x = £, has a 

discontinuity of first kind, the jump being equal to - f -^—, 




TT g ( x ’$) 


[a>(*)] x=£ ’ 


>o(0’ 


Where G| x= . +0 defines the limit G(x, as x -> £ from the right and 
G\ x= ^_ 0 defines the limit of G(x, t) as x from the left. 

(iii) . In each of the intervals [a, and (Y A] the function G(x, , 

considered as a function of x, is a solution of the equation (1): 

L(G) = 0...(4) 

(iv) . The function G(x, £) satisfies the boundary conditions (2) 

V k (G) = 0,k = \, 2, 3,..., n ...(5) 

If the boundary value problem (1), (2) contains only the trivial solution 
m(x) = 0 then the operator L contains one and only one Green’s function 
G(x, i;). 
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"Consider" i/j (3c), wf (jc) ,...,u n (x) be linearly independent solutions of the 
equation L(u) = O'. from the condition m. The known Green’s function 
G(x, t) must have the representation on the intervals [a , and (^, b ]: 

G (x, £) = a x u x (x) + a 2 u 2 (x) +... + a n u n (x) a < x < \ 

And (7(x, ty = l\u x (x)+b 2 u 2 { < x)+...+b n u n {x),£ 1 <x<b, Where 

a u a 2 ,...,a n , b x ,b 2 ,...,b n are some functions of £. From the condition I. The 
continuity of the function G(x, £) and of its first (n- 2) derivatives with 
regard to x at the point x = % yields. 

[Vi (^) + b 2 u 2 (%) + ...+b n u n (^)]-[«!«! ($) + a 2 u 2 (4) +... + a n u„ (4)] = 0 

.. x 

[ Vi' (5) + b 2 u' 2 {^) + ... + b n u' n (4)] ~[a x u[ (4) + a 2 u 2 (§) +... + ='0. 

N. \ 

[ Vf(4) + b 2 u 2 (5) + • • • + KK (^)]-[ a l«r($) + a 2 u 2 =0 



[^«r 2 (^)+*2«r 2 (4)+. --+^«r 2 (4)]-[a lM r^(^) (^)]=o 

Also \ ^ 

[%r‘ {%)+hu n i x (4)+...+6„«r l (y]-[«i«r‘ (y] ■=-^ 

Assume Q(£)=fy.(£)-ajfc(£)>£=l>2,-‘\, then, the system of linear equations 
in Q(^) are obtained. 

c 1 « 1 (^) + c 2 « 2 (^) + ... + c a (^)=o 

c,«i'U)+c , 2 «H§)+...+cx(4)=o 


c,«r ! (5)+c,„r ! (y+... ♦ (41=0 

A«r’ u)+c 2 »r 1 ©+...+c,»r' = --^r . w 

mas] 

The determinant of the system is equal to the value of the Wronskian 
(T(u,,« 2 ,. ,.,u n ) at the point x = £, and is therefore different from zero. 

From the boundary conditions (2), we have 

F t («) = 4 («) + **(«) ...(7) 

Where A k («) = a k u (a) + a x k u (a) + a 2 k u"(a) +... + a k n ~ l u n ~ i (a) 

Bk{ u ) = Bk u (b) + Blu'(b) + B%u"(b)+... + B% l u n 
Using the condition TV, we have 

V k{G) = a x Ak (u l )+a 2 A k (u 2 )+...+a„A k (u n )+.. .+^ (u [ )+b 2 B k (u 2 )+.. .+b n B k (u n )~ 0 
Where k = 1 , 2 ,... ,n 
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Since a k =b k -c k , then we have - |j| 

( 6 ,- c ,) 4 ( u ,)+( 6 2 -c 2 )A k (u 2 ) + ...+(b„-c„)A k (u„) + b i B t (u l )+b 2 B k (u 2 ) + ...+b„B k (u„) = 0 

' . . >■' 
=> h V k{ u \) +b 2 V k( u l) + --- + b n V k{ U n) 

= c l A k (u l )+c 2 A k (u 2 )+...+c n A k (u n ), ...( 8 ) 

Which is a linear system in the quantities b l ,b 2 ,...,b n . The determinant of the 0^ 
system is different from zero: 




^ ’ r‘-’ V 

v? £-i?r*■‘7--' L . ~ 


w 

VM) 

K(u 2 ) - 

V 2 (u 2 ) . 

.... V { {u n ) 

. r 2 («„) 

*0 


KM 

KM •• 

. KM) 


\ 

X V- " 




The system of the equations ( 8 ) contain a. unique'•-solution., in 
MSMi (*)»••• A (S) and since , Xf 

a k (£) = b k(ty~ c k{K)’ ^ f°ll° ws that the quantitie'SC;%(q) are defined 
uniquely. \ ■ 1 

(i) . If the boundary value problem ( 1 ), (2) is self- adjoint, then Green’s 

function is symmetric i.e., 

G(x, t,) = G&x) \\\ 

V V*— 

The converse is true as well. \ • 

,_ ij-.:-. 

(ii) .If at one of the extremities of ah^ interval [a, b\ the coefficient of the 
derivative vanishes e.g., P 0 (a) = 0, then.the natural boundary condition 
for the boundedness of the solution x = a is imposed, and at the other 
extremity the ordinary condition is specified. 

• Particular Case: We shall construct the Green’s function G(x, 
for a given number E,, for the second differential equation 

Z,(«) + <K x ) = 0 • •••(!) X:'; 

/ \ X 

Where L = - 7 - P~\+q ...( 2 ) & 

dx \ ax) 

Together with the homogeneous boundary conditions of the form 

au + p— = 0 .. .(3) 

dx K 1 



The Green’s function G(x, £,) constructed for any point a <x<b contains 
the following properties. 

(i) . G, (^) = G 2 (4) J h follows that the function G(x, is continuous in x 

for fixed \, in particular, continuous at the point x = t,. 

(ii) . The derivatives of G (which are of finite magnitude) are continuous at 

every point within the range of except at x = ^ where it is continuous 
so that 

g;r)-g;k)=--4 
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(iii) . The function G l and G 2 satisfy the homogeneous conditions at the end 

point x = a and x = b respectively. 

(iv) . The functions G x and G 2 satisfy the homogeneous equations LG = 0 l&v' ; v./'.' '.'• 

in their defined intervals except at z = £, i.e., 



LG l =0,x<£, 
LG 7 = 0, x>Z, 




Consider the Green’s function G(x, t) exists, then the solution of the given 7^777'/ }J r y:.yyy 

•"* "* 'V* 

differential equation can be transformed to the relation 




u{x) = \G{x,i,)m d ^ 


. 


-(4) a 




Consider two linearly independent solutions of the'homogeneous equation 
L(u) = 0. Let w = W| (x) and u = u 2 (x) be the nori-trivial 'solutions of the - 

equation, which satisfy the homogeneous conditions at >= a and x-b 
respectively. V V\ \ §m% : - 

Consider the Green’s functions for the problem Fromthe conditions III and 
IV, in the form v \ s \ 


IV, in the form v \ . \ -^1.7-' 7 

G(x^)=\ ClUl ^’ X< ^’ \"'"V- ...(5) 

G[x, v \c 2 u 2 (x),x>^ C u m 

Where the constant C, and C 2 are chosen in :a-manner that the conditions I vi¬ 
and II are fulfilled. \ 

Thus, we have Sife i; - 

-WK r '- S'' - : 

Q« 2 (4)-q«,(^=o a## 


1 L'. * - 

^ 2 m 2 (£) - 00 - ~ ■ ■ *00 

The determinant of the system (6) is the Wronskian W\u x (\),u 2 (t)\ : . 

evaluated at the point x = % for linearly independent solutions k,(x) and MM p 
u 2 (x) and hence it is different from zero W(fy * 0. 


- «i(0 a 2 (E) 

'[MtWWhii .L -< 7 ) ftsPfe- 

“llv “ 21 s) . 

* The conditions that «, (x) and u 2 (x) satisfy the equation L(u) = 0 are ; , v 7; 

• •• 

(pu{)'+qu x = 0, ...(i) ^ 


: The conditions that u, (x) and u 2 (x) satisfy the equation L(u) = 0 are 


(pul)' + qu x =0, ...(i) 

and (/?M 2 y+^ 2=0 ...(ii) Si.^^ 

Multiplying (i) by « 2 and (ii) by and subtracting, we have 

“2 OS)' - «i (M / = 0 => [/>(«1«2 - «2«1)]’ = 0 => «1«2 - «2«l' = ( C //>) • 


•••(ii) wm 

■ - v '- 7' 




5-» U 1 (^) u 2 (^) -U 2 ( 4 ) u l ( 4 ) : 
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From th'e system (6), we have 
Cv=~u 2 (%),C 2 ^~u^) 

Thus the relation (5) reduces to 

G ( x , t ,) = \ , -(9) 

*>5 

This result breaks down iff C vanishes, so that u { and u 2 are linearly 
dependent, and hence are each multiples of a certain non-trivial function 
U{x). In this case, the function U(x) satisfies 
together with the end conditions at x = a, x = b. 

Converse: The integral equation 


u(x) = \G(x,Z,)^)dt, 

a x 

X 

Where G(x, £) are defined by the 
equation 

Z,(w) + <j)(jr) = 0 

Together with the prescribed boundary 
We know that 

i r* b 

u ( x ) = ~r I“i(^)“2( x ) < t , U)^ + J«i(- 


u'(x) = ]«2 ( x)u x (x)u 2 ($)<t>(^)d£, 


u"(x) = ~] « 2 (*)“i U)<K^+j u’{x)u 2 {K)t}{K)dk 


= -^[“2 (*)“l ( x )~ u l (*)«2 (*)]<!>(*) 

Since L(u) - p(x)u n {x) + /?(jt)t/(jt) + g(jt)w(jt). 


£»(*)- “ f { t «2 (*)}“! (5)♦(6)<<5+ J { i «2 (*))“2 ft )4(5)45 


Again, ^(x) and u 2 (x) satisfy L(u) = 0, hence the first two terms vanish 
identically. 

So Lu[x) = -<t>(x) => Z,k(x)+<|>(x) = 0 

Therefore, a function u(x) satisfying (10) also satisfies the differential 
equation (11). 
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Again from (12) and (13), we have 

a 

u X a ) = --^\ u 2fe)*fe) d $’ 

a 

Which shows that the function « defined by (11) satisfies the same 
homogeneous condition at x = a as the function 

Example: Find the Green’s function for the boundary value problem 



1^10 • 




a u ( \ n 

with the conditions u (0) = u (l) = 0 
Solution: The differential equation is given as 

^-^r-u(x)= 0 

dx 2 { ’ ' 

with u(0) = k(1) = 0 


V V 


"’S*. \v 


The general solution of the differential equdtiohisgiVen as 


u(x) = Ae x + Be 
«(0) = 0 =>A- 


A + B = 0 


«(l) = 0 =>Ae + Be l = 0 




A = B=0 


These conditions have only a trivial solution m(x) = 0, hence we can jt. 
construct a unique Green’s function. 

The fundamental system of solutions for the differential equation (1) is given 


u { (x) = sin hx Satisfies the boundary condition u (0) = 0 

u 2 (x) = sin h(x- 1) Satisfies the boundary condition u (l) = 0 

The solution are linearly independent. The value of Wronskian for sin hx 
and sin/i(x-l) at the point x = £; 


nth 


Isinh^ sinh(£ > — l) I 


v J |cosh^ cosh(^-l)| 

=> Pr(^) = sinh^cosh(^-l)-cosh^sinh(^-l) = sinhl 

Here p(x) = 1. From (9), we obtain the Green Function of the form: 






• -v 


K 






sin hx smh{^-\) 
sinM 

sinsin/*(x-l) 


, 0<x<^ 


, ^<x<\ 


- ■ 
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Example: If V(x) has continuous first and second derivatives and satisfies 

d 2 u 

the boundary value problem —-+= 0 with u(0)-u(/) = 0 then u(x) is 

dx 1 

continuous and satisfies the homogeneous linear integral equation 

i 

o 

Solution: The differential equation may be written as 


-+ku = 0 => -—— = —Xu 


By integrating with respect to x over the interval (0, x) two times, we obtain 

- = -x[ufe)d$+C XX'X\ 

dx i v „ \ 

i 

or m(x) = -X J(x-^)m(^)c/^+Cx+Z) \ vv "’'"h ■••(2) 

0 V —VX 

v \ \\ s 

Where C and D are the integration constants,- td .be determined by the 
boundary Conditions. \ \\ 

«(0) = 0 => D = 0 X. 

. \ 

«(/)=o^-x} o (/-^«(^)^+a=o r- 


C = yJ 


v X. 

\>N 


Substituting the value of the constants C and D in (2), we have 


«(x) = -X. j(x-^)M(^)^ + y|x(/-^)«(^ + y|x( / -^)«(^)^ 

0 ^0 x 


“W = ^ Jy( /_ JC )“(4)^4+Xjy( 1- ^)“(^) rf ^ 

O x 

l 

u(x) = ljG(x,Z,)u(£,)dt, 


where G(x, ^) = 


Example: Construct the Green’s function for the boundary value problem 
—- = 0 with the conditions w(0) = w'(l) and w'(0) = «(l). 


Solution: The differential equation is given by 
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$ 


...(4) 


with the conditions i , • r 

“(0) = «'(!) 

«'(°) = u(l) 

The general solution of the given differential equation is 
u(xj = Ax+B 

J The boundary value problem contains only a trivial solution w(x) * 0, hence 
we can construct a unique Green’s function for it. 

Consider the unknown Green’s function G(x, £) must have the 
representation on the interval [0, and (£, l] 

a x x + a 2 ,0<x<^ 

b l x + b 2 ,^<x<l 



G{x,£,) = 




>-( 5 ) 


Where a u a 2 , by, b 2 are unknown functions of £7 The Green’s function 
G(x, £) must satisfy the following properties. , W \ " 


I. G(x, £) is continuous at x = i.e., \ 


\ , 


a & + a 2 ~ b k + b 2 => (^i ~ a i)\ + {b 2 \=> C£ + C 2 -0 ...(6) 

II. The derivatives of G has a discontinuity \of magnitude —7-- at the 

‘ 1 po(v 


point x = E,, where p 0 (%) = co^ficient of highest order derivative in 
the given differential equation (= l)\v._ 


dG 

bx Jx=$+o 


i .-i 


v dx J x= ^ 0 


bi -ai =-l => Q = -1 and C 2 =£, 


...(7) 


III. The Green’s function G(x, £) must satisfy the prescribed boundary 
conditions 


G(<U) = G'(l4 ^ a 2 =^ 

G'(0,q) = G(U) => a l =b l +b 2 

Consider C k = b k -a^=> C x = 6, -a, and C 2 =b 2 -a 2 ...(10) 
From the equations (7) to (10), we have 

b 2~ a 2 =%> h- a \ = _1 

bi+b 2 =a l , a 2 =bi 

By solving these equations, we have 
ci\ =-^ + 2, a 2 =— fj + l, t\ i4 + lj b 2 = + 1 

The required Green’s function G(x, is obtained: 

fH + 2)* + B + l) 0<x<^ 


...( 8 ) 

...(9) 


G{x,$ = < 


,R + 1)* + 1 


£<x<l 












* IV*-' : 






W 4 :" 

■ 

■■ : 






*. ■ : ?. 


. 

■irA'-.l:- ‘.' 7 . •• • - 

h! ... . ' : . ' 

/ ■ . 

— . 




'-J ‘ . ' Z8A/ll, (Fir 

st Floor) Jia Sarai, Hauz.Khasi Near I.I.T.j New I>elhi-110016, PIl: (0U)-265J7527, Cell: 9999183434 & 9899161734,8588844789 

E-mall: infold iDsacademv.com;Website:ww.dipsacademvxom '* £**■:*• ■* 









ASSIGNMENT (SHEET-01) 

1. Consider the initial value problem 5. The Charpit’s equation for the PDE 


—+ 2—= 0, u(0,y) = 4e~ 2y . 
dx dy X ’ 

Then the value of «(l,l) is 

(a.) 4e' 2 
(b.) 4e 2 
(c.) 2e* 

(d.) 4e 4 

2. Let a,be K be such that a 2 +b 2 * 0. Then 
the Cauchy problem a^-+b^- = \\ x,yeR 

u[x,y} = x on K, \ ( 

\ ~~' v O 

ax + by = 1 \ 

% \ v 

(a.) has more than one solution if either a or''-- x \’ 
b is zero \ 


(b.) has no solution 

(c.) has a unique solution 

(d.) has infinitely many solutions 




Which of the following are complete 
integrals of the partial differential equation 

pqx + yq 2 =1? 


(a.) z = ± + ^ + fc 
a x 


(b.)z=JA* 

b x 


(c.) z 2 =4(ax + y) + b 

(d.) (z-b) 2 = 4(ajc + y) 

For an arbitrary continuously differentiable 
function /, which of the following is a 

general solution of z(px-qy) -y 2 -x 2 
(a.) x 2 +y 2 +z 2 =f(xy) 

(b.) (x + y) 2 + z 2 -f{xy) 

(c.) x 2 +y 2 +z 2 =f(y-x) 

(d.) x 2 +y 2 +z 2 =f((x+y) 2 +z 2 ) 


2 2 /V du du 1 

up +q +x+y = Q,p = —,4 =— are given by 
.dx dy 

(a) dx _ dy _ du _ dp _dq 

-l-/? 3 -1 -qp 2 2p 2 u + 2q 2 2 pu 2 q 

(b) dx _ dy _ du _ dp _ dq 

2pu 2 q 2p 2 u + 2q 2 -1 - p 3 -\~qp 2 

x , \ dx _dy _du _ dp _ dq 

. "o~ 


( d ) dx_ = dy _-. du dp dq 

\ 2 q 2pu x + y p 2 qp 2 ' 

Consider the Cauchy problem of finding 

, , , , du du . . 

u = u(x,t) such that — +u — = 0 for 

dt dx 

x £ > 0 u(x, 0) = k 0 (x), xeR which 

choice(s) of the following function for u 0 
yield a C 1 solution k(x,0 for all xeR and 


(a.) «„(x) = — T 
1+x 


(b.) M 0 (X) = X 

(c.) u 0 (x) = l + x 2 
(d.) u a (x) = 1 + 2x 

7. Let u = u(x, t) be the solution of the Cauchy 


du {du\ . . 

ui ~T + =1 xeR > ?>0 

problem dt \ dx) 

u{x, 0) = -x 2 xeR 


(a.) u{x,t) exists for all xeRand t> 0. 

(b.) |u(x,t)I—> co as t-yt* for some t*> 0 
and x*0 

(c.) u(x,t) <0 for all xeR and for all 
t< 1/4 

(d.) «(x,t)>0 for some xeRand 0<t<l/4 
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Assignment Sheet-1 


Let jc = x(s),y = y(s),« = u(s),s eR,be 12. The solution of the Cauchy problem for the 

dz dz 

the characteristic curve of the PDE first order PDE x —+v— = z, on 

/■ \ f \ ' dx dy 

f du)f . . 

Jl — I - w = ° passing through a given D = {(x,y,z)\x 2 +y 2 * 0 ,z> 0 }, with the 

n 2 D •. initial condition x 2 +y 2 = 1,2 = 1 is 

curve x = 0,y = r,« = r ,reR. Then the 

characteristics are given by (a.) z = x 2 + y 2 

(a.) x = 3r(e s -l),7 = |(e- I +l),« = rV 2s (b .)z = {x 2 +y 2 f 

(b.) T = 2r(e' s -l),y = r(2e 2j -l),«=y(l + e _2s ) (c.) z =(l-(x 2 + y 2 )f 2 

(c.) x-2z{e s = — [e s +l),w = rV* t4)..; z = ( x2 + /) 

... , 13>v Consider thefirst order PDE p + q = pq where 

(4) «-,(«--1)., — 2 r[«--r>(2.--1) 

^ ' N s> '"'P s —, q = —. Then which of the following 

The initial value . problem 'V\ \\ dx ^ 

du du ]" v \ \ are correct? 

— + x — = jc, 0 < jc < 1,/ > 0 and u (jc,0) = 2x .. 

dt dx 

(^ TVip PVmmit’e pniiatinne fnr flip ahnvf^ 


(b.) z = (x 2 +y 2 f 
(c.) z=(2-(x 2 + y 2 )) V 
\4j;;z = ( x 2 + y 2 j n 


(a.) A unique solution u(x,t) which 

~^CO OS t ~>CO \ 

V 

(b.) More than one solution \>, 

(c.)A solution which remains bounded as 
t —>co 

(d.) No solution 

10. Let xyu = c, and x 2 + y 1 - 2u = c 2 , where 

c, and c 2 are arbitrary constants, be the first 
integrals of the PDE 

x ( u+ y 1 )j x -y( u+xl )^=( x ‘-y t ) u ' 

Then the solution of the PDE with 
x + y = 0, w = 1 is given by 

(a.) Jt 3 + y 3 + 2xyu 2 + lx 2 u = 0 

(b.) x 3 + yx l + [x 2 + xy j u = 0 

(c.) x 2 + y 2 +2(xy-l)n + 2 = 0 

(d.) x 2 -y 2 -u(x +y-2)-2 = 0 

11. Let z = z(x,y) be a solution of —— = 1 

dxdy 

passing through (0,0,0). Then z(0,l) is 

, (a-)0. 

(b ) L 
(c.) 2. 

(d.)4. 


^ d = ——, q = —. Then which of the following 
dx dy 

are correct? 

(a.) The Charpit’s equations for the above 
PDE reduce to 

dx _ dy _ dz _ dp _dq 

i -q i -p -pq p+q o 

(b.)A solution of the Charpit’s equation is 
q=b, where b is a constant. 

(c.) The corresponding value of p is 


(d.) A solution of the equation is 

z = -^— x+by + a, where a and b are 
b -1 7 

constants. 

14. The Cauchy problem 
xu x +yu y =0 

u(x, j’) = x, on x 2 + y 2 =1 has 

(a.) A solution for all x e R, y e R 

(b.)An unique solution . in 

{(*,y) e tf 2 :(x,y)*(0,0)} 

(c.) A bounded solution in 

{(x,y)e.R 2 :(;t,y)*(0,0)} 

(d.) An unique solution in 


{{x,y)eR 2 i(x,y)*m} 


but the solution is unbounded. 
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15. Let u = u(x,y) be the complete integral of the 

PDE — • — = x y passing through the points 
dx dy 


(0, 0 , l) and |T), 1 , -J in the x-y-u space. 
Then the value of w(x, y) evaluated at (-1, l) 
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19. The solution of the Cauchy problem 
u x + u y = 0, u(x, 0) - e x is 

(a.) u(x,y) = e x+y 
(b.) u(x, y) = e x e 2y 
(C .)u{x,y) = e‘-> 

(d.) u{x,y) = ^[e x * y +e x - y ] 

20. The solution of the PDE*—+ y— = 3u with 

, dx dy 


(c-)2 
(d.) 3 

16. The differential equation 

du du n. '■ 

x—+y — = 2 u \ 

dx dy _ 

satisfying the initial condition " 

y = xg[x),u= f[x) with \ \ 

(a.) /(x) = 2x, g(x) = 1, has no solution ’\ x N „' w 

V ‘ X V"'— 

(b.) /(x) = 2x 2 ,g(x) = l, has infinite number 
of solutions v‘\ V*" 


\ V 

(c.) /(x) = x\g(x) = x, has a unique solution 


initial data u (x, 1) = * (1 - x) is 
. '\(a>) % xy 2 -xfy 

xy 1 -xfyL 

o-\T*) xg-xv \ 

' y 2K>^ The complete integral of the non-linear first 
\ \ '• order partial differential equation 


y 


(l + x 2 )p 2 -qx 2 = 0 is 


(d.) /(x) = x 4 ,g(x) = x, has a unique 
solution. 

17. A general solution of the PDE 
uu I + yu = * is of the form 


(a.) / u 2 -x 2 ,^— =0, where /:r 2 -»R 
V x+u J 

is c‘and V/ * (0, 0) at every point 


(a.) z = ta-v/l + x 2 + - — +c 
2 

a 2 

(b.) z = ± —= +a y+c 
VI+ x 2 

(c.) z = ±-a(l + x 2 f' 2 +^-+c 
3 2 

(d.) z=±a\fl + x 2 +a 2 y 2 +c 


(b.) u 2 = g - +x 2 , g e 0 (R) 

\x+u J 

(c.) /(u 2 +x 2 ) = 0,/eC 1 (R) 

(d.) f(x+y)= 0, /eC'(R) 

18. The partial differential equation 

(“J 4 ^+“ w +(“,) 2 =0 is 

(a.) Linear and of order 4 
(b.) Quasi-linear and of order 2 
(c.) Quasi-linear and of order 4 
(d.) Linear and of order 2. 


22. The characteristics of the partial differential 

. du du 

equation —+ — = u are 

dx dy 

(a.) Exponential curves along which u 
remains constant 

(b.) Exponential curves along which u varies 
linearly with x 

(c.) Straight lines along which u remains 
constant 

(d.) Straight lines along which u varies 
exponentially with x 

23. For obtaining the complete integral of the 
non-linear partial differential equation 

pq(x 4 +y 4 ) = (p + q) , the required forms of 
dz dz • i • 

p-—,q = — , involving an arbitrary 

dx dy 

constant a, are 
(a-) p = (x 4 + a)~ { , q = (y 3 - a)~ l 
(b.) p = (2x A +ay\q = (3y i -a)- 1 
(c ■) p = (-2x 4 +ay\q = (-3/ - a)~ x 
(d.) p = (x 4 - a)" 1 , q = (> 3 - a)' 1 
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24. The general integral of 

x(y-z)^-+y(z-x)^- = z(x-y) is 
ox dy 

(a.) F(2x+3y + z,xyz) = 0 

(b.) F(x + y+z,xyz) = Q 

(c.) F(x-y + z,xyz) = 0 

(d.) F(x+y-z,xyz) = 0 

25. The general solution of 

.dz '& , _ ... 

2— + 3— = cos(3jc-2>’) is 
cbc dy 

(a.) 2z = x cos(3x - 2y) + /( 3x + 2y) 

(b.) 2z = xsin(3x-2y) +f(3x + 2y) 

(c.) 2z = xcos(?>x-2y) + f(3x-2y) 

(d.) 2z = x sin(3x - 2 y) + /(3.x - 2y) V ' '‘7^-^ 

"'"v. 

26. The general integral of z(xp-yq) = y 2 - x 2 is \ 
(a.) z 1 =x 2 +y 2 + f(xy) 

(b.) z 2 — x 2 -y 2 +/(xy) 

(c .) z 2 =-x 2 -y 2 +f(xy) ^ \ 

(d.) z 2 =y 2 -z 2 +f(xy) 

27. A singular solution of the partial differential . . 
equation z + xp - x 2 yq 2 - x 2 pq = 0 is 


(a.)z = - 

y 


(c.)z=4 


(d.) z = — 

X 

28. Let w(x 5> y) = 2/(7)cos(x-2^),(jc,^)6M 2 , be 
a solution of the initial value problem 

2u x +u y =u u{x, 0) = cos(x). Then /(l) is 
equal to 


Assignment Sheet • 


29. The integral surface of the first order partial 
differential equation 

2y(z - 3)—+ (2x - z) — = y(2x - 3) passing 
dy dy 

through the curve x 2 + y 2 = 2 x, z =■ 0 is 
(a.) x 2 + y 2 - z 2 - 2x + 4z = 0 
(b.) x 2 +y 2 -z 2 -2x + 8z = 0 
(c.) x 2 +y 2 -z 2 -2x + 16z = 0 

\ (d.) x 2 + y 2 + z 2 — 2x +8z — 0 

30. The integral surface satisfying the partial 

\>v v differential-..equation —+ z 2 — = 0 and 
. ^ dx dy 

■ v - passing through the straight line x = \,y = z 


( a -) (x-l)z + z 2 = y 2 
„ _ “ (b.) x 2 +y 2 — z 2 =1 

(c.) (y-z)x + x l = 1 
(d.) (jc - l)z 2 + z = y 

31. The integral surface for the Cauchy problem 

dz dz 
—+— = 1 , 
dx dy 

which passes through the circle 
z = 0, x 2 + y 2 = 1 is 

(a.) a' 2 + y 2 +2z 2 + 2zx-,2yz-l = 0 

(b.) x 2 +y 2 + 2z 2 + 2zjc + 2yz -1 = 0 

(c.) x 2 + y 2 + 2z 2 - 2 zx - 2 yz -1 = 0 

(d.) x 2 + y 2 + 2z 2 + 2zjc + 2 yz + 1 = 0 

32. The integral surface satisfying the equation 

y—+x— = x 2 +y 2 and passing through the 
dx dy 

curve x = l-r,y = l + /,z = l+t 2 is 


(d.)- 

2 


(a.) z = xy + 




(b.) z = xy + -(x 2 -y 2 f 
(c.) z = xy + ^(x 2 -y 2 f 

O 

(d.) z = xy + ^-(x 2 -y 2 f 
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' The solution of xu x + yu - 0 is of the form 
(a.) f(ylx) 

(b.) /(*+>') 

(c.) f(x-y) 

(d-) f(xy) 

The characteristic curve of 
2 yu x + (lx+y 2 )u y =0 passing through (0,0) is 

(a.) y 2 = 2(e* +x-l) 

(b-) y 2 = 2(e* -x + 1) 

(c.) /= 2(e x -x-l) 

(d.) y 2 = 2(e x +x + l) 

The initial value problem 
u x +u y =1, «(j,j) = sini',0<s<l has 

(a.) Two solutions 
(b.) A unique solution 
(c.) No solution 

(d.) Infinitely many solutions • 

Let «(x,y) = /(^) + g(y 2 cos(y)) \ 

s 

Where / and g are infinitely differentiable 

\ 

functions. Then the partial differential \ 

equation of minimum order satisfied by, u is - 

\ ^ 

(a.) u xy + xu xx =u x ^ 

(b.) u v + xu a =xu s 
(c ■) U v -XU a =U x 

(d.) u xy -xu a =xu x 

Let u(x, y) be the solution of the Cauchy 
problem 

xu x +u y = 1, u(x, 0) = 2 ln(x), x > 1 
Then u (e, l) = 

(a.) -1 
(b.)0 
(c.) 1 
(d.)e 

Consider the partial differential equation 

du du . . _ . 

— + u — = 0 satisfying the initial 

dt fix 

condition u(x,0) = a + fix . If m(x,/) = 1 

along the characteristic x = t +1, then 

(a.) a = 1,0 = 1 

(b.) a = 2,0 = O 

(c.) a = 0,0 = 0 

(d.)a = O,0 = l 


Consider the partial differential equation 

du du . . _ 

x - y— = u .the characteristic curves for 

dy dx 

the above equation in the (x, yj plane are 

(a.) Straight line with slopes 1 
(b.) Straight-lines with slopes -1 
(c.) Circles with centre at the origin 
(d.) Circles touching y axis and centered on 
x-axis 

If u(x,y) is a solution to the partial 

, . du du .. 

differential equation x - y — = u with 

\0 dy dx 


'u(x,0) = s ilt —x , 


— r then 


V2’V2 


equals 


" <1 ' ) lV 

(b.) -e* 

4 

(c -> t / 1 

<“•) r 2 

The integral surface of the partial differential 

du du . ... , 

equation x — + y — = 0 satisfying the 
dx dy 

condition u(\,y) = y is given by 
(a.) u(x,y) = ^ 

(b.) u(x,y) = -^- 
x + \ 

(c.) 

(d.) u(x,y) = y + x- 1 

The integral surface of the partial differential 

equation x— + y— = 0 satisfying the 

M dx dy 

condition w (l, y) = y is given by 

(a .)«(*,,) = ■£ 

(b .) u (x,y) = -&- 

Jt + 1 

(c.) «(x,,) = -X- 


(d.) u(x,y) = y + x- 1 
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43. Complete integral for the partial differential 
equation z = px + qy - sin (pq) is 

(a.) z = ax + by + sin(ab) r 

(b.) z = ax + by- sin ( ab ) 

(c.) z = «* + >' +sin (6) 

(d.) z = .T + Z>v-sin(<i) 

44. The Cauchy problem u x -u y =2 with the 

Cauchy data on r: 2,?) has 

(a.) One solution 
(b.) Two solutions 
(c.) No solution 
(d.) Infinite solutions 

45. The general integral of the partial differential 
equation (y + zx)z x -(x + yz)z y 

is . - 

(a.) F(jc 2 +y 2 + z 2 ,xy + zj = G , 

(b.) F (jc 2 + y 2 - z 2 , xy + z) = 0 \ 

(c.) F(x 2 -y 2 -z 2 ,xy + z} = 0 ^ 

(d.) F (x 2 + y 2 + z 2 , xy - z) = 0 

Where F is an arbitrary function. 

46. The partial differential equation of the family 
of surfaces z = (x + y) + A(xy) is 

(a.) xp-yq = 0 
(b .) xp-yq = x-y 
(c.) xp+yq=x+y 
(d.) xp + yq = 0 

47. Using the transformation u = — in the partial 

.v 

differential equation xu x =u+yu y , the 

transformed equation has a solution of the 
form w = 


(a.) /(-) 
y 


Assignment Sheet -1 


48. A general solution of the equation 


(b.) f(x + y) 
(C-) f(x-y) 
(d-) f{xy) 


-u(x,y) = e x is 


(a.) u(x,y) = e x f{y) 

(b.) u(x, y) = e~ x f(y) + xe x 

(C.) u(x,y) = e x f(y) + xe- x 

(d.) u(x,y) = e~ x f(y) + xe~ x 

49. The integral curves of the equation 
'v . P+2q = 5z+cosh(2x-y) are given by the 
\-'-.^ ''Hntersectibnpf the surfaces 

V> N S ' '" ',(a.) 2x-jy = c 1 ;[Sz + cosh(2x-y)]e~ 5x =c 2 

i*v vV 

■ • (b.) x-2y = C\,[5z + cosh(x-2y)\e~ 5x = c 2 

H l (c.) 2x - y = q; [z + cosh (2x -y)~\e x = c 2 

(d.) x-2y = q;[z + cosh(2x--y)]e _ * - c 2 

S ‘ 50. The complete integral of the partial 
differential equation 

xp 3 q 2 + yp 2 q 3 + (p 3 + q 3 }-zp 2 q 2 = 0 is z = 

(a.) ax + by+ {ab~ 2 + ba~ 2) } 

(b.) ax-by + [ab~ 2 -ba ~ 2 ) 

(c.) -ax + by + (ba~ 2 - ab~ 2 ) 

(d.) ax+by-[ab~ 2 + ba~ 2 ^j 

51. For arbitrary real valued smooth functions / 
and g, the function u defined as 
u (x,t) = f(x+t)+g{x-t) is a general 
solution of 

(a.)u t +u x =Q 
(b.) u t -u xx = 0 

(C-)«ff+K*= 0 
(d.) u tt -u xx =0 
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ASSIGNMENT (SHEET-02) 


The variables £ and rj which reduce the 

differential equation ^--x 2 ^- = 0 to the 

dx 2 dy 2 

canonical form are 

(a.) Z = y 2 +^x,i = y 2 ~x 

(b.) Z = y+^x 2 ,7 1 = y-^x 2 

(c.) Z=y + x 2 ,7j=y-x 2 
(d.) Z = yyx 2 ,tj = y 2 -x 
If the partial differential equation 

(x —l) 2 -(y-2) 2 u yy +2xu x +2yu y +2xyu=0 

is parabolic in S c R 1 but not in R 2 \S, then 


6. Which of the following is elliptic? 

(a.) Laplace equation 
(b.) Wave equation 
(c.) Heat equation 
(d.) u xx +2u xy -4u yy =0 

7. The equation 

* 2 {y- !) z « -*(y 2 ~') z xy+y(y 2 ~ l ) z yy +z x = 

\ 

is hyperbolic in the entire xy-plane except 


(a.) j(x,y) g R 2 : x = 1 or y = 2| 
(b.) j(x,y) g R 2 : x = 1 and y = 2 


(c.) {(x,y)Gi? 2 :x = lJ 

(d.) [(x,y)sR 2 -.y = 2\ , \ 

i V 

The characteristics curves of the Equation 
x 2 u xx -y 2 u yy -x 2 y 2 +x; x>0,« = (x,y) arc 

(a.) Rectangular hyperbola 
(b.) Parabola 
(c.) Circle 
(d.) Straight line 

If f(x) and g(y) are arbitrary functions, 
then the general solution of the partial 

differential equation u ^-^-——— = 0 is 

dxdy dx dy 

given by 

(a.) a(x,y) = f \x) + g{y) 

(b.) u{x,y) = f{x + y) + g{x-y) 

(c.) u(x,y) = f(x)g(y) 

(d.) u(x,y) = xg{y) + yf(x) 

The number of characteristic curves of the 
PDE 

( x 2 + 2y)u xx + (y 3 -y + u) Uyy 

+x 2 {y - l)w^, +3u x +u = 0 
passing through the point x = 1, y = 1 is 
(a .)0 
(b.) 1 
(c.) 2 
(d.)3 


V\ (a:) x-axis.. 

V s \ (&'•) y - axis 

s ' 

x ->\(c.) A line parallel toy -axis 

. (d.) A line parallel to x -axis 

8. \ The second order partial differential equation 

**. 2 

X ( x-y ) d 2 u , \ ■ ( 2 2 ) d 2 u 

4 dx 2 V ' 1 ’dxdy 


+ cos 2 (x 2 +y 2 )^-^-+(x-y) 


5y 2 v 'dx 


region 


+sin 2 (x 2 + y 2 +u = 0 is 

(a.) Elliptic in tl 
|(x,y):x*y > x 2 +y 2 <|j 


(b.) Hyperbolic in the region 
{(x,y):x*y,^<x 2 +y 2 <^j 


(c.) Elliptic in the 
{(x,y):x*y,J<x 2 +y 2 < 


region 


(d.) Hyperbolic in the region 
j(x, y):x*.y,x 2 +y 2 <^| 

9. The general solution of the partial differential 

equation [D + D' -2)(Z>-D')z = 0 , where 

d d 

D = — and D = — is 

dx dy 


(a.) fi( x+ y) +elx fi( x -y) 
(b.) e2x f 2 (x+y)+f 2 {x-y) 

(c.) e ~ 2x M x+ y) + fi( x ~y) 
(d.) fi( x+ y) +e2x fi( x -y) 
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partial differential equation 


d 2 u d 2 u 


y —-+ x — - = 0 is hyperbolic in 

dx 2 dy 1 

(a.) The second and fourth quadrants 
(b.) The first and second quadrants 
(c.) The second and third quadrants 
(d.) The first and third quadrants 

The general solution of X- + —r = 0 is of 

dx 2 dy 2 

the form 

(a.) u = f(x+iy)+g(x-iy) 

(b.) u = f(x+y)+g(x-y) 

(c.) u = cf(x-iy) 

(d.) u = g(x+iy) 

The second order partial differential equation 

x ^x, 

+ 4 Ufy +(cosx)« xv + e x u x + e y u y = 0 has , ~ 

(a.) Exactly one characteristic curve passing ’■ 
through every point V 

■v. 

(b.) Exactly one characteristics curve passing s 
through the point (0, 0) \ 

(c.) Two distinct characteristic curvesmas§ing"v^- 
through every point \ X 

(d.) Three distinct characteristic cuhfes— 
passing through the point 
The partial differential equation 


Uxx +{x+y)u + 


x 2 + y 2 


Uyy + (sill X ) =0, 


(a.) Elliptic in {(x, y) :xe R ,y e R} 

(b.) Parabolic in {(x, y): x e R, y > 0} 

(c.) Hyperbolic in {(x, y):xy> 0} 

(d.) Parabolic in {(x, y): xy <0} 

The complete integral of the PDE 

d 2 u . d 2 u d 2 u x+v 
dx 2 dxdy dy 2 

involving arbitrary functions 4>! and <t> 2 is 
(a.) ^ (y+x)+<t> 2 (y + *) +^e x+y 

(b.) <)>[ (y+x)+x(j)2 (y+x) +i?L-}-l e x+ y 

(c.) <f>i (y-x) + <j> 2 (y-x)+^e x+y 
(d.) <t>i (y -x) + x «|)2 (y -x) + e x+ >’ 


I 


In the region x > 0 ,y >0 the partial 
differential equation 




2\ d 2 u 


V\ V 


v 'ctc z +2 v ’dxdy 

X 2 )0=» 

(a.) Changes type 
(b.) Is elliptic 
(c.) Is parabolic 
(d.) Is hyperbolic 

Pick the region in which the following 
differential equation is hyperbolic 

:-Xxx+' 2 ^. =» x +u y 

(a.) xy# 1 
'-(b.) xy#0 " 

"(c-) xy> 1 
(d.) xy > 0 


The expression 


dI-d] 


sin(x-y) is equal to 


(a.) cos(x-y) 

X 

(b.) sin(x - y) + cos(x - y) 

x 

(c.) cos(x-y) + sin(x-y) 

(d.) ysin(x-y) 

A general solution of the second order 
equation 4 u^-u^ = 0 is of the form (where 

/ and g are twice differentiable functions.) 

(a.) f(x)+g(y) 

(b.) f(x + 2y) + g(x-2y) 

(c.) /(x + 4y) + g(x-4y) 

(d.) / (4x + y) + g(4x - y) 

The general solution of the partial differential 

equation = x+y is of the form 

dxdy 

(a.) ^xy(x + y) + F(x) + G(y) 

(t>.) ^ xy{x-y)+F(x) + G{y) 

(c.) ixy(x-y) + F(x)G(y) 

(d.) ^xy(x+y) + F(x)G(y) 
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20. The characteristics of the partial differential- 

equation 36^4-/ 4 ^-7/ 3 ^ = 0 , are 

dx 2 dy dy 

given by 

(a.) x+-£=c l ,x-—=c 2 

y y 

/u \ 36 36 
(b.) x+— = c u x-— = c 2 

y y 

7 7 

(c.) 6x+— = c x ,6x—-r = c 2 

y y 

i i 

(d.) 6x + -r = c u 6x -r- = c 2 

y y 

21. The second order partial differential equation 

(a.) Hyperbolic in the second and the fourth 
quadrants 

(b.) Elliptic in the first and the third . 
quadrants \ 

(c.) Hyperbolic in the second and elliptic, in '' 
the fourth quadrant % Xw 

(d.) Hyperbolic in the first and the third" 
quadrants 

22. The partial differential equation 

d 2 u „ d 2 u d 2 u du du . ■ 

x—- + 2 xy -+ y —- + *— +v— = 0 is 

dx 2 dxdy dy 2 dy dx 

(a.) Elliptic in the region jc < 0,y < 0 ,xy > 1 

(b.) Elliptic in the region x > 0,y > 0 ,xy > 1 

(c.) Parabolic in the region x< 0,y <0,xy>l 


(d.) Hyperbolic in 
jt<0,y<0,xy>l 


region 


23. The partial differential equation 

d 2 u d 2 u . . 

—r-y—r = 0 has 

dy dx 2 

(a.) Two families or real characteristics curve 
for y < 0 

(b.) No real characteristics for y > 0 

(c.) Vertical lines as a family for 
characteristics curve for y = 0 

(d.) Branches of quadratic curves as 
characteristics for y* 0 


If u(x,t) satisfy the partial differential 

... d 2 u d 2 u , , s. , c 

equation — -=4 —- then u(x, t) can be of 
dt dx 2 

the form 

(a.) u(x,t) = f{e x ~ 2 ‘) + g(x+2t) 

(b.) u(x,t) = fix 2 -4t 2 )+g(x 2 +4t 2 ) 

(c.) u(x,t) = f(2x-4t)+g(x+2t) 

(d.) u(x,t) = f{2x-t) + g{2x+t) 

(Where / and g are non-trivial smooth 
function). 

The second order partial differential equation 

: ,' U^ +XM W = 0 is 

'x(a.) elliptic for y > 0 
- .(b.) hyperbolic for x > 0 
,. (c.) elliptic for x < 0 
(d.) hyperbolic for x < 0 
Let P(x, y) be a particular integral of the 

partial differential equation 

2 

^r~— = 2y-y 2 ,Then P{ 2; 3) equals 

dx 2 dy 

(a.) 2. 

(b.)8. 

(c.) 12. 

(d.) 10.- 

2 2 2 

T , n „ r . d u „ d u d u .. 

ThePDEis —+2-+—^ = 0is 

dx 2 dxdy dy 2 

(a.) Parabolic and has characteristics 

l,(x,y) = x+2y,x\{x,y) = x-2y 
(b.) Reducible to the canonical form 

d 2 

—^- = 0, where £,(x,y) = x+2y 
d% 2 

(c.) Reducible to the canonical form 

—- = 0, where r\(x,y) = x+2y 
dr\ 

(d.) Parabolic and has the general solution 
u = (x-y)f{x+y)+f 2 {x-y) where 
f,f 2 are arbitrary functions. 

Canonical form of r-t- 0 is/are 

(a.)^ = 0 
dV 2 

(b.)_^!_=o 

dUdV 

(c.) — = 0 ' 

dU 

(d.)—= 0 
dV 
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29. Partial 


Differential 


Equation 


Uxx + Uyy +U 22 + Uyz + U^ = 0 


(a.) Hyperbolic , 

(b.) Parabolic 

(c.) Elliptic 

(d.) Laplace equation 

30. The characteristics of the equation 

J&A 


-y —4 =0 are 

W) 


(a.) x+y ,x-y 


(b.) - £ 

y y 


(C )*y,~ 
y 


(d.) xy , x 2 y 2 

31. The second order partial differential, 
equation. 


. I —. d 2 U d 2 U . I — d 2 U 

(i-'Jxy)TT+ 2 ^T-+( [+ yl x y)—r = ois 

ox oxoy ay 


\ 


(a.) Hyperbolic in the second and the fourth N ‘ 

quadrants. ' \__ 

(b.) Elliptic in the first and the\ tliird 
quadrants. 

(c.) Hyperbolic in the second and elliptic in 
the fourth quadrant. 

(d.) Hyperbolic in the first and the third 
quadrants. 

32. Find the integral surface passing through the 
given curve p 2 + q 2 = 4z where z = x 2 +1 and 
y = 0 then its general solution 

(a.) z = x 2 +(y + 1) 2 if b = a 
(b.) z = x 2 +(y- 1) 2 if b = -a 
(c.) z = x 2 +(y- 1) 2 if b = a 
(d.) z = x 2 + {y + 1) 2 if b = -a 

33. The particular integral of the differential 

equation (£> 3 - D)y = e x + e~ x , D s — is 

dx 


(a.) -^(e x +e~ x ) 


(b.) jx(e x + e- x ) 


(c.) jx 2 (e x + e' x ) 
(d.) \x\e x -e- x ) 


Assignment! 


Solution of 2«:^2ft =;si , , 

dx 2 dxdy dy 2 V ' 


(a.) z = <j> x (y+2x} + <j> 2 (2y + x) —cos(2x +y) 

3 


(b.) z = (j) x (y + 2x) + ^ 2 (2y+ x)--xsin(2x + y) 


(c.) z = (f\ (y + 2x) + ^ 2 (2y + 3x)--xsin(2x + y) 


H (&>) z = <j> x (y + 2x) + (j) 2 (2y + x) --xcos(2x + y) 

The ...general solution of 
'• (p + D'-2) (D'-P')z = 0 is/are 

Ja.) z = e 2x ^(x-y)+^ 2 (x + y) 

(b.) z=e 2x fa(x + y) + <t> 2 (x-y) 

(c.) z = e 2y i/> l (x-y) + <t> 2 (x + y) 

(d.) z = e 2y <j> x (x + y) + <t> 2 (x - y) 

In the (x,r) plane, the characteristics of the 
initial value problem u t +uu x = 0 , with 
u(;c,0) = x,0<x< 1, are 

(a.) Parallel straight lines 

(b.) Straight lines which intersect at (0,-1) 

(c.) Non-intersecting parabolas 

(d.) Concentric circles with centre at the 
origin 

Which of the following is NOT elliptic? 

(a.) Laplace equation 
(b.) u ix +2u xy +4u yy =0 

(c.) Harmonic function 
(d.) None of these 

The Partial differential equation 
y iu xx-(x 2 - 1)^=0 is 

(a.) Hyperbolic in {(x,y): y < 0,|x| > l} 

(b.) Parabolic in {{x,y ): y > 0,U| > l} 

(c.) Hyperbolic in {(x,y ): y > 0,|x| > l} 

(d.) Hyperbolic for all pair ( x,y) 
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39. Characteristic curves of y 2 r-x 2 t = 0 are 
(a.) x 2 +y=c l ,x 2 -/ =c 2 

(b.) x 1 + y 2 =CpX 3 —y 2 =c 2 
(c.) x 2 +y 2 =c„-x 2 +y 2 -y = c 2 
(d.) None of the above 

40. The characteristics curves of the equation 
X l u xx - y 2 u yy = x 2 y 2 + x; x > 0, u = u (x, y) are 

(a.) Rectangular hyperbola 
(b.) Parabola 
(c.) Circle 
(d.) Straight line 


I 
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The partial differential - equation 
xU a +2xyV xy +yU yy -xU y -y\J x = 0 is/are 

(a.) Hyperbolic in region xy< 1 

(b.) Hyperbolic in region xy> 1 

(c.) Elliptic in region xy < 1 

(d.) Elliptic in region xy> 1 


A general solution of the equation 

du(x,y) , . _ x . 

—— —+u(x,y) = e is 
dx 


\W 


x f{y) 

xb )'« (*> y)'-< x f{y) + xe x 

\ 

Nc.) u (x, y) = e x fiy) + xe~ x 
(d.) u(x,y) = e~ x f{y) + xe~ 
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Given 4 points A, B, C and D (in that order) 
along a parallelogram whose sides have 

slopes ±- in the (x, t) plane and u is a 
c 

solution of the wave equation u n -c 2 u xx = 0. 
Suppose u(A) = ^,u(C) = ^,u(B) = j. Then 
u(D) equals 


(b.)- 

12 


«L)1 

12 


. .. / \ —j=e 4 ‘ t>0,xe E\ 

The function u (x, t) = < y jt \ -- \ 

0 t<0,xs'R, 

is a solution of the heat equation in \ '"''x \ 

\ 

(a.) {(x, ():xeR, t e R} \^w 

(b.) {(x,t):xsR, t > 0} but not in the set 
{(x, t): x e R, t < O} 

(c.) {(x,t):xeR, teR}\{(0,0)} 

(d.) j(x,r):xeR,reR} 


4. The diffusion equation 


d 2 u _ dm 


u = u(x, t ), u( 0,0 = 0 = u(K, t), 

u(x, 0) = cos x sin 5x admits the solution 


p-y or _ 

(a.)-^sn^x+e 20 * sin4x 


-ibt 

(b.)-|^sin4x+e 20 'sin6xj 


V \x '-201 

\ •. / .V < 


>(c.) sin3x + e l5( sin 5x1 

\ — 2 V J 


Consider the boundary value problem 

U XK +Uyy=0,xe (o.njje (0,it), 

«(x,0) = w(x,it) = u(0,y) = 0 If an additional 
boundary condition « JC (7i,y) = siny is 
satisfied, then u(x,n/2 ) is equal to 

(a.) |(««-*")(<,' +«-) 

n(e n +e' x ) 

W M 

( e x- e -x) 

^ (e*+<r n ) 

(d.) 


{d.) ---j^sin5jc-Ke 20/ sirijcJ 

\ \ \ 

5. It is required to solve the Laplace equation 

d 2 u d 2 u n n n ' ■ . - . 

— r- + —- = 0, 0<x<a,0<y <b, satisfying 

- Sx 2 dy 

the boundary conditions 

u(x,0) = 0, «(x,f>) = 0, «(0,y) = 0 and 

u(a,y) = /(y). If c n ’s are constants, then the 

equation and the homogeneous boundary 
conditions determine the fundamental set of 
solutions of the form 


, \ , \ v • . nKX • mny 

(a.) u[x,y) = 2^c n sinU——sin—=- 

«=l b b 


si \ / , t-. . nnx . tiny 

(b.) u[x,y) = 2_jC n sin——sin—— 

n= i b b 


, s / v v . nnx . , nny 
(c.) u(x,y)= 2_c n sm—smh-f- 

«. l b b 


(d.) u(x,y) = Y i c„smh^-siah^- 

n= 1 b b 


Let u(x,t) be the bounded of—~^-x- = 0 
V ' . . dt dx 2 

e 2x -l 

with n(x,0) = ———.Then lim u(l,t) equals 
e x +1 

(a.)-1/2 
(b.) 1/2 
(c.)-l 


mmMMmmmmmm 
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7. The vertical displacement «(x, i) of an 

infinitely long elastic string is governed by 

the initial value problem 

d 2 u . d 2 u ' / 

—r -4 —-oo< a: <oo, t>0, u(x, 0) = — 

dt 2 dx 2 

and—(x, 0) = 0 The value of u (x, t) at x = 2 
dt 

and t = 2 is equal to 
(a.) 2 

(b-)4 

(c.)-2 

(d)-4 

8. Let u(x,t) be the solution of u lt =u xx ; 

0<x<l,f >0 ,m(x,0) = x(1-x), 

(\ n 

u, (x, 0) = 0. Then u is 

v2 4 ) 

(a.) — 

16 

(b> 4 
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11. For the diffusion problem 

= u/(0<x<n,t >0) , u(0,l) = 0 

«(jt,f) = 0 and w(x,0) = 3sin2x, the solution 

is given by 

(a.) 3e~‘ sin 2x 

(b.) 3c -4 ' sin2x 

(c.) 3e” 9, sin2x 

(d.) 3e~ 2t sin 2x 

12. Let PQRS be a rectangle in the first quadrant 
whose adjacent sides PQ and QR have slopes 
1 and -l respectively. If u{x,t) is a solution 

\^*Kr° and 

Vv, \nf-Rj = 1/2, then u(S) equals 

i v 

A ^>.)2 

N>s ~'< b ) 1 
\"\ C'-(c.) 1/2 


\ N\ _ u . d 2 u d 2 u 

\sjl3.~ Consider the wave equation — = 4— 

^ dx 1 

0<jc<7i,t>0, with w(0, t) = u(nj) = 0 , 

* 

x u (jc, 0) = sin x and — = 0 at t = 0 . Then 

^ dt 
* ( K 7C A . 

w U IS 

U 2 ) 

(a.) 2 
(b.)l 
( c ) 0 
(d-)-l 

14. Solve the heat equation 
u t -u xx =0,0<x<n,0<t <o ° 

u (x, 0) = sin x, 0 < x < it 

w(0, t ) = n(jt, t) = 0, t > 0 Then u 

(a.) Is unbounded in (0, rc)x(0, oo) 

(b.) Takes both positive and negative values 
in (0, 7t)x(0, oo) 

(c.) Is negative in (0, ti)x(0, oo) 

(d .) u(x,t)<e~‘ , for all 

(x, t)e(0, 7t)x(0, oo) 

15. Let u(x,t) be the solution of the initial value 


(d.) —1/2 


, \ ' 

9. Consider the Neumann problem.. 

u xx + u yy ~ 0 < x < n, -1 < y < 1 

u x (0, y) - u x (jr, y) = 0. The problem admits 
solution for 
(a.) a = 0, p = l 

(b.)a = -l,p = | 

(c.)a = l,P = -| 

(d.) a = 1, p = -n 

10. Consider the one dimensional heat equation 

du d 2 u . ■ . ., / ... i 

— = —— ,0<x< 1, t> 0 with the initial 
dt dx - 2 

condition «(jc, 0) = 2 cos 2 tcc and the boundary 

conditions — (0,/) = 0 = — (1,/) . The 

dr } dt y } 

temperature u(x,t) is given by 
(a.) u{x,t) = [-e~ 4n * cosItvc 
(b.) u(x 1 t) = l + e~ 4n / cos2tcc 
(c.) u(x,t) = \-e~ 4n * sin2nx 
(d.) u(x,t) = \ + e~ 4n 1 sin27Dc 


equation 


»i d 2 u d 2 u _ , . 

problem —--- = 0; w(jc,0) = sinjc; 

dt dx 

-1 • Then u(k,kI 2) equals 

(a.) ji/ 2 
(b.) 1/2 
(c.)-l 

(d.)l 
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Let: u(x, t) be the solution of. the one 20." 
dimensional wave equation 

u„ -41^ = 0, -00 < X < 00,/> o 

“M) = j 16 *’ [ X ^ 4, and 

[ 0, otherwise , 

.. (. \X\^2, 


u t (x,0) = 


0, orten** F ° f 1 < 1 < 3 > ”(2, 


( a -) j 16- ( 2_2 0 2 ] + '^[ 1 -min{l,;-i}] 

(b.) I[32-(2-20 2 -(2 + 20 2 ] +; 

( c .)I[ 32 _(2-2;) 2 -(2 + 2 f ) 2 ] + i 

(d.) - 16—(2—2f) j+-[l-max{l-f,-l}] 

Let u(x, t) be the solution to s the~wave._. 
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The solution Of the’ Cauchy problem 

Uyy(x,y)- u x X (x,y) = Oi 

«(x,0) = 0,iq,(x,0)-x is u(x,y) = 

(a •)- 

y 

(b.) x.y 
(c.) xy+- 


■\ \ \. 

2f;N .^ I£>{jf,.0 is the D’Alembert’s solution to the 

\\ \“ .V.. d 2 U d 2 U m _ n ...;a 


vW 




equation 


3 2 u, . 5 2 m, K '' s ' 

&2 (^0=-r(^0. \ 


\ r)xLij Cr-U . i 

wave equation,-= ——, xe R, / > 0, with 

^ dt 2 dx 2 

the condition «(x,0) = 0 and -^(x,0) = cosx, 


n \ 

\\ then u 0, — I is 


du \ «v \~—- 

h(x, 0) = cos(5tuc), — (x, 0) = 0 . Then, the V v 22. The function u(r, 0) satisfying the Laplace 

Ot \ x w 


value of u (1,1) is__ 

A function u(x,<) satisfies the wave equatioif\''-~' 


d 2 u _d 2 u 
dt 2 dx 1 


,0<x<l,t>0 


1 ^ 1 (. 1 


If “[2 ,0 J"4 , “i 1 ’2j“ 1 and “(°’2 J = 2 


8 2 u 1 du 1 d 2 u A ■ 2 

equation — 7 +- — +-^—- = 0, e<t <a 

Sr 2 r dr r 2 50 2 

subject to the conditions 
u{e, 0) = 1, «(e 2 ,0) = O is 


(a.) In 


then u 1 is 

12 J 




(b.)ln 


(c.) In 

r 


(°>7 


Zfr-e 2 


Let u(x,y) be a solution of Laplace’s 
equation on x 2 + y 2 <1 If 

oinAi—/ SdlG ^ <> S 0 £ « 


, „ . Sins lor U < 0 < IT i/uimuu.; - 

u(cos0,sin0) = l “ Then 

[0 for it < 0 < 2 ji by 

w(0,0) equals 0 

V ; 4 (a.) 1 + r 2 cos0 

(a.) 1 /ti 

( b ) 2 In (b.) l+lnr+rcos20 

(c.) 1/(2 ji) _ (c.) 2r 3 cos 2 0 

(d.) tc/ 2 (d.) l-r 2 + 2r 2 cos 2 0 ;r,f * ! * 
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(d.) £ - 2 sin " 0 

n-1\&~€ j 

23. The solution of Laplace’s equation 

XiX-L^OIn the unit disc with 
dr* r dr r * a) 2 

boundary conditions «(1,0) = 2cos 2 0 is given 

by 

(a.) l + r 2 cos0 
(b.) l+lnr+rcos20 
(c.) 2r 3 cos 2 0 


(d.) 1-r 2 +2r 2 cos 2 0 
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24. Consider the boundary value problem: 
u a +u yy = 0 in Q = |(x,y);x 2 +y 2 < lj with 

— = x 2 +y 2 on the boundary of Q ( — 
dn dn 

denotes the normal derivative of u). Then its 

solution u(x,y) 

(a.) Is unique and is identically zero 
(b.) Is unique up to a constant 
(c.) Does not exist 
(d.) Is Unique and non-zero 

25. Let u be a solution of the initial value 
problem 

-r -—7 = 0; w(x,0) = x 2 ,^-(x,0) = 0 .Then 
dt 2 dx 2 v & 

m( 0,1) equals , _ 


29. Let G(jc, y) be the Green’s function of the 
boundary value problem 

[(l+jc)ir']' + (sinx)« = 0, xe[0,l]u(0) = w(l) = 0 
Then the function g defined by 

g(x) = G^x,^j,x€[0,l] 

(a.) is continuous. 

(b.) is discontinuous at x= ~. 

(c.) is differentiable. 

\ ^Xd^does not have the left derivative at x =-. 


30. x Cdnsider ''the boundary value problem 

C\. ' u 'xx +u yy = 0,xe'('0,Tc),y e(0,n), 

''m (x, 0) = u (x, ji) = u (0, y) = 0 Any solution of 

\ \ 

N this boundary value problem is of the form 


\ X 


i N. \ 

26. Let u(x,t), xeR,t>0, be the solutioh,of the. 

initial value problem u tt = u a , «(x,0)^-,.. 
m,(x, 0) = 1. Then w(2,2) is equal to 

27. The solution of the initial value problem 

u„ = 4« n , t > 0, -oo < x < oo satisfying the 

conditions m(x,0) = x, u,(x, 0) = 0 is 


(a.) sinhnxsinny 


(b.) cosh nx sin ny 


(c.) ^a n sinhnxcoshny 


(d.) £a„cosh nx cos ny 

nel 

31. Let u(x,t) = e mx v(t) with v(0) = 1 be a solution 

du <5 3 k . 
to — = —rthen 
dt dx 3 


(a.) u(x,t) = e‘ 


_ ICO(X-G) /) 


(b.) u(x,t) = e K 


(c.)2x 
(d.) 2t 

28. Which of the following, concerning the 
solution of the Neumann problem for 
Laplace’s equation on a smooth bounded 
domain, is true? 

(a.) Solution is unique 

(b.) Solution is unique up to an additive 
constant 

(c.) Solution is unique up to a multiplicative 
constant 

(d.)No conclusion can be drawn about 
uniqueness 


(c.) u(x,t) = e‘ 
(d.) u(x,t) = e‘ 


_ Jo>(x+m t) 


/Gf(x-t) 


_ . , . , . du d u 

The partial differential equation—= — r+w 

dt dx 2 

, dv d 2 v 

can be transformed to — = —t For 
dt dx 2 

(a.) v = e~‘u. 

(b.) v = e'u. 

(c.) v = tu . 

(d.) v = -tu 
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34. 


35. 


36. 


Which of the following satisfies the heat 
equation without source terms and with 
diffusion constant in one space dimension? 


(a.) sin 


v 4 A 


(b.) e 1 sin x 

(c.) x 2 -t 

_-x 2 /4< 

(d.) 


i 

Let u(x,t) be the 
u tl -u JCX =\,xeR,t> 0, with 

u (.t, 0 ) = 0 , u t (x,0) = 0,xeR. 

Then u ( 1 / 2 , 1 / 2 ) 

(a.) - 
8 


(b.) -- 
8 


solution of 


(C >4 


If n(x,i) satisfies the wave equation 


< 4 > -i 



37. ....... Let a, b, c be continuous .functions defined 

on R 2 . Let V x , V 2 , V 3 be nonempty subset of 
R 2 such that V x \JV 2 (jV 3 = R 2 and the PDE 
«(*> y) u xx +b(x, +c(x, y)u n , =0 is 

elliptic in V x , parabolic in V 2 and hyperbolic 
in V 3 , then 

(a.) V x , V 2 and V 3 are open sets in R 2 
(b.) V x , and V 3 are open sets in R 2 
(c.) V x , and V 2 are open sets in R 2 
, (d.) V 2 , and V 3 are open sets in R 2 

38. '"'•Consider the heat equation 

it, t> 0 , with the boundary 
' A conditionsw( 0 , t)'= 0 ,«(n, t) = 0 for t > 0 , and the 

initial condition u(x, 0 ) = sin x. Then j is 

Let u(x,t) be the solution of the initial 


boundary value 

d 2 u d 2 u . 

—7 =— 7 , 0<x< od, r>0 


problem 


dx L 


(kx\ 

«(x, 0 ) = cos —, 0 < 
2 j 


X<oo 


d 2 u 2 3 2 « „ n 

~T = c — T-,xeR,t> 0 


df 


dx L 


sin- 


roc 


with initial 

0<x<c , 

, and 


^-(x, 0 ) = 0 , r >0 Then 

(a.) The value of u(2, 2 ) = -1 
(b.) The value of u(2, 2 ) = 1 


conditions u(x,0) = 

0 Elsewhere 
u, (x, 0 ) = 0 for all x, then for a given t > 0 , 

(a.) There are values of x at which u(x,t) is 
discontinuous 

(b.) u(x,t) is continuous but u x (x,t ) is not 
continuous 

(c.) u(x,t),u x (x,t) are continuous, but 
(x,t) is not continuous 

(d.) u(x,t) is smooth for all x 

The second order PDE - yu^ + x 2 y = 0 is 

(a.) Elliptic for all x e R, y e R 

(b.) Parabolic for all x e R, y e R 

(c.) Elliptic for all x e R, y < 0 

(d.) Hyperbolic for all x e R, y < 0 


(c.) The value of u 


I'l 


\l2 


(d.) The value of \j = \ 

40. Consider the Laplace equation in polar form: 

* 0; o< rS 4,0se<2« 
Sr 2 rdr r 2 30 2 

satisfying u(a, 0 ) = /( 0 ), where / is a given 

function. Let o be the separation constant 
that appears when one uses the method of 
separation of variables. Then for solution 
u(r, 0 ) to be bounded and also periodic in 0 
with period 2n, 

(a.) a cannot be negative 

(b.) a Can be zero, and in that case the 
solution is a constant 

(c.) a can be positive, and in the case it must 
be an integer 

(d.)The fundamental set of solutions is 
{l, r n sinnO, r n cos/i©} , where n is a 
positive integer. 
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The solution of the differential equation 44. Let the heat 

. ■ - ■ du d 2 u d 2 u dV . „ . 


equation 


—j = fix); x e(0, 1) y(0) = y(l) = 0is given 
dx - " 

l 

by yix) = J G{x, where 


(a.) G(x, %)= 


*U-i); x<t 

^(x-l); x>L 


S 2 (*-D; *>$ 


(c.) G{x,t,)=; 


(d.) G(x,J;) = 


x{?-\); X<\ 

%{x 2 -l); x>% 

sinjt(^-l); x<^ 
sin^(jc-l); x>i; \ 


admit an exponential function 
exp(i(£-x + wt)) as its solution, where k is a 
non zero constant real vector, and w is 
constant. Then the solution 

(a.) Remain constants on certain planes in 


(b.) Repeats itself after a certain length L 

K.^ (cQ Has, in general, an amplitude decaying 
V " ^ 'exponentially with time t 

' (dr) Is bounded uniformly for x e R 3 for a 
V< fixed t \ 

V \ 

N '4§f Vs "Let u(x,t) satisfy for xeIR,i>0 




Let u be a solution of the heat equation \ 
u, — u rY — 0 0 < x < n and (>0 
u(0,t)-u(n,t) = 0 t>0 > Then \ 

«(x,0) = sinx + sin2x, 0<x<7tl \ 

' 

(a.) u(x,t) —> 0 as t ->» for all x e (0 ,tQ n 

(b.) t 2 u{x,t) -> 0 as t -> oo for all x €.(0, nj v '" * 

(c.) e 2 u{x,t) is a bounded function for 
xe(0,7i),t>0 

(d.) e 2l u(x,t) -> 0 as t -> oo for all x e (0, n) 

Let u(x,t) be the solution of the equation 




d 2 u du „ d 2 u . , . „ . . 

—-+—+2—- = 0. A solution of the form 
dt 2 dt dx 2 

u = e a v(t) with v(0) = 0 and v'(0) = 1 


d 2 u du 


which tends to zero as t-> oo and 


a TT 111VU vvllviU VV £4V1V Ui7 l r W U11U 

dx 2 dt' 

has the value cos (x) when t = 0 then 

(a.) « = where a nA 

are arbitrary constants 

(b ) u = Xr=i a « sin ( /u+ ^ ) e_ " 2 ' where a nA 

are non-zero constants 

(c.) « = ^* =1 a n cos(nx+d„)e" n< where a n are 


not all zero and b„=0 for n > i 


(d.) u = Yf tl= \ a n cos(nx + d„)e " ( where 
a, # 0 , a n =0 for n> 1 , and b n =0 for 


(a.) Is necessary Bounded 
(b.) Satisfies \u(x,t)\<e‘ 

(c.) Is necessarily unbounded 
(d.) Is oscillatory in x 

46. The Green’s function G(x, Q, 0 < x, £ < 1 of 
the boundary value problem 
y" + Xy = 0, y(0) = 0 = y(l) is 

(a.) Symmetric in x and C, 

(b.) Continuous at x = C, 

01 *(*.01 _ , 


«/(*, ol <*(*.01 


47. ThePDE 

u xx +u yy +Xu = 0, 0 < x, y < 1 

«(x,0) = «(x, 1) = 0, 0 <x<l * has 

“(0, y) = «(1, y) = 0, 0 < y < 1 

(a.) A unique solution u for any X e R 

(b.)Infinitely many solutions for some kR 

(c.) A solution for countably many values of 
X v 

(d.) Infinitely many solutions fQr all X e R. 


8A/11, (First Floor) jia Sarai, Haiiz Khas;TSear LCT., New DelhM10016, Ph.: (dfl)-26537527, Cel!:'9999183434 & 9899161734,858*844789 

<__WeBstesxtSfKdipsati^smiiMiiiil’ ((iS; '_____— 



















REGtNERATI NG LOGICS 


An ISO 9001 : 2008 Certified Institute 


Assignment Sheet -3 


48. A bounded solution to the partial differential 

du d 2 u . 
equation — = —-+e is 
H dr dx 2 

(a.) u(x,t) = -e~* 


(b.) u(x, t) = e x e ‘ 


(c.) u(x,t) = e x +e‘ 

(d.) u(x,t) = x-e~‘ 

49. If the initial value problem for partial 

,. cc . . . du d 2 u 

differential equation —+—- = 0; 

dt dx 2 

u(x, 0) = sin(rct), has a solution of the form 
u[x, t) = 4>(r)sin(7ac),then 

(a.) <j) is always negative. 


(a.) a >0and6>0 
(b.) a>0and6 = 0 
(c.) a = 0and6 = 0 
(d.) a <0and2> = 0 

Consider the boundary value problem (BVP) 
»'' = -/,«(0) = t/'(l) = 0 on [0, l] , Where 


, du u d u 
i =— « s — 7 
.. dx dx 2 


u = —— u" s —- . Assume fix) is real- 


(b.) <(> is always positive. \ 

\ A' 

(c.) <|> is an increasing function. 

\ x 

(d.) (j> is a decreasing function. - . \—- 

50. Let u be of the boundary value problem 


Define for 


u"+-u' = mte{ 0,1) T 

t > l 

u\Q) = a,u(\) = b 

x 2 +y 2 <1, v(x,y) = u{jx 2 + y 2 


g(*,y) = f N x2 +y 2 h then v is a solution of 


v X x +v yy=g in \{x,y)'x 2 +y 2 
the PDE 1 | if 

v(x,y) = 0 on{(x,y);x 2 +y 2 =1^ 


Vqittied continued function on [0, 1]. Then, 
which of the fpllowing are correct? 
fa.) The Green’s function 

G(x,C,),(x,Q€[0, l]x[0, l] , for the 
above BVP is 

„ „ [x forO<x<C 

(b.) G(x,( 3 ) = \ Both G and 

[(, for C, < x < 1 

are continuous on [0, l]x[0, l] with 


—7 having a discontinuity along x = C, 
dx 

(c.) G(x, Q satisfies the homogenous 
equation u" = 0 for 0<x<(, and 
(,<x<\ 

(d.) The solution of the given BVP is 

"W=J 0 V(;)-';+jy (;)«<?• 
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